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EXPLANATION OF REVISIONS TO ASAP 



Since the issuance of the original ASAP report in December 
1974, and distribution of the source program, feedback from many 
users and direct assistance from Dr. Robert Bevensee of Lawrence 
Livermore Labs, University of California enabled the preparation 
of this revised version. During the past 2 years, no additional 
discrepancies have arisen and it is now felt that the code can be 
safely called "revised". 

The nature of the improvements was as follows: 

1. To correct the manner in which structures were 
treated when elevated over a ground plane. 

2. Improvements and corrections in calculating finite 
ground plane effects. 

3. Text corrections in explanation of ground effects. 

4. Instructional changes in the DESCRIPTION, FREQUENCY, 
and CHANGE cards. 

5. Corrected sample problem outputs. 

The user should note that some facilities will experience 
printing errors when NEAR FIELDS are called for. Statement 58 in 
the MAIN program should be suitably rewritten if that occurs. 
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!• INTRODUCTION 



Although many thin-wire computer programs have been 
developed for the purpose of analyzing antennas and 
scatterers, few of these programs have been directed toward 
the student of electro-magnetic theory. The majority of the 
programs are directed to the engineer or advanced student 
for the purpose of analyzing designed structures or 
verifying experimental data. 

The purpose of the study is to develop a computer 
program by modifying an existing computer code which can be 
utilized as an educational method to develop insight into 
radiating structures by the beginning student of 
electro-magnetic theory. 

The modified Ohio State University Antennas-Scat terers 
Analysis Program (OSUMOD or ASAP) is directed toward the 
beginning student who does not yet have the expertise 
necessary to manipulate the input data for proper execution 
of the larger more comprehensive analysis program. Even 
though ASAP is small in core requirements and is fast in run 
time, it is capable of analyzing structures to assist the 
engineer with design problems. 

Since the resulting program, ASAP, is primarily directed 
toward students, the program has been limited to structures 
which contain less than 50 monopoles (segments) , no longer 
than one-fourth of a wavelength, and which have less than 51 
nodes (intersections and endpoints) . If a ground plane, 
either perfect or finite is present; the stated limits above 
are halved due to the generation of an image structure. 
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II« ORIGINAL PROGRAM 



A. THEORY 

Reference 1 presents the electro-magnetic theory for the 
analysis of antennas and scatterers in an isotropic, linear, 
and homogeneous ambient medium. The analysis is performed in 
the frequency domain with an excitation caused by either a 
generator or an incident wave. 

In the analysis, a piecewise-sinusoidal expansion is 
used for the current distribution. The matrix equation 
Z I = V is generated by enforcing reaction tests with a set 
of sinusoidal dipoles located in the interior region of the 
wire. Since the current distribution has the same form as 
the expansion mode, this formulation is known as the 
"sinusoidal reaction technique". 

B. COMPUTER PROGRAM 

Reference 2 presents the computer program corresponding 
to the theory presented in Ref. 1. 

1. Input Format 

In the program, the input data must specify the 
frequency, wire radius, wire conductivity, the parameters of 
the exterior medium, coordinates of the points to describe 
the shape and size of the wire configuration, a list of the 
wire segments, and the indicators for the various outputs. 
Table 1 is the input data necessary to analyze a half-wave 
dipole. 



2. Output Format 

In the original form, the only outputs which could 
be requested by the input data stream are the following: 
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a. Antenna Problems 





(D 


Current Distribution on the Structure. 




(2) 


Input Impedance. 




(3) 


Radiation Efficiency. 




(4) 


Near-Zone Field. 




(5) 


Far-Zone Field. 


b. 


Backscatter ing Problems 




(1) 


Absorption Cross Section. 




(2) 


Scattering Cross Section. 

• 




(3) 


Extinction Cross section. 


Scattering 


(4) Complex Elements of the Polarization 
Matrix 


c. 


Bistatic Scattering Problems 




Echo 


Area. 



Table 2 is an example of the output data available for 
data of table 1. 

3. LIMITATION 

Although the program can analyze a structure with up to 
50 segments, 55 points and 60 dipoles modes; it can not 
analyze a structure in the presence of a finite ground 
plane. 
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AN EXAMPLE OF THE OUTPUT DATA FOR THE ORIGINAL PROGRAM 



Ill* MODIFIED COMPUTER PROGRAM 



A. Input Format 

As illustrated in table 1 the format for the input data 
cards is not self explanatory. This format can be determined 
by referring to the FORMAT statements of the program of Ref. 
2. Since the modified program is directed toward the 
student, the input data format was changed to allow free 
format. Reference 2 was written in a form which permitted 
modifications to allow flexibility in specifying input data 
for the analysis program. Appendix B, titled "User's 
Manual", discusses the input data cards necessary for proper 
execution of an analysis problem. Appendix B is self- 
contained and may be used independently of the remainder of 
this document. 

B. Output Format 

In the original computer program, the absence of labels 
encumbered the output data and lessened the usefulness of 
the program. To improve the usefulness of the modified 
version, detailed labels were added to the output data. As 
with the input data. Ref. 2 was written in a form which 
enabled modification to allow more specific output data for 
the analyzed problem. With the addition of the polar 
plotting package, the far-zone electric field intensity 
polar radiation and reradiation patterns can be plotted. A 
sample problem can be found on page 120 in Appendix B, 

User's Manual. 

C. Finite Ground 

To enable the student or the engineer to have an 
improved analysis program, the finite ground effects were 
added to ASAP. The theory corresponding to the ground 
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effects, which utilize Fresnel reflection coefficients, is 
discussed in Appendix A, titled "System Manual". Also 
discussed in Appendix A is the modified computer program and 
the corresponding theory. The electro-magnetic theory was 
developed in Refs. 1, 2, and 3; and it is restated with its 
corresponding computer code to assist in the understanding 
of the methods applied. Appendix A is self-contained and 
may be used independently of the remainder of this document. 



IV* CONCLUSION 

The addition of ground effect techniques to the original 
program did not alter the accuracy or the computational 
capabilities of the program. The ground effect techniques 
utilized the results of the original program and modified 
these results to account for the effects of the presence of 
the finite ground. 

To verify the numerical results of ASAP, the input 
impedances of both a horizontal and a vertical dipole were 
compared to the solutions of the exact form of the 
Sommerf ield ' s equation. As can be seen in table 3 the 
finite ground treatment of ASAP agrees favorably with 
Sommerf ield ' s solutions. The ASAP finite ground results are 
also in excellent agreement with the previous computer 
solutions of Refs. 4 and 5. 



V* RECOMMENDATIONS 



Although the program is a general analysis tool for 
students, several future modifications will enhance the 
program as a design tool for engineers. These items include: 
varying the wire radius on the structure; incorporation of 
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non-radiating elements such as transmission lines; varying 
the wire insulation radius, conductivity, and dielectric 
constant; and a geometry generation package such as dipole 
array or helix. One major change that would both improve the 
speed and reduce the core requirement is that of symmetry. 

No attempt was made to utilize the symmetry in the 
admittance matrix when the ground plane is present. If 
symmetry were applied, the structure size limit with the 
ground plane present would be approximately that of the 
structure without the ground plane. 



VERTICAL 01 POLE 
FREQUENCY 3 MHZ 
LENGTH .5 WAVELENGTH 

RADIUS .005 METERS 

DIELECTRIC CONSTANT (RELATIVE) 
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CONDUCTIVITY HEIGHT/WAVELENGTH 

.1 'll 

.30 

.35 

.45 

.001 .25 

.30 
.35 
.45 

.00001 .25 

.30 
.35 
.45 



ASAP 



123. 75* J 
98.62* J 
87.60* J 
78 • 69* J 
107.18*J 
91. 80* J 
8 5 • 15* J 
0O.26*J 
103. 20* J 
91. 83* J 
85.90* J 
86 .78*J 



68.30 
38.26 
35.64 
41.79 
55.36 
40.18 

39.30 
43.52 

57.99 

41.99 
40.14 
43.23 



EXACT* 



126. 5*J 
100. 2 * J 
88 • 50* J 
79 • 2 l * J 
119. 4*J 
94 . 30* J 
85 • 4 1 * J 
80 • 1 7* J 
1 15. l*J 
94 .09* J 
86 .0 1* J 
80.72*J 



83.89 

49.52 

46.52 
52.60 
71.46 
49.02 
48.93 

54.83 

73.89 
50.69 

49.83 
34.63 



HORIZONTAL DIPOLE 
FREQUENCY 3 MHZ 
LENGTH .5 WAVELENGTH 

RADIUS .001 METERS 

DIELECTRIC CONSTANT (RELATIVE) 



10 



CONDUCTIVITY HEIGHT/WAVELENGTH 



.1 



.001 



.00001 



.5 
.3 
.1 
.5 
.3 
• 1 
.5 
.3 
.1 



ASAP 



84 • 20* J 
132. 45* J 
39. 46* J 
88. 09* J 
1 1 7 . 69* J 
71 .92* J 
90 . 37* J 
1 15. 7» *J 
70 . 07 * J 



23.69 

46.09 

84.99 

31.80 

42.39 

71.07 

31.33 

45.06 

63.07 



EXACT* 



87 • 74* J 41.77 
136. 0*J 67.91 
40 . 54* J 89.71 
91 • 41* J 
120. 3*J 
78 • 1 7*J 
93 *8 5*J 



50.65 

62.69 

69.91 



'iD * - _ - 

l£f.?*J 



6 7*J 



60.81 



* COURTESY OF LAWRENCE LIVERMORE LABORATORY 



TABLE 3 
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SYSTEM MANUAL 



INTRODUCTION: The Ant ennas-scatt er ers Analysis Program 

(ASAP) for thin wire structures in a homogenous conducting 
medium performs a frequency domain analysis of antennas and 
scatterers. The program is applicable in the presence of a 
ground either perfect or finite. This appendix will describe 
the computer program which accomplishes this. Although the 
program was written for the IBM 360 computer system it can 
be executed on another system with minor modifications. 

A piecewise-sinusoidal expansion is used for the current 
distribution. The matrix equation ZI = V is generated by 
enforcing reaction tests with a set of sinusoidal dipoles 
located in the interior region of the wire. Since the test 
dipoles have the same current distribution as the expansion 
modes, this may be regarded as an application of Galerkin’s 
method. Rumsey's reaction concept was most helpful in this 
development, and therefore the formualtion is known as the 
"sinusoidal reaction technique". 

The main routine and each subroutine is discussed 
separately in this appendix. The writeups for the 
subroutines are arranged alphabetically by subroutine name 
after the main program. Each of the discussions includes the 
purpose of the subroutine, brief description, and a listing. 
After the subroutine writeups is a table of the more common 
symbols used in this program. 

The input data and program limits are discussed in 
detail in the next appendix titled "USERS MANUAL". 
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GROUND EFFECTS: In the modified antenna analysis 

computer program finite and infinite ground effects were 
added by using the reflection coefficient technique. The 
method in which this technique was used required the 
generation of an image structure. In this section the 
reflection technique will be discussed in detail. 

In order to apply ground effects to the electric field, 
the field for the image structure was first calculated as if 
a ground were not present. Then, the field was decomposed 
into parallel and perpendicular components. (A parallel 
component is the component which is parallel to the plane of 
incidence. A perpendicular component is one which is 
perpendicular to this plane. The plane of incidence is the 
plane containing the normal to the reflecting surface and 
the incident ray.) 



Z 1 
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Consider an image monopole with the electric field in 
the 1 direction. The ray, 6, is a vector which is 
perpendicular to 1 and passes thru the point of interest. To 
apply reflection technique, the plane of incident must be 
found. It is advantageous to define a new coordinate system 
(a//3/Z) where a and are parallel to the xy plane with a in 
the plane of incident and (3 perpendicular. 



If the direction cosines (cosx, cosy, and cosz) are 
known, it can be shown that the components of the field in 
the a/3 (xy) plane have the following relationship: 








— 




— — 


1 1 










E 




cos 4> sin <j) 




E 










X 


1 










E 




sin <J> -cos <(> 




E 


___ 








L Y 



where $ = arctan (cosy/cosx) . 



Now the reflection coefficients for the interface can be 
applied as: 



I I (R) 

E = R E 

I I 



I I 



±(R) 



= R E 



where R and R will be defined later in this section. 

I I 1 

Applying the matrix equation above yeilds: 
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* (R) 

E 

X 




COS 


* 


sin <|> 




1 1 (R) 

E 


<R) 

E 

y 




sin 


4> 


-cos <f> 




1 (R) 

E 



(the square matrix is unique, in that, the inverse is equal 
to the original matrix) . Since the image direction is 
opposite to the original monopole, that is. 



(1 x z) . . = - (I x z) . , 

original image 

the 2 component of the field, which is in the plane of 
incident, is given by: 




From electro- magnetic theory the reflection coefficients 
for the fields in medium (1) at the interface with another 
medium (2) are defined as: 

for perpendicular 

cos 0 - y 0 ' - sin 2 0 

R = ... 

H cos 0 + y & ' - sin 2 0 

and for parallel 

S' cos 0 - y &' - sin 2 0 

R = ===== 

V &• cos 0 + y &• - sin 2 0 

where 0 is the angle of incident as measured from the normal 
to the interface and 



&' = (6 + o /jco) / (& + a /ioi) 

2 2 11 



where the subscripts correspond to the mediums above. 
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To determine the relationship between R , R and R , R 

I I _L V H 

a perfect ground (S = 0, a =oc ) was investigated. 

r 



limit R - - 1 
H 



limit R = +1 
V 



But, for a perfect ground the contributions to the field 
from the image monopole would be equal to the field of tne 
original monopole but opposite in sign due to the chosen 
reference direction. 



11 (H) 11 I I 

E = R E = - E 

i I 



R = - E 

JL 



therefore 



R = - R 

I 1 v 



R = R 
1 H 



In summary, the contribution to the electric field of a 
monopole over a ground plane at a given point is given by: 



(R) (R) (fi) (R) 

E = E cosx + E cosy + E cosz 
x y z 



where 



(B) 

E = R E COSX + (R - R ) E COSX COS 2 <t> 

X _L I I _L 



+ (R - R ) E cosy sin <t> cos <|> 

1 1 _L 
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E cosy - (R 



where 



E 



(R) 

y 



R 

i I 



i I 




E cosy cos 2 i> 



+ (R - R ) E cosx sin <|> cos <J> 

I I _L 



(R) 

E = - R E cosz 

z I I 

is the field without the ground plane present and 



R 

1 1 


- S' cos 0 + 


y &• • 


- sin 2 


S' cos 0 + 


y & » - 


sin 2 0 


R 

± 


cos 0 - y S* 


- sin 2 


0 


cos 0 + y S ' 


- sin 2 


0 



&• = & - j (a/S co) 

r 0 
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MAIN 



PURPOSE: to control the input, output, and the flow of 

calculations. 

METHOD: The main program controls the flow of the 

required calculations by calling only a few subroutines. 
These subroutines in turn call other subroutines which 
actually do the required calculations. The order of the 
calling sequence is diagramed after the listing for the main 
program. 

The DIMENSION statements at the beginning of the main 
routine provides the required storage for a wire structure 
with up to 50 segments, 60 nodes and 60 dipoles without the 
presence of a ground plane. If a ground plane is present 
one-half of the reserved storage is required for the image, 
therefore a wire structure with up to 25 segments and 30 
nodes can be analyzed. 

NM denotes the actual number of monopoles (segments) , 

INM is the corresponding dimension, and the dimension for 
CG, VG, and ZLD is twice INM. The second subscript for MD 
always has a dimension of 4 to correspond to the number of 
segments meeting at a given node. 

N denotes the number of simultaneous linear equations 
and ICJ is the corresponding dimension. The dimension for C 
is (ICJ * ICJ + ICJ)/2. 

In the statements above statement 4, the initial 
conditions and defaults are established. After calling 
subroutine READ to determine the input parameters, the IF 
statements output the parameters to be used for the 
calculations. In the DO LOOP ending at statement 7, the 
the input data of the structure geometry is stored in order 
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to recall if the structure is to be moved for ground plane 
calculations . 

After the image structure is generated and structure 
location is moved, subroutine SORT is called to determine 
the dipole modes. Prior to calling SGANT, the load and 
generator information is established. 

Subroutine SGANT is then called to calculate the 
elements of the impedance matrix. If FEEDS or GENERATORS 
are specified by the input data stream, subroutine GANT1 is 
called to solve for the current distribution due to these 
forcing functions. 

In the DO LOOP ending with statement 29, subroutine 
GNFLD is called to calculate the near-zone field for the 
current distribution of the subroutine G A NT 1 . 

The subroutine GFFLD is called for the far-zone field of 
the current distribution of the subroutine GANT 1 in the DO 
LOOP ending at statement 35. The subroutine GFFLD is called 
again in DO LOOPS ending at statements 42 and 51, if 
bistatic and backscattering calculations are requested by 
the input data stream. 

CALLS TO: GANT1 

GFFLD 

GNFLD 

POLPRT 

READ 

SGANT 

SORT 



24 



on 



0 I HE NS ! ON X( 60) « 
DIMENSION 11(60) 



ji, Y (60 I * Z(60), XG( 60 ) « YG(60), ZGI60I 

---30), 12(60)* 13(60), J A( 63 ) , J8(60), KFLAG4 30 ) 

DIMENSION CPHI (500) , C THE T ( 500) 



uincniiun ikmusuuj . 

DIMENSION 0ATYK 360I. 0ATY2(360), 0ATY3( 360) r 0ATY4(360) 
DIMENSION 0(50). I A( 50 ) t 10(50), 1SC(50), MD150,4), NOlSO) 
01 MENS I ON L ZO ( 60 ) , KGEN(60) 

COMMON I WL 

0 1 MENS I ON XNP ( 5 0), YNP(50), ZNP(50) 

;gmplex c( 1030) 

- - Tl($ 



COMPLEX COATK 500).CDAT2(500),CDAT3(500) . 
COMPLEX C J (60) , EP ( 60) • EPP ( 60 ) , E T ( 60) ,ETH 
COMPLEX CGO(50). SGO(60) ,CG( 1 05 ) , VG( 1 00 ) • l 
COMPLEX VOLT(60l .ZLLD(60T 

Complex eppSiEptS.etp$,etts,ex,ey.ez 
COMPLEX EP2,EP3,EP4,eAA,ETA,GAM,Yll,Zli,Z$ 
OATA PI ,TP/3.14l59a6. 28310/ 

OATA EQ,UO/0.854E-l2,1.2566E-6/ 

1 NGEN - -1 
I GAO * -1 
LOAO ■ -1 

0 M « *1 
1CAR0 » 0 
AM - *1 

1 FLAG « 0 

VOL TCI) - (l.,0.) 

HGT - 0. 

NM ■ 0 
NP « 0 
MSG * 0 
SIG2 « -1. 

T02 « -1. 

SIG3 - -I 



CDAT4I500) 
60) 
n 



ZLOUOO) 



ER3 

T03 

CMM 

ER2 

FMC 

INM 

1CJ 



i 
0. 
50. 

ioo. 

50 

60 



WRITE ( 6,741 

DO 2 1-1,30 
2 KFLAG(I) - -1 



DO 3 J* 1 , I NM 
I SCI J I • 0 
VG(J) - (.0,.0) 



0001 

0002 

0003 

0004 
0C05 
0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 
0027 
0020 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 
0030 

0039 

0040 

88i! 

0043 

0044 

0045 

0046 

0047 

0048 



ZLO( J ) - ( .0, .0) 
JJ - J* INM 
VG(JJ) - (.0..0) 
(. 6, .6 



3 ZLO(JJ) - 



NFFP 
N0IP 
N0AP 
AFFP 
AFFT 
A 0 I P 
ABIT 
A6AP 
A0 AT 
STEP 
KNM 



0 

0 

0 

1000 . 

1000 . 

1000 . 

1000. 

1000 . 

1000 . 

1. 



CALL RE AO ( IA. 10, 10 ISC, ICARO, IGA IN, 1 GAD, INEAA.INT, I SCA T. IWR , I FLAG , 
1KFLAG,KGEN,L0A0,LZ0,MSG, N0AP.N0 1 P, NF FP, NGEN ,NM. NP,A0 AP, AB AT , A FFP. A 
2FFT 1 A&IP 1 AbIT,AM.0M,CMM 1 EA2,tR3,ER4|FMC,HGT #PHAF , PHA I • PHI F , PHI I ,PH 
3SF , PHSI ,THAF.THAIfTHlF,THII • THSF , TH§ I , S f G2 * 5 1 G3, i I G4 , T02 , T03 , VOLT , 
4X, XNP, Y,YNP.Z,ZLLD,ZNP, STEP) 

WRITE (6,56) 

IF ( MSG.LT.l ) 30 TO 5 

IF (MSG.EQ.l) WRITE (6.70) KFLAGUO) 

IF (IFLAG.EQ.4) GO TO 1 

5 IF (IFLAG.EQ.5I STOP 

IF (AM.LT.O) WRITE 16,127) 

IF (AM.LT.O) GO TO 6 

IF ( (NM.GT.O) .ANO.(NP.GT.O) ) GO TO 7 
WRITE (6.116) 

6 IF (IFLAG.EQ.1) GO TO 1 
MSG - 2 

GO TO 4 

7 WRITE (6,114) 

WRITE (6,113) 

WRITE (6,1121 

F (KFLAG( ll.EQ.l 
(KFL AG( 2) • EQ. 1 
( KFL AG ( 3) • EQ • 1 
(KFLAG( 20) .NE 



F 
F 
F 
F 
F 

F ( KFL » G( 5 ) • EQ . 1 
F (KFLAG( 6) .EQ. 1 
F (KFLAG( 7) • EQ. 1 
F (KFLAGI0I.*E.1 
F (KFL AGj 9^|E9* 1 



I 



WRITE 
WRITE 
. WRITE 
1) WRITE 



(6,03 
(6,84 
( 6 , 05 j 
(6,07) 



FMC 

AM 

CMM 



F CKFLAGI 



WRITE (6,86) 
WRITE (6,00) 
WRITE ( 6 , 09 ‘ 
WRITE (6,90 
WRITE (6,91 
WRITE (6,92 
WRITE i^93 



EQ * i ) WAITE 




, $ l^ 



0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 
0050 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 
0067 
0060 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 
0070 

0079 

0080 
0001 
0002 
0003 
0084 

0005 

0006 
0087 
0080 
0009 

0090 

0091 

0092 

0093 

0094 

8r 



96 



25 



(KFLAG411l.EJ.il WRITE 
(KFLAG(26).NE.l) WRITE 



(6,95) 

( 6 , 122 ) 



T 03 



WRITE (6,1231 
WRITE (6,125) 
WRITE (6,1241 
WRITE (6,126) 



ER4,SIG4 

HGT 



l^r'w Abl | tnlc* 1 ) MK1 1C (Ollcc J 

m GRD.GT.l) .AND.(KFLAG( 25) .EQ.l) ) 

GR0.EQ.1) . AN0.(KFLAG( 25). EQ.l) ) 

(( IGRO.GT.l) . AND. (KFL AG ( 25 I .EC.l) I 
( ( I GRD • GT . 0 ) .ANO. ( KFL AG ( 25) . EQ. 1) ) 

(KFLAG( 21) .EQ.l ) WRITE (6,121) INT 
WRITE (6,111) 

IF (KFL AG( 12 ) .EQ • 1 ) WRITE (6.96) C 1 , K A ( I ) ,X ( I A( I ) ) , Y( I A( II),Z.(KA€ 1 
l)),I8(l),xU8(l)),Y(I8(I)l,Z(IB(l)),I*i,NM) 

WRITE (6.111) 

IF (KFLAG(24) .GT .0) WRITE (6,119) ( L Z0( I I , ZLLD( I ) , I- 1 , LO AO ) 

IF (KFL AG( 14) .GT .0) WRITE (6,118) (LZO( I ),ZLL0(I),I*1,L0A0) 

WRITE (6,111) 

IF (KFLAG(23).GT.O) WRITE (6,120) (KGEN( I ) , VOLT( I ) , I *1 , NGEN) 

IF ( KFL AG ( 1 3 ) • GT .01 WRITE (6,97) ( KGEN4 I ) , VOL T( I ) , f « 1 , NGEN) 

WRITE (6,111) 

WRITE (6,114) 

WRITE (6,98) 

WRITE (6,112) 

IF (KFLAG4 22I.NE .1) WRITE (6,110) 

IF (KFLAG(I5).EQ.l) WRITE (6,99) 

IF (KFLAG( 16I.EJ.1 j 



IF (KFLAGU7I.EQ. 



WRITE (6,100) 
WRITE (6,101) 



PHAI, PHAF, THAI, THAF, STEP 
PHI I.PHIF ,THI I,THIF*STEP 



IF (KFL AG ( 16I.E0.1) WRITE (6,1021 PHSI , PHSF ,THSl , TH$F. ST EP 
“ (KFLAG( 19). EQ.l) WRITE (6.103) (XNP( I I , YNP( 1 ) , ZNP( I ) ,1-1, I NEAR) 



IF _ 

F ( AFFP.LT.500. 

F (AFFT.LT. 500 
F (A8AP.LT. 500 
IF (A0AT.LT. 500 
IF (A8IP.LT .500. j 



ITE (6,105) AFFP 
WRITE (6,104) 



WR 



3SSTI If* 

WRITE 
WRITE 



109) 

6, 108) 

IF (A8iT.LT. 500.' I WRITE ItltoZI AO I I 
F (( 18 ISC. GT. 01 .ANO. 4ISCAT.LT. Oil WRITE (6,73) 
F (KFL AG(4) .LT. 1 ) GO TO 129 



AFFT 

A8AP 

ABAT 

A8IP 

A8IT 



00 128 1*1, INM 

128 ISC(I)*1 

129 FHZ * FMCM.E6 
OMEGA « TP*FHZ 

IF * * * 

IF 

IF ( S IG3.LT. 0. ) 
“ T03.LT ‘ 



(S IG2.LT. 0. I EP2*ER2*E0*CMPLX(l* 1 -T02) 
(T02.LT. 0.) EP2 * CMPLX(ER2*E0,-$IG2/0MEGA) 



If 

If 

a TA 



(iGRO.GTdi EP4 
(IGRO.GT.l) ‘ 
(KFL AG( 21 ) 



EP3« 
EP3 ■ 



ER3*E0*CMPLX 
CMPLX4 ER3*“ 
CMPLX ( ER4* 
ERR * EP4/6P3 
GT.O) WRITE (6,121) 



A RTIU0/EP3) 

AM *_0MgGA*CSQRT(-U0*EP3) 



S*-£lG3/&M|GA) 
0,-SIG4/OMEGA) 

INT 



0097 

0090 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 
0109 



811 



8 1 



0 

1 

2 

3 

4 

5 

6 
7 
6 

x 9 

20 

w 

4 

25 

26 
7 



IF (KFLAG(12).NE.l) GO TO 9 



29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

140 

141 

142 

143 
0144 



NPG 
NMG 

00 8 I 
XG( I) 
YG( I 
"“II 



NP 

NM 



8 ZG( I) 



1 1 , NPG 

i i!" 



ill! 



9 00 10 
X ( 1 I 
I) 



10 



Y(I 

Zd) 



1*1, NPG 
* XG( I) 

« YG( I) 

« Z G ( I ) 



NP 

NH 

WL 



NPG 
NMG 
0 

F (1GR0.LE.0 
ET UP IMAGE 
ZM IN * 2(1) 

K - 0 



) GO TO 15 
FOR GROUND PLANE 



IF (Z(I).LT.ZMIN) ZMIN*Z ( I ) 

00 11 1*1, NP 
Z(I) * Z(1)*HGT 

IF (Z( II.GT.l.E-60) GO TO ll 

1 WL * 1 WL + 1 

11 CONTINUE 

IF (ZMIN.GE.O.O) GOT 01 2 
WRITE (6,117) 

IF ( IFL AG.EQ.l ) GO TO 1 
IF (IFLAG.E0.2I STOP 
MSG - 2 
GO TO 4 

12 00 13 J *1 , NM 
K * J ♦ NM 

I A( K I * I A ( J ) 



13 



IF (IA(J).GT. IWL ) IA(K)*IA( JI4NP-IWL 
I B( K) « I8(J)*NPMWL 



IWLP 



I WL M 



00 14 I*IWLP,NP 
J * I ♦ NP- I WL 

x(j) « x(h 

Y(J) * Y(I) 



BHI 



26 



14 Z( J) « -Z( I) 



15 



16 



17 



18 

19 

20 

21 



-I WL 
( U 



i 12, 

16 



[ 3 , J A, J8 * MD , NO, NM , NP »N, MAX * MI N, ICJ, INN) 



GO TO 
STOP 



1 



GO TO 1 
STOP 



KNM « NMM 
NM * 2* N " 

NP * 2* NP- 

CALL SORT (lA,13,Ili 
IF (MAX.LE.4) GO TO 
WRITE ( 6.71) 

IF (IFLAG.EQ.1) 

IF (IFLAG.EQ.2) 

MSG * 2 
GO TO 4 

IF (MIN.GE.l) GO TO 17 
WRITE (6,72) 

IF (IFLAG.EO.l) 

IF (IFLAG.EQ.2) 

MSG * 2 
GO TO 4 
WRITE (6,56) 

IF ( MAX.GT.4.0R.MIN.LT.1.0R.N.GT.ICJ) GO TO 54 
112 - l 

IF ( LOAD. GT *0 ) GO TO 19 

DO 18 1*1, NM 
ZLO(I) * (O.,0.) 

IF (NGEN.GT.OI GO TO 21 

DO 20 I *1 . NM 

V G < I ) * (0.,0.) 

KN • NM 

IF ( I GRD .GT. 0 I KN * NM/2 
J * 1 

ANTENNA CALCULATIONS 
IF (LOAD.LE.O) GO TO 24 
IF (KFLAG(24).GT.O) GO TO 22 

00 23 J*l , KN 

22 00 23 1*1. LOAD 
LZD( I f 



23 



IF 


(( I A ( J ) • 


EO. 


If 


(KFL AG( 2 


4) . 


If 


( (KFLAGl 


14 


if 


« I KF lag 


24 


CONTINUE 





) .GT.O) ZLO( K) ■ ZLLD( I ) 

4>.GT.0). ANO. j IGRD.GT .0) ) ZLO(J»KN)*ZLO(J) 
4 1 *GT <, 0) * AND* ( I GRO * GT •Oil ZLO(K+KNi*ZLO(K) 



019 3 

0194 

0195 

0196 

0197 

0198 

0199 

0200 
0201 
0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 

o°ill 

0213 

0214 

0215 

0216 
021 ? 
0218 
02 19 
0220 
0221 
0222 

0223 

0224 

0225 

8 m 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 



24 IF (NG6N.LT. 0) GO TO 27 
KN • NM 

IF (IGRO.GT.O) KN * NM/2 
IF (KFLAG(23).GT.O) GO TO 25 

00 26 J* 1 , KN 



25 



00 26 I-l.NGEN 
K * KGEN( I ) 

IF (( I A ( J) .EQ.K) .AND. (KFLAG( 131. GT. 0)1 VG< J ) * VOLT ( I ) 
IF (KFLAGI 23) .GT.O) VG( K ) *VOLT ( I I 

. _ _ , G f .C 



IF ((KFLAGU3). GT.O). AND. ( I GRO . 



r.O) ) VG< J*KN)— VG( J) 

IF (( IGRO.GT.O) .ANO.(KFLAG( 23) .GT. oil VG(K*KN)*-VGiK) 
26 CONTINUE 



27 



28 



CALL SGANT ( I A , I 6 , INM , INT . I SC , 1 1 1 1 2, 1 3, JA, JB, MD, N. NOtNM.NP , AM , 8M ,C 
• CGDfCMM P D.EP2,EP3-fcrA,FHit GAM ,SGD,X,Y,Z,ZLD,ZS,ERR* IGRDi 
IF (N.GT.O) GO TO 28 
IF (IFLAG.EQ.2) STOP 
MSG * 2 

IF (IFLAG.EQ.l) GO TO 1 
GO TO 4 

IF (NGEN.LE.O) GO TO 36 
WRITE (6,75) 

WRITE (6,76) 

WRITE (6,77) 

WRITE (6 t 82) 

^ ‘ 8, MD. N, ND,NM,AM,C,CJ,CG 



29 



30 



CALL GANfl ( IA,I8,INM,IWR, II, I2 t I3,;i2,JA.J8 t MD.N, 

L , CMM, D.EFF , GAM.GG.CGO ,SGO ,VG,Yllt21l«ZLD , ZS , IGRD) 

WRITE (6,57) E FF f GG, Z 1 1 
NEAR FIELD 
IF ( INEAR. LE.O) GO TO 30 
WRITE (6,75) 

WRITE (6,78) 

WRITE (6,77) 

DO 29 1*1. INEAR 
XP * XNP( 1 ) 

YP « YNP( I ) 

ZP * ZNPi I ) 

CALL GNFLD ( I A , I 8 , I NM , 1 1 . 1 2 . 1 3.M0, N, NO, NN, A M, CGO, SGO ,6 T A , GAM, C J , D , 
IX,Y,Z,XP,YP.ZP.EX,EY.EZ, iGRO.EftRT 
WRITE (6,56) XP,YP,ZP 
WRITE (6,59) EX,EY,EZ 
CONTINUE 

FAR FIELO 

IF ( IGA IN.LE • 0) GO TO 36 



0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 
02 59 
0260 
0261 
0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 
02 72 

0273 

0274 

0275 

0276 

0277 

0278 

mi 

0281 

0282 

0283 

0284 

0285 

811* 

0286 



27 



00 31 1*1, 
DATY1 ( 1 ) * 
0ATY2(I) • 
0ATY3C I ) • 
31 DATY4(I) « 



360 

0 

0 

0 

0 



WRITE (6,75) 

WRITE (6,79) 

WRITE (6,77) 

WRITE (6,02) 

INC * 0 
NPL * -l 

IF (KFLAG( 16) .EQ.l J WRITE (6,69) 
IF (NFFP.E0.1) GO TO 32 
NPHA - ( PHAF-PHAI )/STEP+l 
NTHA « (THAF-THAtWSTEPM 
GO TO 34 

32 IF (AFFT.GT.500.) GO TO 33 
NPL * L 

NPHA - 360 
NTHA ■ 1 
PHAI - 0. 

THAI - AFFT 

step - l. 

GO TO 34 

33 NPL - 2 

1 

360 
AFFP 
0. 

1. 



NPHA 
NTHA 
PHAI 
Thai 
step 

34 PH * PHAI-STEP 
00 35 K* I , NPHA 
PH - PH ♦$ t E P 
TH « THAI - STEP 
00 35 
PHS PH 

phsth - 0 

TH * THfSTEP 



35 1-J,NTHA 



IF ( ( IGRO.GT.O) . ANO . ( (TH.GT .901. ANO. (TH.LT .270) ) ) GO TO 35 

2GPP,GTf .PH,SG0,SCSP,SCSt,SPPH,SPTH,STPH,STfH,TH,X,Y,2,2LD,ZS,ETA,G 
3 AN , ER R* IGR D) 

ETMAG « CA0S1ETTS) 

EPMAG - CAESIIPPS) 

IF(ETHAG.GT.l.E-32) PHSTH* 57. 295779*ATAN2 ( A IHAGC ETTS)»REAL(ETTS)) 



0209 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 
0290 

0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 
0300 

0309 

0310 

0311 

0312 

0313 

0314 

0315 

0316 

0317 
0310 

0319 

0320 

0321 

0322 

0323 

0324 

0325 

811 $ 

0320 

0329 

811 ? 

0332 

0333 

81 !$ 

0336 



IF( EPMAG. GT. i.E-32 ) PHSPH* 5 7 . 295779* AT AN 2 ( A IM AG( EPPS ) , RE AL ( E PPS ) ) 

IF (NPL.EQ.l) OATYl(K)«EPMAG 

IF (NPL.EQ.l) 0ATY2(K) *ETMAG 

IF (NPL.EQ.2) OATYK I )*EPMAG 

IF (NPL.EQ.2) OAT Y2 ( I ) •ETMAG 

IF (KFLAG(16).NE,i) GO TO 35 

WRITE (6,60) TH, PH, GTT, GPP, ETTS, ETMAG, PHSTH, EPPS, EPMAG, PHSPH 

35 CONTINUE 

WRITE (6,56) 

IF (NPL.LE.O) GO TO 36 
CALL POLPRT ( 1 , DATY1 ) 

CALL POLPRT (2.0ATY2) 

BACK SCATTERING 

36 IF ( ISCAT.LE.O) GO TO 54 
WRITE (6,75) 

WRITE ( 6,80) 

WRITE (6,77) 

WRITE (6,02) 

L - 0 
NPL - -1 
INC - l 

IF (N0AP.EQ.1) GO TO 37 
NPHI - (PHlF-PHl I l/STEPM 
NTHI • (THIF-THllWSTEPfl 
IF (IWR.LE.O) WRITE (6,62) 

GO TO 39 

37 IF (A0AT.GT.5OO. ) GO TO 30 
NPL - I 

NPHI * 360 



NTHI • 1 
“II * 0. 



PHI 
THII 
STEP 
GO TO 39 



ABAT 



30 NPL 
NPHI 
NTHI 
PHII 
THII 
STEP 
39 PH - 



* 360 

- A0AP 

* 0. 

PHh-STEP 



DO 42 K* 1 , NPHI 

‘ StE 



PH • PH*St EP 
TH • THII-STEP 

00 42 1*1, NTHI 



0304 



28 



TH * th*st ep 

IF ( ( IGRD.GT.O) . AND. ( <TH.GT.90I . ANO . < TH . L T . 270) ) ) GO TO 42 

CALL L GFrLO < I A, I B, INC, INM.IWR.il, 12^3* njtWD.N, NO, NH, AH. ACSP,ACST 
l ,C,CGD 1 CG,CJ|CMM,D 1 ECSP,EC$T,EP,eT,EPP, ETT, EP PS, E PTS.E TP$ , ETT S.GG * 
2 GPP, Gif .PH, sdo,s£s£,SCSt ,SPPH,SPTH,STPM, STTH,TH»X,Y,2,ZLO,ZS,ETA,G 
3AM, ERR, iGRD) 

If (NPiIltIo? WRITE *6,63) PH,TH,SPPH,SPTH,STPH,STTH,ACSP,ACST,ECS 



IP.ECST. SCSPt SCST 
40 cphiilI - PH 



CTHET(L) 
CDATKL) 
CDAT2IL) 
C0AT3IL ) 
C0AT4IL l 



TH 
EPPS 
EPTS 
ETPS 
ETTS 



IF INPL.NE.il GO TO 4i 
OATYKKI » CABSIEPPS) 
0ATY2IK) « CABSIEPTSJ 
0ATY3IK) * CABS! ETPS) 
0ATY4IK) * CA8SIETTS) 

GO TO 42 

41 DATYK I) =* CABS! EPPS) 

0ATY2I I ) - CABS! EPTS) 

0ATY3I I I - CABS! ETPS) 

0ATY4I I ) * CABS! ETTS) 

42 CONTINUE 

WRITE (6,82) 

IF (NPL.LE.O) GO TO 43 
CALL POLPRT I 7 , OA TYl ) 

CALL POLPRT (8,DATY2) 

CALL POLPRT (9.0ATY3) 

CALL POLPRT U0.DATY4) 

IF (KFLAGI 17 ) .NE • 1 I GO TO 45 

43 WRITE (6,64) 





DO 


44 


IfcUsI CPHUIl.CTHETdl 


44 


WRITE 




BISTATIC SCATTERING 


45 


IF 


(IBISC.LE.O) GO TO 54 




WRI 


ITE 


(6,75) 




WRI 


ITE 


<6,81) 




WRj 


[TE 


(6,77) 




WRI 


ITE 


<6,021 




WRj 


[TE 


<6, 61} CPH1 (LI ,CTHET( L) 




WRI 


TE 


<6,821 



0305 

0386 

0387 

0388 

0389 

0390 

0391 

0392 

0393 

0394 
U395 

0396 

0397 

0398 

0399 

0400 

0401 

0402 

0403 

0404 

0405 

0406 

0407 

0408 

0409 

0410 
04U 

0412 

0413 

0414 

0415 

0416 

0417 
0410 

0419 

0420 

0421 

0422 

0423 

8 * 5 * 

0426 

0427 

0428 

0429 

0430 

0431 

0432 



,u iu 

>1 )/STEP*l 

Imi'th 



46 



L * 0 
INC - 2 
NPL « -1 

IF (NBIP.EO.l) GO TO 46 
NPHS « ( PHSF-PHS ' ‘ 

NTHS « (THSF-THiL 
IF (IWR.LE.O) WRI 
GO TO 40 
IF ( AB I T • GT *500. I GO TO 47 
NPL - l 
NPHS « 360 
NTHS * I 
PHSl « 0. 

THSl « ABIT 
STEP * 1. 

GO TO 48 



47 NPL 
NPHS 
NTHS 
PHSI 
THSl 
STEP 

48 PH • 



* 1 

- 360 

- A6IP 

- 0 . 

« 1. 

PHSI-STEP 



ANO • ( (TH.GT.90) • ANO. (TH.LT .270) ) ) GO TO 51 



00 51 K-l.NPHS 
PH * PHfSTEP 
TH * THSl-STEP 
IF ( ( IGRD.GT.O) , 

00 51 I * 1 i NTHS 
TH • TH^STEP 
L • L ♦ 1 

CALL GFFLD ( IA,I B.INC . INH.IWR.I 1.12,13, I 12,H0.N,N0.NH, AH.ACSP.ACST 
l.C.CGD.CG.CJ.CMM.O.ECSP.ECST.tP.if ,EPP,ETT.EP PS, EPTS. ETPS, ETtS.GG, 
2GPP.GTT .PH.SGD,SCSP,SCST,SPPM,SPTH,STPH,STTH,TH,X,Y,i,2L0,2S,ETA,6 

3AM, ERR, IGRO) 

IF ( IWR.GT.O) GO TO 49 

IF (NPL.LT.O) WRITE (6,63) PH, TH , SPPH, SP TH, STPN, STTH 



49 CPH((L) - PH 
CTHET(L) « TH 

EPPS 
EPTS 
ETPS 
ETTS 

IF (NPL.NE.l) GO TO 50 
OATYKK) « CABS(EPPS) 
DATY21 Kl - CAB St EPTS) 
DATY3 J K ) - CABSJeTPS) 



COAT I ( L 
CDAT2(L 
C0AT3(L 
C0AT4I L 



0433 

0434 

0435 

0436 

0437 

0438 

0439 

0440 

0441 

0442 

0443 

0444 

0445 

0446 

0447 

0448 

0449 

0450 

0451 

0452 

0453 

0454 

0455 

0456 

0457 
0450 

0459 

0460 

0461 

0462 

0463 

0464 

0465 

0466 

0467 

0468 

0469 

0470 

0471 

0472 

till 

0475 

0476 

0477 

0478 

till 
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CABS(ETTS) 
2) GO TO $1 
CABS* EPPS) 
CA8S1EPTS) 
CABSt ETPS) 
CABS(ETTS) 



0 ATY4( K ) ■ 

50 IF (NPL.NE 
0 ATY 1(1) * 

0 ATY2 C I ) - 
0ATY3U) « 

0ATY4(I) * 

51 CONTINUE 

WRITE ( 6*82) 

IF (NPL.LE.O) GO TO 52 
CALL POLPRT I 3 « OAT Yl ) 

CALL POLPRT 14,DATY2) 

CALL POLPRT 15,DATY3) 

CALL POLPRT (6,DATY4) 

IF ( KFL AG( 18 I • NE • i ) GO TO 54 

52 WRITE (6,66) 



53 

54 



00 53 l * 1 » L 
WRITE (6,65) 



:PHI( I) ,CTHET1I) , COAT 1<I)«CDAT21I)«C0AT31I) , C0aT4( I ) 



55 



IF (IFLAG.EQ.l) GO TO 1 
IF (IFLAG.EQ.2) STOP 

KKFLAG* 0 
K JFLAG-0 
KMFL AG S D 
KNFLAG* 0 

IF (KFLAG( 13) .GT .0) KKFLAG* 1 
IF (KFL AG( 23) .GT .0) K JFLAG* 1 
IF IKFLAG( 14) .GT.O) KMF LAG* 1 
IF (KFLAG( 24) .GT.O) KNFL AG* 1 
00 55 1*1,30 
KFLAG( I ) - -1 

1 



1 



KFLAG<8) 

KFL AGl 2 0 
KFL AG ( 2 6 ) * I 
IF (KKFLAG. GT. 
IF (KJFLAG 
IF (KHFLAG 



IF (IFLAG.EQ.6I 
GO TO 4 



.GT.O) KFL AGC 131*1 
.GT.O) KFL AG( 23 I * l 
.GT.O) KFLAGI 14) * 1 

' ms) 

6,115) 



IF ( KNFLAG.GT • 0 ) KFLAG<24)« 

If (iflag.eq.S) write 



_ ( 6,68 

WRITE (6,11 



FORMAT ( 1H0) 

FORMAT 1 10X, * THE RADIATION EFFICIENCY IS • , F15. 7// 10X, • THE TIME-AV 
1ERAGE POWER INPUT IS • ,F 15 . 7// 10X , • THE ANTENNA IMPEOANCE IS • ,F15. 



0516 



0528 



,E10.4,2(3X,3(E10.4,2X),F6.1,1X 



27,* ♦J , tF15.7//l 

58 FORMAT ( 10X f • IHc NEAR-FIELO ELECTRIC FIELD INTENSITY AT THE OBSERV 
l AT I ON POINT ',E11.5, a ,',Ell.5, »,',E11.5,» (X,Y,Z RESPECTIVELY) IS: 
2 ' / / ) 

59 FORMAT l 2DX , • EX* * , F 1 5 . 7 , • ♦ J * . F 1 5 . 7/ 20X , * EY *• , FI 5 . 7 , • fJ’,Fl5.7/20 
1X.*EZ** ,F15.7,* ♦ J’,F15.7////) 

60 FORMAT ( 3 X , F5 . 1 , 2X , F 5 . I , 3X , E l6 . 4 , 2X , C 
1) ) 

61 FORMAT ( T41 • * FOR BISTATIC SCATTERING THE I NC I OENT • /T41 , • PL ANE WAVE 
I IS PHI* 1 ,F$.l.* THETA** ,F5.1///> 

62 FORMAT (• I NC l OENT • ,T2 7 , « ECHO AR§A SIGH A* ,T66, ’ A8S0RPT ION * ,T90. • EX 
l T INC T ION* ,Til4,* SCATTERING’/' PLANE* ,T2 5,* ( INC I OENT -SCAT TER EO) 1 
24X,3( 5X, 'CROSS SEC T ION* ,6X) /• WAVE • ,52X ,3 ( 10X , • FOR* , 1 1 X) /• PHI 

3 THETA* ,3X,* PH I- PHI*, 3X,* PH I -THETA *,4X, • THETA-PHI * , 2X, • THETA- THETA 
4* ,3( 5X, *PHt* , 7X, * THETA* ,4X) ) 

63 FORMAT ( l X , 2 J F 5 . 1 , l X ) , 1 6 ( E 1 0. 4 , 2 X ) ) 

64 FORMAT ( T54, ' 8ACKSC ATT ER ING* /• I NC I OENT • , T 3 7 , • EL ECTR 1C FIELO POLAR 
1 I Z AT I ON SCATTERING MATRIX*/* PLANE *, T49 ,*( INCI DENT- SC ATTEREO )’/ 3X 
‘ ' WAVE* ,T23, * PHI -PHI* , T49 , * PHI - THE TA ■ , T 7 5 ,* THETA-PHI *, T1 02 ,*THETA- 



C i WAVE ilij* rni-rni f I ' rril - I nt Ifl , I I 3 f 
3THETA* / • PH ( THETA* ,3X,4( 3X. 'REAL* ,8X, * IMAG 
65 FORMAT ( IX , 2 ( F 5 . I , 1 X) , 2X. 4( E 1 1 . 5 . 2X , £ 1 1 . 5 , 3 X ) ) 



8X) ) 



0 3 ru Kn« I liA,£ir3. I , 1 A * fCA;*»(Cit.3,cA,Cil.3,3AII 

66 FORMAT ( T54, * B I S T A T IC * /T3 7, * ELE C TR I C FIELO POLARIZATION SCATTERING 

1 MATRIX*/* OBSERVATION* ,T56,*( I NC I DENT- SC AT TERED j */* P01NT*,14X. 

2 • FHI-PHl * ,T49, * PHI-THETA* , T76 t * THETA-PHI • , T 1 0 1 , • THETA-THE TA * / • P 
3HI THETA* j,4X,4( 3X . * RE AL’ , 8X. » )MAG*,8X) ) 

67 FORMAT (• 68ERSV AT f ON * , T2 t , • £ CHO AREA SIGMA*/* POINT* ,T25. *( INCI 

1 DENT-SCAT TER EO I * /* PHI THETA* , T 1 4 , • PH I -PH I * , T24 , • PHI - THETA • ,T37, 

2 *THET A-PHI * ,T48,*THETA-THETA* i 

68 FORMAT ( 1 H l ,5X , * CfiNT I NUE EXECUTION WITH THE FOLLOWING AOOITIONS AN 

10/0R CHANGES*//) 

69 FORMAT ( 54X* • ELECTR IC FIELO l N TE NSI TY • / 5X , ' OEGREE S * , 1 1 X , • POW ER GA I 
IN* .20 X, * THETA*,42X, 'PHI* /3X,* THETA • ,3X,* PHI *.7X, *THETA*,8X,* PHI • , 1 
2X,i(8X, • REAL’ ,8X. *]MAG* , 8X , *MAGN* ,5X, ’PHASE*) ) 

70 FORMAT ( 10X, ’♦♦♦♦♦ERROR IN OATA CARO NUMBER *,I2,* EXECUTION STOP 
1 PED******* ' ) 

71 FORMAT ( 40X, • ♦ 

1UR ♦ * /40X , • ♦ 

2 RO • * /40X , ’ * 

3 *• ) 

72 FORMAT (40X, ’♦ 

IS ♦ • /4DX , • ♦ 

2R0 /40X , ' ♦ 

3 ** ) 

73 FORMAT 130X,*A 8ACKSC A TTE R I NG CALL MUST BE INCLUOEO FOR A BISTATIC 



A WIRE SEGMENT MAYNOT BE SHARED 8Y MORE THAN FO 

01 POLE MOOES CHECK DESCRIPTION OATA CA 

EXECUTION STOPPEO 

AN ISOLATEO WIRE MUST HAVE AT LEAST TWO SEGMENT 

ANO THREE POINTS CHECK DESCRIPTION OATA CA 

-TOP - "* 



EXECUTION STOPPEO 



1 CALL*//50X, 'REQUEST IGNORED*/////) 

74 FORMAT l • l • , T50 . 3 7 ( ’ •' ) /T50 , • , T86 • *♦ • / 



1 T50,** 

2 T50, 



! DUf J I \ w l f I 5U| * , I OO , ' 

OHIO STATE UNIVERSITY 
ANTENNA ANALYSIS PROGRAM 



*•/ 

♦•/ 



0529 

0530 

0531 

0532 

0533 

0534 

0535 

0536 

0537 

0538 

0539 

0540 

0541 

0542 

0543 

0544 

0545 

0546 

0547 

0548 

0549 

0550 

0551 

0552 

0553 

0554 

0555 

0556 

0557 

0558 

0559 

0560 

0561 

0562 

0563 

0564 

0565 

0566 

0567 

0568 

0569 

0570 

0571 

0572 

0573 

0574 

0575 

0576 
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3 T50* ' * MOOIFIEO FOR USE AT *'/ 

4 T50 » ' * NAVAL POSTGRAOUATE SCHOOL *•/ 

5 T 50» * * 3 SEPTEMBER 1975 *'/ 

6 T50.'*'|T86,'*'/T50#37('*')) 

75 FORMAT ( 1 1 ' , T50 , 29( • * ' ) /T50 , ' * • ,T78 . ' * ' 1 

76 FORMAT (T50,'*' , 1 I X , 'ANTENNA • v Tt8.'*' ) 

77 FORMAT ( T50, '*’ , 8 X , * C A LCULA T I ON S • , T78 * * * • / T 50 , * • ' , T 78 , • * - / T50 * 29 ( • 
1*')) 

78 FORMAT ( T50. 1 *' , 9X , ' NEAR F I E LO ' . T78, • * • ) 

79 FORMAT ( T 50 , 1 • • , 9X , • F A R F I ELO' , T78 ,* *• ) 

80 FORMAT ( T 50, 1 * ' , 7 X , • B ACK SC A T TE R I NG * ,T 78 , • * • ) 

81 FORMAT (T50,'*' , 4 X , ' 8 l S T A T I C SC ATTER I NG * ,T78t ' ** ) 

82 FORMAT (////) 

83 FORMAT ( T 30, ’ FRE QUENC Y ( MH2 I • , T 8 1 . E 1 I . 5 I 

84 FORMAT (T30,'WIRt RAOIUS ( ME T E R S ) 1 . T 8 1 . E 1 1 . 5) 

85 FORMAT (T30,'WlRE CONDUCTIVITY ( M EGAMHO S /ME T E R ) • , T 8 1 , E 1 1 . 5 1 

86 FORMAT (T30,'WIRE INSULATEO I NO/ Y E S I * , T 8 5 , ' YE S * ) 

87 FORMAT (T30,'WIRE INSULATEO < NO/ YE S) 1 ,T 8 5 , ' NO* ) 

88 FORMAT CT30, ’INSULATION RAOIUS l ME TERS > * ,T 8 1 , E 1 1 . 51 

89 FORMAT (T30,» INSULATION CONDUCTIVITY (MHOS/ HE TER) ' «T81«E11.5) 

90 FORMAT (T30, ’ INSULATION DIELECTRIC CONSTANT ( REL AT I VE ) • , T8 1, E 11 .5 * 

91 FORMAT <T30,' INSULATION LOSS T A NGE NT ' . T 8 l.E 1 1 . 5 ) 

92 FORMAT (T30, 'EXTERIOR MEO I UM • , T 8 1 , ' FREE SPACE') 

93 FORMAT (T30, ’EXTER IOR MEOIUM CONOUCTIVITY ( MHOS/METE R) '. T8 1 • E 1 1 . 5 ) 

94 FORMAT ( T30, ’EXTERIOR MEOIUM OIELECTRIC CONSTANT ( RELAT 1 VE ) • , T81 ♦ 

1 Ell. 5) 

95 FORMAT (T30, 'EXTER IOR MEOIUM LOSS TANGE NT ' ,T8 1 , E 1 1 . 5 ) 

96 FORMAT ( T50 , ' W I RE STRUCTURE ’ / / T 2 0 , ' SEG' , 4X , 2( • NOO E ’ , 19 X , ' LOC A T I ON • 
1 t l8X) /T21. 'NO. • ,3X.2( ’ NO. ' ,9X. ' X’ .13X, * Y' , 1 3X , ' 2 ' , 7 X) / ( T2 1 , 1 2, 5X , 
22 ( I2,5X T Eil.5,4X, E11.5,4X,Ell.$, IX))) 

97 FORMAT (T50, 'ANTENNA F EE 0 S' / T40 , ' NOOE • , l 6X , ' VOL T S' /T4l » ' NO. • , 12X , 
l ’REAL' , 7X,* IMAGINARY' /(T41 , I 2 , 6X , 2( 4X , E 1 1 . 5) I ) 

98 FORMAT (T50,'*«, 6X,'0UTPUT RE QUE STEO' ♦ T 78 , ' * ' ) 

99 FORMAT ( T 30, ’ STRUC TURE CURRENTS') 

100 FORMAT (T30»'FAR FIELOS FOR PHI VARYING FROM' , IX , F5. 1, ' TO ' # F5.1t 
1 'ANO THETA VARYING FROM ',F5.1i' TO ',F5.1/ 

2T50 , ' IN STEPS OF '.F5.1.' OEGREE S. • ) 

101 FORMAT (T30,’ B AC K SCAT TER I NG FOR PHI VARYING FROM ',F5.1*' TO ',F5. 
lit* ANO THETA VARYING FROM ',F5.1,' TO ',F5.1/ 

2T50* ' IN STEPS OF '»F5.1.' DEGREES. ') 

102 FORMAT ( T30,' 81 STATIC SCATTERING FOR PM I VARYING FROM '*F5.1#’ TO 
1 ' » F5 • 1 , • ANO THETA VARYING FROM '*F5.1t' TO ',F5.1/ 

2T50 , ' IN STEPS OF '.F5.1,* DEGREES.') 

103 FORMAT (T30,'NEAR MELOS FOR FOLLOWING POINTS ( X • Y • Z ) • / 501 T40 , 3( E 1 

LD THETA-' .F5. 1) 



1 1 • 5, 5X) II 

104 FORMAT (T30,'PLOT FOR FAR FIELD THET A-' *F5. 

105 FORMAT (T30,'PLOT FOR FAR FIELD PHI-'.fS.D 

106 FORMAT (T30,'PLOT FOR BISTATIC SCATTEftlNG-F 



G-FOR THETA-',F5.U 



0577 

0578 
05 79 

0580 

0581 

0582 

0583 

0584 

0585 

0586 

0587 

0588 

0589 

0590 

0591 
059? 

0593 

0594 

0595 

0596 

0597 

0598 

0599 

0600 
0601 
0602 

0603 

0604 

0605 

0606 

0607 

0608 

0609 

0610 
0611 
0612 

0613 

0614 

0615 

0616 

0617 

0618 

mi 

0621 

0622 

0623 

0624 



107 FORMAT ( T 3 0 , ' PL OT FOR 8ISTATIC SCATTERING FOR PHI-',F5.I) 

108 FORMAT (T30, 1 PLOT FOR 8 ACK SCAT T E R ING THE T A* • , F5 . [ ) 

109 FORMAT (T30t*PL0T FOR SACKSCATTERING PHI-'.F5.1I 

110 FORMAT ( T 30, • NO OUTPUT OR PLOTS REOUESTEO'I 

111 FORMAT (//) 

112 FORMAT <T50,'*' , T 78 , ' * • /T50 . 29 ( • *» )} 

113 FORMAT ( T 50 , • • ' » 8X, 'INPUT DATA ',T78,'*') 

114 FORMAT (T50t29(*** I / T 50 ? ' ♦ ' * T7b • ' * ' ) 

115 FORMAT (lOXj’SINCE THIS OATA SLOCK OOES NOT HAVE A TERMINATION CAR 
10 A change Carj IS ASSUMEO') 

116 FORMA! (//10Xi40< ' *’ ) / 10X , ' THE DESCRIPTION ANO THE GEOMETRY OF THE 
1 STRUCTURE'/IOX* ’MUST BE STATEO IN THE FIRST OATA 8L OCK . * / 10X , • *** 
2* EXECUTION STOPPEO ***') 

117 FORMAT (// lOX.'NO PART OF THE WIRE STRUCTURE CAN LIE 8EL0W THE GRO 
l UNO PLANE.' /10X, ’♦♦♦♦EXECUTION STOPPEO****') 

118 FORMAT ( T50, • STRUCTURE LOAO S’ / T 40 , ' NOOE • , 1 6X , ' OH MS • / T4 1 . • NO . • , 1 2X 
1 .'REAL' ,7X, ’ I MAGI NARY' /( T4 1 , I 2 , 6X , 2( 4X , E 1 1 . 5 ) ) i 

119 FORMAT (T50, ' STRUCTURE LO AO S * / T i9 , * SE GHENT ' , 1 4X , ' OHMS ’ / T4 1 , • NO' , 1 2 
lX.'REAL' . 7X, ’ IMAGINARY' /(T41, I 2.6X,2(4X, El 1 •511) 

120 FORMAT ( T 50 ♦ ' ANT ENNA FEEOS' / T39 ,' SEGMENT 1 4X . ’VOLTS' /T41 , 'NO. '. 12 
IX .'REAL' ,7X. ' IMAGINARY' /( T41 , 1 2«6X,2(4X,Eil .51)1 

121 FORMAT (//no, 'THE NUM8ER OF INTERVALS FOR CALCULATING THE ELEMENT 
1 S' /T30* • I N THE IMPEOANCE MATRIX WITH S I M PSONS-RUL E INTEGRATION IS* 
2./T30.13,’. IF CLOSE 0 FORM INTEGRATION IS REQUIRED SET INT-O'///) 

122 FORMAT (T30,'GROUN0 PLANE ( NO/YE S) ' ,T85 ♦ ’ NO ' ) 



123 FORMAT (T30, 'GROUND PLANE ( NO/ Y E S ) • ,T8 5 , ’ Y E S • ) 

124 FORMAT (T30,'GR0UN0 OIELECTRIC CONSTANT (RELATIVE) ',T81,E11.5/ 

1 T30, ' GROUNO CONOUCTIVITY I MHOS/ME TE R ) ' , T8 I . E 1 1 . 5 ) 

125 FORMAT (T30, 'GROUNO PLANE • .T83 .' PERFECT • I 

126 FORMAT j T30.' ANT ENNA HE IGHT (METERS) ',T8l«Ell*5l 

127 FORMAT ( // 10X .40 ( • *• 1 / 10X , • THE WIRE RADIUS MUST BE STATED' /10X,40( 



l'*'l) 

END 



0625 

0626 

0627 

0628 

0629 

0630 

0631 

0632 

0633 

0634 

0635 

0636 

0637 

0638 

0639 

0640 

0641 

0642 

0643 

0644 

0645 

0646 

0647 

0648 

0649 

0650 

0651 

0652 

0653 

0654 

0655 

0656 

0657 
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BLNK 



PURPOSE: to compress data to the left by removal of 

the blank spaces on the input data cards. 

METHOD: A (I) character is compared to the blank; and 

if it is true, the A (1+1) character is shifted to the A (I) 
position. 



CALLED BY: READ 



CALLS TO: NONE 



SUBROUTINE BLNK (A) 
DIMENSION A ( 80 I 
DATA BLANK/' '/ 

K - 0 

oo t i*i,eo 

J * t-K 
AIJ) - Alt) 

l IF iAI l ) • EQ. BLANK ) K-KM 

IF IK.EQ.O) RETURN 
A<8l-K) • BLANK 
RETURN 
ENO 



1 

2 
3 
A 

5 

6 
7 
0 
9 

10 
11 
12 
13 
1 A 
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CBES 



PURPOSE: to calculate the quantity B01 where 

B0 1 = J (z) / J (z) . 

0 1 

METHOD: If the absolute value of the argument for the 

Bessel function is less than 12, B01 is calculated via the 
power series expansion for the Bessel function in the DO 
LOOP ending at statement 3. If greater than 12, the 
asymptotic expression is utilized at statement 4. If the 
magnitude of the complex part of the argument for the Bessel 
function is greater than 20, B01 is set to (0.,-1) . If the 
complex part of the argument is negative, the sign of E01 is 
changed prior to returning to the calling program. 

CALLED BY: SGANT 

CALLS TO: NONE 



subroutine cees a, son 

COMPLEX ARG.CC.CS.fX 

COMPLEX 601.2. f ERMJ.TERHN.M22*, JNI2) 
OATA PI /3.1* 159/ 

IF <CA8SI2).GE. 12.0) GO TO * 

FACTOR * 0.0 
TfckHN = (0.,0.) 

HZ 24 * -0.25*Z*Z 



T ERMJ « ( 1 .0,0.0) 




9 
1 0 


00 3 NP = l » 2 




11 


N = NP — 1 




12 


J N ( N P J * TERHJ 




13 


H * 0 




l* 


H = H*l 




15 


TERHJ * TERHJ*HZ24/FLOAT(M*(N*M) ) 


16 


JN(NP) = JN { N PI ♦■TERHJ 


17 


IP (NP.NE.l) GO TO 


2 


18 


FACTOR = FACTOR* i . 


0/FCOAT(H) 


19 


TERHN = TERMN*TERMJ*F ACTOR 


20 


ERROR = CABS(TERHJ) 


21 


IF < ERROR. GT . 1.0E- 


10) GO TO l 


22 


TERHJ = 0 . 5* Z 




23 

2^ 


801 * JN(i)/JN<2) 




25 


RETURN 




26 


Y *= AIHAG(Z) 




27 


IF (A8S( Y) .GT.20. 1 


GO TO 5 


28 


ARG *= ( .0, 1. )*Z 




29 


EX *= CEXP(ARG) 




30 


CC * EX M ./EX 




31 


CS - ( . 0*- 1 • I *( EX- 


i./EX) 


32 


80 1 * (C§*CC)/(CS- 


CC) 


33 


RETURN 




3* 


B01 » (.0,-1.) 




35 


IF (Y.LT.O.) 801 - 


(.0,1.) 


36 


RETURN 




3 7 


END 




38 
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DSHELL 



PURPOSE: to calculate the mutual impedance term 

contributed by the dielectric insulation on the surface of a 
thin wire. 



METHOD: The contribution to the impedance matrix is 

calculated utilizing the equation below 



z 

mn 




- &) In (b/a) 
j w&& 

2 



/ 



m, n 



F* (1) F' (1) dl , 
m n 



where z is defined in subroutine SGANT, & is the 
mn 2 

dielectric constant of the insulation, b is the outer radius 

of the insulation, a is the inner radius, & is dielectric 

constant of the external medium, and F is the sinusoidal 

expansion function. 



CALLED BY: SGANT 



CALLS TO: NONE 



SUBROUTINE OSHEL L I AM, BM, OK ,CGOS , SGOS.E P2. EP , ET A ,GAM,P 1 1 ,P1 2 I 
complex cgds.sgds, ep2,ep,eta,gam,pu .pi 2 .fio.csr 

DATA PI /3.IA159/ 

GD - GAM*DK 

C ST • ( EP2-EP)*ETA*AL0G(BM/AM)/(4.*Pl*EP2*SGDS*$GDSI 
Pli • -CSr*<GD*SGDS*CGDS» 

P12 - CST*(GD*CGDS*SGDS) 

RETURN 

END 



3 

4 

5 

6 
7 

e 

9 
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EQUAL 



PURPOSE: to determine position (location) of the equal 

symbol on input data card. 

METHOD: The character search begins in the column 

passed to the subroutine. On returning to the calling 
program r the argument passed is the column following the 
equal symbol. 

CALLED BY: READ 

CALLS TO: NONE 



SUBROUTINE EQUAL <N) 
INTEGER A « EQULS 
COMMON /A/ A(80l 
OATA EQULS/* *• / 

K - N 

00 1 I*K$ SO 

n * m 

IF (A(I).EQ. EQULS) GO TO 2 

1 CONTINUE 

N - l 

2 RETURN 
END 



0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 
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EXP J 



PURPOSE: to calculate the exponential integral with 

complex limits. 

METHOD: The exponential integral is defined as: 



W 12 




dv = E^ (VI) - E i (V2) + j2nti 



9 



where the integration path is the straight line from V 1 to 
V2 on the complex v plane and 

E t (z) = 

z 




The integration path is a horizontal line in the w plane 
or an inclined straight line from VI to V2 the v plane. The 
integer n is zero unless this path intersects the negative 
real v axis at a point between VI and V2. When there is 
such an intersection, 

a) n = 1 if In (VI) > Im (V2) 

b) n = - 1 if Im (V 1 ) < Im(V2). 

The term j2nn is caluclated below statement 12. 

CALLED BY: GGMM 



CALLS TO: NONE 
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SUBROUTINE EXPJ <V1,V2,W12) 

COMPLEX EC,fcl5,S.T f UC.VC,Vl,V2,Wl2,Z 
01 ME NS ION VI 21), W ( 2 1 ) , 0(16), E(l6) 

DATA V/ 0.2 220466 76 00,0. 11009321601 ,0.2992 7363601 , 0 . 5 77 5 14366 0 1 ,0.9 

10374674601 .0. 15902074602.0.933070126-01 , 0 . 4 9269 1 74 EOO, 0 . 1 2 1 5 5 95 46 0 

21.0. 22699495601 . 0 . 366 76 Z2 76 0 1 , 0 . 5*2 5 3 36 6 EO 1 , 0 . 7 5 6 59 1 62 EO 1 , 0 . 1 0 1 20 2 
32 0EO2,O. 13 1302 32 E 02,0. 1 66 5440 0E 02 , 0 . 20 7 7647 9E 02 , 0 . 2 5 623 09460 2 , 0 . 3 1 
*4075 19E02, 0. 305 3O603EO2 ,0. 48026006602/ 

OAT A W/ 0.45 896*60 EOO, 0.41 700083 6 00,0. 1 1 3 3 7 3 30EOO . 0 . I 03 99 1 976 -01 ,0 . 
126101 72 06-03, 0. 09054 791 E-06,0. 2 1 0 2 340 7 EO 0 , 0 . 3422 101 7600,0.263027 50 
2 E 00,0. 1 264 2 5 826 00 , 0.402 06 06 5E- 01 , 0. 0 56 3 0 77 06-02 , 0 , 1 2 12 4 3 6 IE -02 , 0 . 1 
3 1 16 74406-03,0. 64 5 9926 7E -05 ,0.2 2 263 169E- 06,0.4 2 2 74 3O4E-O0 ,0.3921097 
436-10,0. 14565 l 52E- 12 ,0. 140302706-15,0. 16 O05949E- 19/ 

DATA 0/0.22495042602,0.7441 1560602,-0.41 431576603 ,-0.7 8 75 4339 E02.0 
1. 11254/446 02,0. 1 6 02 176 IE 03, -0.23 06 21 9 56 0 3,-0.5009 46 0 76 03,-0. 6 040 78 

254602.0. 12254 770E 02,-0. 101619766 02 ,-0.4 72 19 59 1E01 ,0.79 72 9681 601,-0 
3 . 210695 74 EO 2 , 0 . 2 2 0464 90E0 t ,0.09 720244EO1/ 

OATA E/0. 211031076 02,-0.3 7959 707EO3,-O.9 74 0922O6O2,O.12 9OO6 7 2 6O3.O 
1 .179492 266 02,-0. I 29 1 0 9 3 1 E 03 .-0 . 5 5 705 5 746 03 . 0 . 13 52 4 &01E 02 .0 . I 4 696 7 2 
2) 60 3,0. 179495 2 06 02 , -0 . 3 29 0 1 6 1 46 00 , 0 . 3 10 2 00 36E 02 , 0. 8 1 65 76$ 7E 0 l , 0 . 2 2 
3 2 3696 IE 02,0. 3 91 24 09 26 02 , 0 . 8 1 6 36 79 9E0 1/ 
l * VI 

00 12 J I M= 1,2 
X - REAL! Z) 

Y = A 1 M AG ( Z ) 

E 15 - < .O..OJ 
A0 - CABSi Z) 

IE (A0.EQ.O.) GO TO 11 

IE (X.G6.0. .AN0.A8.GT • 10. ) GO TO 10 

YA - A 6 S ( Y ) 

IE (X.L6.0. .ANO.YA.GT. 10. ) GO TO 10 

IE (YA-X.G6.17.5.0R.YA.GE.6.5.0R.X+YA.G6.5.5.0A* X.G6.3.) GO TO 2 
IF (X.LE.-9.) GO TO 6 
IE (YA-X.GE.2.5) GO TO 7 
IE ( X* Y A . GE • 1 • 5 f GO TO 3 
N * 6 • ♦ 3 • ♦ A8 
E 1 5 - l./(N-l.)-Z/N**2 
IN* N-l 

E 1 5 - 1 ./(N-l • )-Z*El5/N 
IE (N.GE.3) GO TO 1 
5-C * 



E 15 - Z *E 1 
GO TO l l 
J1 - 1 
J2 - 6 
GO TO 4 
J1 - 7 
J 2 - 21 



-CMPlX(.5772l6*AL0G(A8)*ATAN2(YvX)) 



1 

2 

3 

4 

5 

6 
7 
0 
9 

10 

11 

12 

13 

14 

15 

16 
\ 7 
10 

19 

20 
21 
22 

23 

24 

25 

26 
27 
20 

29 

30 

31 

32 

33 

34 

35 

36 

37 
30 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 



S * 
YS 



( . 0 ,. 0 ) 

* Y*Y 



00 5 l ~ J l , J2 

XI = VC l > ♦ X 

CE * W(l)/(Xl*Xl*Y$) 

5 S - S *C M PL X ( X l * C 6 ,-YA*CE) 

GO TO 9 

6 T 3 3 X * X- Y * Y 
T 4 * 2 • * X * Y A 

T 5 = X* T 3- YA ♦ T4 
T 6 = X* T4*YA*T3 

UC * CKPLXtDC 11 ) *0(12 J *X*0( 13)*T3*T5-E( 12I*YA-E( 13J*T4«E(ll l*E( 12J 
1 • X ♦ 6 ( 13)*T3*T6*0( 1 2 ) * Y A *0 ( l 3 1 * T 4 ) 

VC = CMPLX(0(l4l*0(l5)*X*D(16)*T3*T5-E(15)*YA-E(16I*T4 t E(l4l*E(15) 

1 * X ♦ E ( 16 )*T3*T6*0( 15 ) * Y A *0( 1 6 ) * T 4 ) 

GO TO 8 

7 T 3 = X * X- Y * Y 
T 4 = 2 . * X* YA 

T 5 * X* T3-YA * T4 
T 6 * X*T4*YA*T3 
T 7 * X*T5-YA*T6 
T 0 » X * T6f YA*T5 
T9 ' X*T7-YA*T0 
T 10 * X*Tft*YA*T 7 

UC = CNPLX(D( 1) *0(2)*X*0(3)*T3*0(4)*T5*0(5)*T7*T9-(E(2)*YA*6( 3 ) * T 4 
1*6(4) *T6*6(5)*T0) ,6 (1) *E(2)*X*E(3)*T3*6( 4 ) * T 5*E ( 5 ) * T 7* T 1 0* ( D ( 2) *Y A 
2*D(3)*T4*0(4)M6*D(5 ) *T8 ) ) 

VC = CMPLX(0(6) *0(7 )*X*0(8)*T3*D(9)*T5*D( 10)*T7*T9-(E( 7 J*YA*E(8)*T 
l4*E(9)*T6*E(lOi*T0),E(6)*E(7)*X*E(8)4T3*E(9)*T5*E(lO)4T7*TlO*(O(7) 
2*YA*O(0)*T4*O(9) *T6*0(10)*T8)) 

8 EC * UC/VC 

S * EC/CMPLX(X, YA ) 

9 EX - EXP(-X) 

T * 6X*CMPLX(C0S(YA) f -SlN( YAH 
E 1 5 * S*f 

IF (Y.IT.O.) 615 * CONJG( 615) 

GO TO 1 1 

10 E 1 5 - .409319/(7*. 193044)*. 421031/(Z*1. 02666)*. 147126/(Z*2. 56708)* 
1 • 206335E- 1 / ( Z* 4. 9 003 5 ) * . 1 07401E-2/( Z* 0. 1 8 2 1 5 ) * . 1 5 865 46 - 4 / ( Z* 12.734 
22)*. 317 03 IE- 7/( 7*19.39 57) 

E 1 5 - E15*CEXP(-Z) 

11 IE (JlM.EQ.il W 1 2 - E 1 5 

12 Z - V2 

Z - V2/V1 

TH - AT AN2(A1NAG( Z ) • RE AL( Z) ) -AT AN2 ( A 1 M A G ( V 2 ) , RE AL ( V2 )) ♦ A T AN2 ( Al MAG 
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72 

73 

74 
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76 

77 

78 

79 
00 
01 
02 

03 
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07 
00 
09 

90 
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96 



1 (VI) ,RE AMVl )) 

A6 * A0 S ( TH ) 

IE (A0.LT.1.) TH • .0 

^ 2 T uSN Ml2 - Et5tCMPLX< '°* IH? 

END 
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98 

99 
1 00 

IS) 

103 

104 
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GANT 1 



PUEPOSE: to consider the wire structure as a 

transmitting antenna and calculate the input impedance and 
current distribution. 

METHOD: If a wire antenna is driven by a voltage 

generator v located at one of the current sampling points 
i 

1 and if displacement currents are neglected. Ampere's law 
i 

yields 



V = v F (1 ) 
m i m i 

where F is the sinusoidal expansion function. Thus, the 

excitation voltages V will vanish everywhere except where 

m 

v is not zero, 
i 



The DO LOOP ending with statement 2 uses the delta-gap 
model defined above to determine the excitation voltage 
CJ (I) for all the dipole modes. These are stored 
temporarily in CG(I). Then subroutine SQROT is called to 
obtain a solution of the simultaneous linear equations. 

SQROT stores the solution (the loop currents) in CJ (I) . 

In the DO LOOP ending at statement 6 , the complex power 
input and input impedance (s) are calculated. The 
time-average power input (PIN) is the real part of the 
complex power input. 

Subroutine RITE is called to make the transformation 
from the loop currents to the branch currents. If IWR is a 
positive integer, RITE will write out the list of branch 
currents. 
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Finally, GANT 1 calculates the 
calling subrouinte GDISS to obtain 
dissipated in the lumped loads and 
conducting wire. 



radiation efficiency by 
the time-average power 
the imperfectly 



CALLED BY: MAIN 



CALLS TO: 



GDISS 



RITE 



SQROT 



SUBROUTINE GANT l I I A, 1 8, INN, IMR, 1 1 , I 2, 1 1 , 1 1 2, JA. JB. NO. N, NO.NN. AH.C 
l.CJ.CG.CMH.D.EFP ,GAN,GG.CGD.SGD.VG.YlI,2ll.2 LD.*S. IgRD I 
COMPLEX YY , CGEN 

COMPLEX Cl il.CJIl l,CGD(II,SGD(ll .VGIU.UDI II , Y II . Ult IS ,GAN, CGI 1 I 

o i mens i on ot I T. iAln. ibiiT. jaIu, j&iiT 
DIMENSION Hill. 12(11. lilil. MOl INN.41 « NDI1I 
COMMON I ML 



DO 3 1*1, N 
CJIII • 1.0, .01 
K • JAI II 



DO 2 XtC* 1.2 
KA * I A ( K J 
KB « I B ( K ) 

IF (KB.EQ.I2U)) GO TO 1 
IF (KB.EQ.IlUl) FI — l. 
CJU) * C J ( I) +F I *VG( J J ) 
GO TO 2 

IF ( KA • EQ. 13 ( II ) FI — l. 
JJ • K+NM 

CJU) - C J (I ) +F I *VG( J J ) 

K - J8( I ) 



3 CONTINUE 



00 4 I* I • N 
CGU ) - CJU) 



CALL SQROT (C.CJ ,0 ,1 12 ,N) 

112 * 2 

Y1I « ( .0*.0> 

NNN * N 

IF (IGAO.GT.O) NNN « (NUWU/2 



DO 6 I * 1 * NNN 
NN ■ I Ai JB ( I ) ) 

CGEN*CGEN/2. 



0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

m 



YY « CJ ( I ) +C0N JG (CGEN) 

IF (CABS(YY) .Lt • I •E**20) GO TO 5 
Zll - ( l. W)*(CAB$(CGEN)**2) 



WRITE (6.81 NN » Z 1 1 

5 Yll - Y 1 i ♦YY 

6 CONTINUE 
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IF ( IWR.GT.O) WRITE (6.7) 

CALL RITE (U.IBf INN, IWR, 1 1 , 12, 1 3* MD.ND. NM, CJ.CGt IGRD) 
GG - REAL(Yll) 

IF UGRO.GT.O) GG-2.*REAL(Yil) 

PIN « GG 

CALL GOISS ( AN.CGt CNN.D.DISS.GAH.NH.SGOf 2L0.ZS) 

PR AO ■ PIN-DISS 
EFf « 100.*PRAD/PIN 
RETURN 



7 FORMAT ( 50X, 9 ANTENNA BRANCH CURRENTS*) 

6 FORMAT ( lOXt • THE INPUT IMPEDANCE AT NODE ’»I3t 9 IS*,F15.7,* ♦ J*, 
IF 15*7/ / ) 

ENO 
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GDISS 



PURPOSE: to calculate the time-average power 

dissipated in the imperfectly conducting wire and in the 
lumped loads. 

METHOD: The time-average power dissipated by the wire 

is calculated in the DO LOOP ending at statement 1 utilizing 
the equation below: 

1 

B f * 

P = s_ I I dl 

d 2na J 

where R is the surface resistance of the wire and a is the 
s 

radius of the wire. 

The power dissipated by the lumped loads is calculated 
by the DO LOOP ending at statement 3. If the wire is 
perfectly conducting, CMM <0, the first calculation is 
by-passed. 

CALLED BY: GANT1 

CALLS TO: NONE 
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SUBRCUTIN6 GOISS < AM, C G , CMM t O t 0 I SS. G AM, N H. SGD , ZL 0 , Z S ) 

COMPLEX CGU) .SGO< U tZLDU) »CJA,C J6,GAM, zi 
DIMENSION D ( I ) 

DATA P 1/3.14159/ 

0 1 SS * .0 

If icmm.le.o.i GO TO 2 

ALPH * REAUGAMJ 
BETA * AIMAG<GAM) 

RH * REAL! I S I / 4 4,*Pl*AM) 

DO I K*1,NM 
DK - DIK) 

DEN * C AB S ( S GO I K ) )*+2 
EAD * EXP< At PH* DK I 
C AD * < EADM./EADI/2. 

CBD « COS(BETA*DK) 

SAD * DK 

IF (ALPH.NE.O.) SAD*( EAD-l ./EADI/(2.*ALPH» 

SBD = DK 

IF IBETA.NE.O. ) SB 0* S l N < B E T A*OK » / BET A 
FA * RHM SAD*CAD-SBD*CBD»/DEN 
FB * 2.*RH*<CAD*S8D-SAD*CeD>/DEN 
CJA * CGIKI 
L * K* NM 
CJB * CG(l) 

1 DISS * DI SS*FA*( CABS! CJA) **2* , CABS(CJB1**2 )*FB*(REAL(CJA l+REALICJB ) 
1 * A I MAGI CJA)*AIMAG(CJ3 ) ) 



2 00 3 J * I » NM 
K * J*NM 

3 DISS ■ D1SS+REAL(ZLD( J >)» < C ABS ( C GU )> ** 2 > + REALI ZLD< K H * < C A8S( CG ( K ) 
I > ♦ ♦ 2 1 

RETURN 

ENO 



l 

2 

3 

4 

5 

6 

7 

8 

9 

10 

I l 

12 

13 

l* 

15 

!* 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 
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GFf 



PURPOSE: to calculate the far-zone field of a 

sinusoidal electric monopole. 

METHOD: If an electric line source has length d and 

endpoints at (x ,y ,z ) and (x ,y ,z ), then the coordinates 

111 222 

of any point on the source are 

x = x +1 cosx 

1 

y = y ^ +1 cosy 

z = y + 1 cosz 
1 

where cosx, cosy, cosz are the direction cosines of the 1 
axis, and 1 is the distance along the source measured from 

the endpoint (x ,y ,z ). Let the current distribution on 

11 1 

the monopole be 



I sinh y (d - 1) + I sinh yi 

1 ( 1 ) = 1 1 

sinh yd 

where I and I are the endpoint currents. The far-zone 

1 2 

field of this source is 

E = (cosx cos9 cos<J> - cosy cose sin<j> - cosz sine) E 

1 



E = (-cosx sin<j> + cosy cos(J>) E 

e 1 



where 



E = 
1 



-yr 

rje 



4nr (1-g 2 ) sinh yd 



+ 




ygd y f<*> 

[ (e - g sinh yd - cosh yd) I e 

yf<2 ) 

+ g sinh yd - cosh yd) 1^ e ] 
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f < 1 > = x sin© cos* + y sin© sin<|> + z cos© 

1 1 1 

f< 2 > = x sin© cos$ + y sin© sin<t> + z cos© 

2 2 2 

g = cosx sin© cos<t> + cosy sin© sin<J> + cosz cos© 

and (r,©,#) are the spherical coordinates of the observation 
point . 

In this subroutine the range dependence has been 
suppressed. The far field vanishes in the endfire direction 
where GK = 0. If a ground plane is present (IGRD > 0) the 

equation above is decomposed into the x # y, and z components 

and the reflection coefficients are applied before E and E 

& <j> 

field components are returned to the calling program. 

CALLED BY: GFFLD 

CALLS TO: HONE 
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SUBROUTINE OFF ( X A , YA . Z A. XB , YB , Z B , Or CGO » SGO * CTH* STH ,CPH, SPH,GAN,ET 
l A* ET l ,ET2 *EPl » E P 2 » I G RO iEArI 
CCHPLEX ERR.RV.RH.RR.EXtEY, EZ.EE 
CGMPLEX ETI.ET2 .EkI, §P2 ,GAM, Ei A 
COMPLEX GO.CGO.SGO.EGO 
COMPLEX EGFA,E&FB,£gGD,ESA,ESB 



COMPLEX CST 
FP • 12.56637 
XAB * XB-XA 
YAB - YB-YA 
ZAB « ZB-ZA 
CA - XAB/O 
CB ■ YAB/D 
CG • ZAB/D 

G * ICA*CPH*CB*SPH)*STH*CG*CTH 
GK ■ l.-G*G 
ETl ' * ‘ 

ET2 

EPl 

E P2 ■ ( .0, .01 
IF IGK.LT..OOI) GO TO 3 
FA • IX A*CPH*YA*SPH|*STH*ZA*CTH 
FB * (XB*CPH*YS*SPH ) *STH*Z B*CTH 
EGFA * CEXPI GAM*F A ) 

§GFB ■ CEXPI GAM* F B ) 

:GGD » CEX PI GAM* G*D I 



L . 

1 .0. .01 

( .0,.0I 

1 .0. .0) 



:st 

ESA 



eta/igk*sgd*fpi 



CST* EGF A* ( E GGD-G*SGD-CGD I 

( l./EGGD*G*SGO-CGDl 

GO TO 2 



- I-1..0I 

IIGRO.EQ.L) GO TO l 
• CSORTIERR- STH*STHl 



E SB * C ST*EGFB* ( 

IF I IGRD.LE.OI 
RV - I-l.tOI 
RH - I-l.tOI 
IF 

RV - -fERR*CTH-RR J /fERA*CTH*RRI 
RH ■ ICTH-RR l/ICTH+RR) 

EX ■ CA*ESA 
EY - CB*ESA 
EZ - CG*ESA 

EE » (EX*SPH-EY*CPH|* IRH-RV) 

EX • EX*RV*E E* SPH 
EY • EY*RV-EE*CPH 
EZ « -EZ*RV 

ESA — EX*CA-EY*CB*EZ*CG 
EX ■ C A*ESB 

II : nnh 

EE • I EX*SPH-E Y*CPH) * IRH-RV I 



0001 

0002 

0003 

000* 

0005 

0006 
0007 
OOOB 
0009 

881 ? 

0012 

0013 

001* 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 
0023 
002* 

0025 

0026 
0027 
002 B 

0029 

0030 

0031 

0032 

8812 

0035 

0036 

0037 

8811 

00*0 

00*1 

00*2 

00*3 

00 ** 

00*5 

m 

00*8 



EX ■ EX*RV*EE*SPH 
EY • EY*RV-EE*CPH 
EZ - -EZ*RV 

E SB — EX*CA-EY*CB*EZ*CG 

2 T - ICA*CPH*CB*SPHI*CTH-CG*STH 
P - -CA*SPH*CB*CPH 

ETl ■ T*ESA 
ET2 - T*ESB 
EPl - P*ESA 
EP2 • P*ESB 

3 CONTINUE 
RETURN 
END 



8858 

0051 

0052 

0053 
005* 

0055 

0056 

0057 

0058 

0059 

8 060 
061 
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GFF1D 



PII8P0SE: to calculate the far-field for the thin wire 

structure. 

METHOD: The far-field for the structure is calculated 

from the loop currents. The loop currents are either the 
currents produced by the transmiiting antenna calculations 
of subroutine GANT1 or the currents produced by an incident 
plane wave. 

If the incident field is generated by a distance source 
with spherical coordinates (r^,0^,4>^) , the excitation 

voltages induced by a incident plane wave are 




E = E exp (Y r-r ) 
i 0 0 

where E is a vector constant, r is a vector from the 
0 0 

coordinate origin to the distance source, and r is the 
radial vector from the origin to the observation point. 



The field E is generated by test dipole m when 
M 

radiating in the homogeneous medium. Using the vector 

potential, the field at the distance point (r ,0 , ) is 

0 0 0 



E 

m 



-/ r 

jwu e 0 



4nr 

o 

m 




exp (Y rf ) 
o 



dl 
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where the radial component is to be suppressed. From the 
above equations, 

4nr /r 

V = 2- e ° e e 

m jwu ° m 

If an antenna gain calculation is desired, INC is set to 
zero. PH and TH denote the spherical coordinate direction 
of the distance observation point. The phi-polarized (EPPS) 
and the theta-polarized (ETTS) components of the electric 
field intensity are returned to the calling program. 

If INC = 1, a backscattering calculation is desired. In 
this case PH and TH denotes the incident angles for the 
incident plane wave. These are also the spherical 
coordinates of the distance source. The outputs returned to 
the calling program include absorption, extinction, and 
scattering cross section for each polarization; scattered 
electric field; and echo areas. 

If INC = 2, a bistatic calculation is desired. In this 
case PH and TH denote the spherical coordinate of a distance 
observer. Since this calculation uses the induced loop 
currents (EP and ET) , a backscattering call must preceed 
this calculation. The outputs returned to the calling 
program consist of the scattered electric field components 
and echo areas. 

EPP(I) and ETT (I) denote the phi-polarized and 
theta-polarized far-zone fields of dipole mode I with unit 
terminal current. In a backscattering situation, the 
excitation voltages EP (I) and ET (I) are obtained by 
multiplying EPP and ETT by the constant CJI. Then calls are 
made to SQROT which stores the solution (the induced loop 
currents) in EP(I) and ET(I) . RITE is called for the branch 



46 






currents CG (J) , and GDISS is called for the time-average 
power dissipated in the imperfectly conducting wire and the 
lumped loads. This power is denoted PDISS and TDISS for 
phi-polarized and theta-polarized incident waves, 
respectively. 

In scattering problems, the incident plane wave has unit 
electric field intensity at the origin. GGG denotes the 
time-average power density of the incident wave at the 
origin. ACSP and ACST denote the absorption cross sections 
for the phi and theta polarizations. 

PIN and TIN denote the time^average power input to the 
wire structure, delivered by the equivalent voltage 
generators VP and VT at the terminals. PIN and TIN apply for 
the phi and theta polarizations, respectivity. The 
time-average power input is regarded as the sum of the 
time-average power dissipated and the time-average power 
radiated or scattered by the wire. ECSP and ECST denote the 
extinction cross sections and SCSP and SCST denote the 
scattering cross sections. 

The distance field is calculated in the DO LOOP ending 
with statement 7 for scattering situations, and in the DO 
LOOP ending with statement 9 for the antenna situation. 

The radar cross sections (echo areas) SPPM, SPTM, STPM, 
and STTM , are defined as 



2 2ar 

o = limit 4nr e S / S 
r ~oD s i 



where S and S denote the time-average power densities in 
s i 

the scattered and incident fields evaluated at the origin. 



For an antenna, the following definition is employed for 
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the power gains: 



2 2ar 

G (8,4>) = limit 4nr e S(r,Q,4>) / P 

p r — ~ i 

where P , GG, denote the time-average power input and 
i 

S(r,8,<|>) is the time-average power density in the radiated 
field. GPP and GTT denote the power gains associated with 
the phi-polarized and the theta-polarized components of the 
field, respectively. 

The use of the variables JFLAG and KFLAG are described 
in subroutine SGANT. 

CALLED BY: MAIN 

CALLS TO: GDISS 

GFF 



RITE 



SQROT 
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SUBROUTINE OF FLO ( I A, IB, INC. INM. l WR, I 1 , I 2* l 3, 1 12 « MD, N* ND,NM, AN* ACS 

3. ETA, GAM, ERR, I GRD) 

COMPLEX ERR 

COMPLEX CJI, ETI , ET2,EPl|EP2,EPPS.ETTS,EPTS.ETPStZS,VP,VT 
COMPLEX CT 1 1 , C J I 1 1 , EPC 1 1 , E T I D . E PP( 1 1 , E T T ( I) , ZLO( 1 1 



complex eta, gam 

0 1 MENS I ON I A ( 1 ) 

DIMENSION 0(1), 

OATA PI ,TP/3. 14159,6. 28318/ 
-4.*PI/(ET A*&AM J 



Ji 1 1 * EPt 1 J.ETl IJ.EPPUl 



13(11* NO( II* HD ( 1 NH*4 1 



CJI 

GGG 

THR 

CTH 

STH 

PHR 

CPH 

SPH 



REAL( I ./ETA) 
.01 745 33*T H 
COS( THR) 

SIN( THR) 

.01 74533*PH 
COS( PHR ) 

S1N( PHR ) 



DO 1 I * I , N 
E TT I I ) « ( . 0 , . 0 ) 

EPPU) » (.0,.0) 



DO 3 K- l » NM 
KA • IAIK) 

K8 • 1 8 ( K ) 

NGRD • IGRD 
IF (K.LE.NM/2) IGRD— 1 



CALL GFF (X(KA) , Y( KA ) . Z I K A I . X< *B ) , Y ( KB ) t Z(KB) .0 ( K ) , CGDOO , SGD (K I ,C 
ITH.STH,CPH,SPH,GAM,ETA,ET1 ,ET2*EPI,EP2, IGRD,§ftR) 

IGRD * NGRD 
NDK - ND( K ) 



DO 3 11*1, NDK 
1 - MD( K » III 
FI « 1 
IF (KB 
IF (KB 



EQ. 12(D) 
EQ.I1( I[ 



i"!ll : mi'll:, F l: £tl 

GO TO 3 

IF (KA.EQ. 13(11) 



GO 

F 1 1 



TO 

1 . 



teem - 



ETT 

CON 



111 

TIN 



INUE 



1 

2 

3 

4 

5 

6 
7 

a 

9 

10 

11 

12 

13 

14 

15 

16 
17 

1 $ 

20 

n 

n 

25 

B 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 
99 
40 

il 

43 

44 

45 

46 

47 
46 



EPPS ■ ( .0,.0) 

ETTS - (.0,.0) 

IF ( INC.EQ.O) GO TO 8 
IF ( INC .E0.2) GO TO 6 

DO 4 I * l * N 
E T I I I « ETT I I ) ♦C J I 
E P ( I I - EPP(1)»CJI 

CALL SCROT (C*EP,0,I12,N) 

(12 * 2 

CALL SORDT ( C , E T . 0. 1 1 2 , N ) 

IF (IWR.GT.O) WRITE (6,10) PH, TH 
IF (IWR.GT.O) WRITE (6,11) 

(Ritl 

. j G ,L “ 

IF ( IWR.GT.O) WRltE 
CALL RITE (1A, 18 



CALL RITE ( IA, 18. INM, IWR.tl .12, l 3 ,M0, NO, NM. EP.CG, IGRO) 
CALL GO (SS ( AM,c6,CMM,D f PDl $ ,GAM , NM , SGO * Z L& , Z $) 



L / , r L 

_ ( 6 , 12 ) 

u all kiil i i a, i o . 1 NM, IWR, t 1 • 12 , I 3, MD.ND, NM, ET .CG, IGRD) 
CALL GOISS ( AM,CG,CMM,D,TDI S,GAM,NM,$GD,ZLO,ZS) 

AC SP * PO ( S/ GGG 
AC S T • TO ( S/GGG 
P ( N - .0 
TIN * .0 

DO 5 1*1, N 
VP « CJ I * EPP ( ( ) 

VT « C J I*ETT( I I 

P (N * P (N+REAL (VP*CON JGIEP( I ) ) ) 

T (N - T (N*REAL(VT*CONJG(ET( I j j j 

ECSP * P ( N/GGG 
ECST « T 1 N/GGG 



SCSP 
SC ST 



ECSP-ACSP 
EC ST- ACS T 



6 EPTS ■ (.0,.0) 

E TPS - (.0,.0) 

DO 7 I-I.N 

EPPS - EPPS»EP( I ) *EPP( I 
EPTS « EPTS*EP I I • F T T ( I 
ETTS • ETTS^ETi I J • E TT I ( 

7 E TPS • E TPS ♦£ T ( 1 )*EPP( ( 

S PPM ■ 2.*TP*(CA8S(EPPS)4*2 
SPTM • 2.MP* CA8S( EPTS 
S TPM ■ 2.*TPMCA8S(ETPsI**2 
STTM • 2.MP*(CA8S( E T T S ) • *2 
RETURN 

8 DO 9 I-I.N 

ETTS - E T T S*C J ( I )*ETT( ( ) 

9 EPPS ■ E PPS*C J ( 1 ) • fc PP ( I ) 

APP * C A8 S ( E PPS ) 

ATT - CABS(ETTS) 

GPP « 4 . ♦ P I • APP* APP*GGG/GG 
G TT - 4 . 4P I *ATT* ATT*GGG/GG 
RETURN 

10 FORMAT ( 1 OX* * BRANCH CURRENT S AS SOCI ATED WITH PL ANE-WAVE SCATTERING 
1 FOR THE INCIDENT ANGLES, PHI-’, F5.lt* * N 2rJ? I J ?c A ’ E I * ^ 

11 FORMAT ( 44X , * CURRENT S INDUCED 8Y THE PHI POLARUEO WAVE* ) 

12 FORMAT (44X, ’CURRENTS INDUCED BY THE THETA POLARIZED WAVE ) 

END 
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63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

?* 

93 
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95 

96 

97 

98 

99 
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GGMM 



PURPOSE: to calacate the mutual impedance between two 

filamentary monopoles with sinusoidal current distribution. 

METHOD: As stated in subroutine SGANT, the mutual 

impedance of coupled dipoles may be expressed as sum of four 
monopole-raonopole impedances. This subroutine calculates the 
mutual impedance with closed-form expressions in terms of 
exponential integrals. 



For skew monopoles it can be shown that the 
monopole-monopole mutual impedance is given by: 



t -z z 

i+j n mm 

z = (-1) B [ e (F -e G + e G) 
ij jl 12 22 



-t -z z 

n mm 

— e (F — e G + e 0)1 
j2 11 21 



where m= 2/i, n = 2/j and 



16 n sinh d sinh d 
1 2 



The functions F are defined by: 

ik 



qz cos V 
i 

F =2 sinh d e E (R + gz cos IJ - qt) 

ik * i i 



where 

being 



k 

q = (-1) , 

considered. 



d and d are the lengths of the monopoles 

1 2 



The functions G are defined as follows: 

ik 
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G = E (R 

ik 2 


+ qz 

2 


+ g ' t 


- jq" 


) + E (R + qz 
2 2 


+ q ' t 


♦ jq") 


- E (r ^ 


* qz 

1 


+ g't 


- jq" 


) - E (R + qz 
1 1 


+ q't 


+ jq") 


i 

where q = (-1) , 


q* = 


k 

(-D , 


and 


q" = qb + q'c 


with 




b = c cos 9 and 


c = d/sin 9 . 


The 


angle 9 is the 


angle 


formed 



by the apparent intersection of the two monopoles. This 
will be discussed later in detail. 



ik 

In the above equation for G , t denotes the position of 

an observation point somewhere on monopole 2. R and R are 

1 2 



the distances from the endpoints of monopole 1 to this 
observation point. Finally, the E functions are defined as 
follows: 



E (a + jg") 



jKq" 

e 



a 




dw 



where a and q" are real quantities with dimensions of 

length, a is a function of t, a = a (t ), a = a(t ) and 

1 12 2 

/ = jw yu&. The integral above is evaluated by subroutine 
EXP J. 




PARALLEL 



SKEW 



To explain the input data for GGMM, refer to the above 
figure. If the monopoles are parallel, then the new 
coordinate system is defined such that the new z axis is 
parallel to the monopoles. The coordinate origin may be 
selected arbitrarily. Si and S2 denote the z coordinates of 
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the endpoints of the test monopole, T1 and T2 are the 
coordinates of the endpoints of the expansion monopole, and 
D is the perpendicular distance (displacement) between the 
monopoles. The mutual impedance of parallel monopoles is 
calculates in the last part of GGMM below statement 5. 

For skew monopoles, let the test monopole s lie in the 
xy plane and the expansion monopole t in the plane z = D. (D 
is the perpendicular distance between the parallel planes.) 
If the monopoles are viewed along a line of sight parallel 
with the z axis, the extended axes of the two monopoles will 
appear to intersect at a point on the xy plane. Let s 
measure the distance along the axis of the test monopole 
with the origin at the apparent intersection. Si and S2 
denote the s coordinates of the endpoints of the test 
monopole. Similarily, let t measure the distance along the 
axis of the expansion monopole with the origin at the 
apparent intersection. T1 and T2 denote the t coordinates 

of the endpoints of the expansion monopole. Let s snd T be 

unit vectors parallel with the positive s and t axes, 

respectively. Then CPSI = s*T = cos (I. The monopole lengths 

are d and d . 
s t 

The output data from GGMM are the impedances P11, P12, 
P21, and P22. In defining these impedances, the reference 
direction is from Si to S2 for the current on monopole s, 
and from T1 to T2 for the current on monopole t. In the 
impedance P , the first subscript is 1 or 2 if the test 

ij 

dipole has terminals at Si or S2 on monopole s. The second 

subscript is 1 or 2 if the expansion dipole has terminals at 

T1 or T2 on monopole t. The monopole lengths d and d are 

s t 

assumed positive in defining the input data CGDS, SGD1 and 
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SGD2 



For parallel monopoles, CPSI = 1 or -1. Si, S2, T1, and 
T2 are cartesian coordinates for parallel monopoles and 
spherical coordinates for skew monopoles. For skew 
monopoles, the radial coordinates Si, S2, Tl, and T2 tend to 
infinity as the angle V tends to zero or n. Therefore, if 
the monopoles are within 4.5° of being parallel, they are 
approximated by parallel dipoles. 

CALLED BY: GGS 



SGANT 



CALLS TO 



EXP J 



SUBROUTINE GGMH C S 1 , S2 , T 1 , T 2, 0, CGOS . SGOI . SGD2 ,CPS 1 . E TA, GAM , P 11 ,P 1 2 
1 ,P2I,P22> 

DOUBLE PRECISION Rl ,R2 ,DPQ, S 1 S . TS 1 . T S2 . S T l , ST2, CO . BO .C PSS, SK . TL 1 1 T 



2 

3 

4 

5 

6 
7 
B 
9 

10 




LUnHLt A t be 1 c * c I i bMU I f ur ( 

DATA PI /3. 14159/ 

D SO * D*D 
SGDS * SGD I 

IF (S2.LT.SU S3DS - -SGDl 
SCOT * SGD 2 





DO 1 K«l,2 
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R 2 * OS QR T <DPQ*SI$*T$2-ST2I 
EK « E8 

00 3 K * 1 • 2 
FK * ( - 1 |**K 
SK * FK*SD I 
EL • EC 

DO 2 L * 1 1 2 
FL * < - 1 I * *L 
EKL * EK*EL 
XX » FK * 80 + FL*CO 
TL1 « FL*TD1 
TL2 = FL*T02 
R R 1 = R 1 ♦SKfTLl 
RR2 = K2*SK*TL2 

CALL EXP J <GAM*CMPLX<RR1,-XX),GAH*CMPLX(RR2.-XX|,EXPA) 
CALL EXPJ <GAM*CMPLX<RRI,XXI ,GAM*CMPLX<RR2,XX) .EXP0I 
E 4 K *L I = E4K,LI*Fl*4EXPA*EKL*EXP8/EKLI 
EL * I./EC 



EK 



1. /E8 



ZO » SDI*DPSI 
ZC * ZO 

EGZ1 - CEXP<GAM*ZC) 

RRI * R 1 ♦ Z 0- T 01 
RR2 • R2*ZD-TD2 

CALL EXPJ IGAMARR1 ,GA M*RR 2 , EXPB I 
RRI - Rl-ZDUDI 
RR 2 = R2-Z0UD2 

CALL EXPJ < GAM* R R I ,GAM*RR2,EXPA) 

F < I * I I * 2.*SG0S*EXPA/£GZ t 
F < I » 2 I • 2.*SGDS*EXP8*EGZ l 

4 SI - S2 

C ST « ETA/<16.*PI*SGDS*SGDTI 

P A 1 • CST*4<F<I.1)*E<2*2I*E$2-E<1*2I/ES2I*ET2M-F<1*2)-E<2*1I*ES2* e 

1 Pll , i , Csf*U-F?{,ll-E<2,2)»ES2^E<i*2l/ES2l*ETl^<Fll.2l*E<2*ll*ES2- ? 

lEU.il/ES2l/ETll 

P21 - CST*< 4-FI2.1 1-E(2,2I*ES1*E< l,2l/ESn*ET2*«F«2,2)*EI2*l I*ES1- 
1E41,1I/ES1I/ET2I 

E?1 1 ,/ESII/eIi I 1 2 f 2l*ESl-E< 1,2 1 /ESI 1*E Tl* 4 -F ( 2 ♦ 2 1 -E 1 2 * 11 *ES1 ♦ 
return 

5 IF < CPS 1 •LT« 0* I GO TO 6 
TA • T 1 
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51 

52 

53 
5* 

55 

56 

57 

58 

59 

60 
61 
62 
63 
6* 

65 

66 
67 
66 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

7 

8 

89 

90 
9 1 

92 

93 

94 

95 

96 



T 8 * T 2 
GO TO 7 

6 TA = - T 1 
T 0 - -T 2 
SGOT * - SGOT 

7 SI - SI 

00 9 1=1,2 
TJ * TA 

DO 8 J= 1 , 2 
Z1J • TJ-SI 
R ■ SQRT4DSQ*ZIJ*ZIJI 
W * R*Z I J 

IF (Z1J.LT.0.I W * DSQ/4R-ZIJI 
V * R-Z < J 

IF 4Z1J.GT.D.) V * DSQ/4R*ZIJ) 

IF U.EQ.J) VI • V 

IF < J.EQ. II W1 = W 

EGZ < I ,JI * CE XP 4 G AM*Z I J I 

8 TJ • T8 

CALL EXPJ 4 GAM* VI ,GAM*V,GP< I II 
CALL EXPJ 4 GAM* Ml , GAM*W.GN4 I I I 

9 SI - S2 

PU - CST*<GH<2^*EGZ< ? 2,2?*GP<2)/EGZ<2,2l-CGOS*CGN< 1 |*EGZ< 1*21 *GP( 1 
l P12 G i < cifiI-GM<2»*EGZ <2 , 1 1 -GP< 2 ) /EGZ < 2, 1 1 *CGDS*< GM< 1 1 *EGZ 1 1 , 1 1 *GP< 
’ 1 P2i e * Z Csf*iGMU)*EGZ< 1.2)*GP( 1 )/EGZ( 1»2)-CG0S*< GH ( 2 1 *EGZ <2,2)*GP<2 

1 P22 G ^ ( c4t*I-GH<1 I* EGZ <1*1 1-GPU I/EGZ4 1, 1 I *CGOS*< GM 4 2 )*EGZ C 2 , 1)*GP< 
1 2 I /EGZ ( 2*1)1) 

RETURN 

END 



97 

98 
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100 
101 
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104 
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107 

108 
109 
1 10 

111 

113 

114 

115 
1 16 
1 17 
118 

119 

120 
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123 
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GGS 



PURPOSE: to calculate the mutual impedances between 

two filamentary monopoles with sinusoidal current 
distributions. 

METHOD: The monopole-monopole mutual impedance as 

defined by SGANT is calculated using the equations defined 
in subroutine GNF. The endpoints of the axial test monopole 
s are (XA,YA,ZA) and (XB,YB,ZB) , and the endpoints of the 
expansion monopole t are (X1,Y1,Z1) and (X2,Y2,Z2). DS and 
DT denote the lengths of monopoles s and t , respectively, 
CAS, CBS and CGS are the direction cosines of monopole s, 
and CA, CB and CG are the direction cosines of monopole t. 

The effects of ground for vertical co-linear monopoles 
are applied in a slightly different manner than mentioned 
previously. As with self impedance calculations, the test 
monopole and the expansion raonopole are laterally displaced 
by the wire radius. This lateral displacement is used to 
determine the angle of incident. This technique is applied 
at statement 8. 

If INT = 0, GGS calls GGMM for the closed form impedance 
calculations. Otherwise GGS calculates the mutual impedance 
via Simpson’ s-rule integration with the following number of 
sample points: IP = INT + 1. If the monopoles are parallel 

with small displacement, GGS calls GGMM to avoid the 
difficulties of numerical integration. 

Since the point (X,Y,Z) of subroutine GNF lies on the 
expansion monopole t, T is the integration variable and is 
measured from (X1,Y1,Z1). Cl is the current at T for the 
mode with terminals at (X1,Y1,Z1), and C2 is the current at 
T for the mode with terminals at (X2,Y2,Z2) . C denotes the 
Simpson ' s-rule weighting coefficient. 
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Below statement 7, GGS performs some analytic geometry 
in preparation for calling GGMM. The remainder of this 
section is concerned with this preparation. 



Let s denote a unit vector in the direction from 

(XA,YA,ZA) toward (XB,YB,ZB). Also let f denote a unit 

vector from (X1,Y1,Z1) toward (X2,Y2,Z2). Then s.t = 

cos 0 = CC where 6 is the angle formed by the axes of the 

two monopoles. Let monopole s lie in one plane P and 

s 

monopole t lie in another parallel plane P^. CAD, CBD and 

CGD are the direction cosines of the unit vector 3 = f x s / 
sin 0 which is perpendicular to both planes. To obtain the 
distance DK between the two planes, a vector If is 

constructed from (XA,YA,ZA) to (X1,Y1,Z1) and take DK = 

E .3. 

i i 

A line is constructed from (X1,Y1,Z1) to the test 
monopole, such that the line is perpendicular to the test 
monopole. SZ denotes the s coordinate of the intersection of 
this line with the test monopole, and the cartesian 
coordinates of this intersection are XZ, YZ, and ZZ. The 

direction cosines of s x 3 are CAP, CBP, and CGP. 

From the point (X1,Y1,Z1) in plane P^, a line is 
constructed perpendicular to the point (XP 1 , YP1,ZP 1 ) in the 
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plane P . This line is 
s 

Let 1 represent a vector 

Pi denotes U (s x 3) . Si 
GGMM. 

CALLED BY: SGANT 

CALLS TO: GGMM 



parallel with 3 and 
from (XZ,YZ,ZZ) to 
and T 1 are defined 



has length DK 
(XP1 , YP 1 r ZP1 ) 
in subroutine 



SUBROUT INE GGS I XA .YA , 2A.X6, YB.2 B ,X1 .Y 1 . 21 , X2, Y2 , 22, AN.OS, CGOS, SGO 
1 S.DT.SGOT, I NT ,ETA,GAM,Pll , Pl2, P2 l , P2$ . E ftR, IGRD) 

COMPLEX EXI,EYl,EX2,EY2,En,EZ2 

COMPLEX P1I,P12,P2I,P22«EJA,EJB*EJ1«EJ2,ETA,GAM,CI,C2»CST 
COMPLEX E GO* CGDS • S GOS » SGOT ,ER1 .ER2.ET1.ET2 

complex err 

COMPLEX EE.EXX.EYY 
COMPLEX PP.PX.PY.PZ 

COMPLEX RRl.RR2.RR3 ,RR4 ,RH1 , RV 1 , RH2 ,RV2 , RH3 ,RV3, RH4, RVA 
DATA F P/12.56637/ 

CA * IX2-XD/0T 
CB • IY2-Y11/DT 
CG « 1Z2-ZD/DT 
CAS » i XB-XAl/DS 
CBS ■ IYB-YAJ/0S 
CGS ■ i ZB-ZAI/DS 

I (Cg!lE?!o 03KAND?(^I.CE*.003I.ANO*( IGRD. GT. OH GO TO 1 
IF lABSICCI.GT. ,997) GO TO 6 
l SZ - »XI-XA)*CASMy1-VA|*CBS*IZWAI*CG$ 

IF I INT.LE.O) GO TO 7 
INS « 2*1 INT/2 ) 

IF I INS.LT.2I INS ■ 2 
IP « INS* 1 
OELT ■ OT/INS 
T ■ .0 

OSZ • CC*0ELT 

Pil - 1 .0, .0) 

P12 - 1 .0, .0) 

P2I • < .0,.0I 
P22 - i .0,.0I 
AMS » A M* AH 
SGN - -1. 

00 5 IN ■ I , IP 
’ SZ-DS 

XXZ ■ Xl*T*CA-XA-SZ*CAS 
YYZ • Yl*T*CB-YA-$Z*CBS 
III ■ Z l*T*CG-ZA-SZ*CG$ 

RS - XXZ**2*YYZ**2*ZZZ**2 
R1 • SQRT IRS*ZZl**2» 

EJA - CEXPI -G AH*R 1 ) 

EJl - EJA/Rl 

R2 • SgRTIRS*ZZ2**2) 

EJB • CEXPI -GAH* R2 ) 

EJ2 * EJB/R2 



0001 

0002 

0003 

0004 

0005 
5Q06 



0008 

0009 

8oll 

0012 

0013 

0014 

0015 
00 1 o 
0017 

ooie 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 
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ERl 
ER2 
FAC 
IF ( 

Ell 
ET2 
IF ( 

RV1 
RH1 
RV2 
RH2 
IF ( 

XG1 
YG1 
2G1 
XG2 
YG2 
ZG2 
RG1 
RG2 
TTI 
TT2 
CTHI 
SSTH1 

CTH2 - COS ( TT2 1 
SSTH 2 - SIN(TT21*SIN(TT21 
CSORTI ERR-SSTH 1 I 
(CTHl-RRll /(CTHURR1I 
-(ERR4CTHl-RRl|/(ERR*CTHl4RRll 
CSQRT( ERR-SSTH21 
(CTH2-RR2)/(CTH2*RR2l 
-( ERR*CTH2-RR2)7(ERR*CTH2*RR2l 
S0RT( (XB-XAlMXB-XAl ♦ ( YB-YA )*( YB- YA ) ) 
0 
0 

IF (RG.LT.l.E-32) GO TO 3 

'PH ' ‘ 

PH 

EXX • ET1 ♦CAS 
ETl*CBS 

C EXX*SPH*EYY*CPH) *(RH1-RV1 ) 
EXX*RVUEE*SPH 
EYY*RV1-66*CPH 
-EU*RV1*CGS 

-EX1*CAS-EYI*CBS*EZ1*CGS 
ET2*CAS 
ET2*CBS 

( EXX*SPH-EYY*CPH ) • (RH2-RV21 



E JAMGDS* ZZ 1 *EJ l*CGOS-ZZ2*E J2 
- E JB*SGDS»ZZ2*EJ2*CGOS-ZZt*EJl 

RSIgT.AHS) FAC • (CA*XXZ4CB*YYZ*CGMZZl/RS 
CC*( EJ2-EJ14CGD$1*FAC*ERI 
CC*( EJI-EJ2*CGOSUfAC*ER2 
GRO.LT. 01 GO TO 4 
(- 1..01 
(-1. tO) 
l-l.*oJ 
C-l.fOl 

GRO.EQ.ll GO TO 2 
XlfT*CA-XA 
YlfT*CB-YA 
Z1M*CG-ZA 
XUT*CA-XB 
Y 1 ♦T*CB- YB 
2 UT*CG-ZB 

SQRT( XGlMGUYGl^YGl) 
SQRT(XG2*XG2*YG2*YG2) 

A T AN ( RG1 / ZG1 I 
ATAN(RG2/ZG2) 

> COS ( TT 1 1 

S l N( TTl 1 ♦ SI N( T T1 1 



RR| 
RH1 
R VI 
RR2 
RH2 
RV2 
2 RG - 
CPH 
SPH 



EYY 

EE 

EX 1 

EY1 

EZI 

ETI 

EXX 

EYY 

EE ' 



0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 
005B 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 
0067 
006B 

0069 

0070 

S°o?i 

0073 

88?t 

0076 

0077 
007B 
0079 
OOBO 
GOBI 
00B2 
00B3 
00B4 
00B5 
00B6 
00B7 
00B8 
00B9 
0090 
009! 

0092 

0093 

0094 
009 5 
0096 



EX2 - EXX*RV2*EE»SPH 


. 0097 


EY2 ■ EYY*RV2-EE4CPH 


009B 


EZ2 » -ET24CGS*RV2 


0099 


ET2 — EX2*CAS-EY2*CBS*EZ2*CGS 


0100 


C ■ 3.f SGN 


0101 


IF ( IN.E0.1.0R.IN.EQ.IP1 C*l. 


0102 


EGO « CEXP<GAHMDT-T)) 


0103 


Cl - C* (EGD-1./EGDI/2. 


0104 


EGO - CEXP(GAM4T) 


0105 


C 2 - C*( EG 0-1. /E 60 1/2. 


0106 


Pll - P1UET1*C1 


0107 


P 12 * P12»ET14C2 


010B 


P21 • P2UET2*C1 


0109 


P22 • P22*ET2*C2 


0110 


T - TfOELT 


0111 


SZ * SZ40SZ 


0112 



5 SGN - -SGN 

CST - -ETA*0ELT/(3.*FP*SG0S*SG0T) 

Pil • C ST*P1 1 
P 12 • CST*P12 
P 21 • CST 4P21 
P22 * CST*P22 
RETURN 

6 SZI - < Xl-XAi*CASMYl-YA!*CBS*IZl-ZA!*CGS 

0R1 - 5QRT((X1-XA-SZ1*CASI**2*IY1-YA-SZI*CBS) **24 ( Z1*Z A-SZI*CGS)** 

1 112 » SZ1*DT*CC 

0R2 - SQRTU X2-XA-SZ2*CAS)**2*( V2-YA-SZ2*CBS)**2MZ2-ZA-SZ24CGSl*4 
121 

000 • (DR1+DR2I/2. 

IF (000.GT.20.4AH.ANO.INT.GT.01 go to i 
IF (OOO.LT. AH) 000 - AH 

CALL GGHH («0*OS*SZ1 t SZ 2 9 000 tCGOS « SGOS * SGOT * 1. «ETA «GAM«P 1 1«P!2«P21 



1 * P22 1 

IF (IGRO.LE.l) RETURN 
IF llGRO.GT.il 



GO TO 8 



7 SS « 

C AO 



SQRT(l.-CCACC) 

(CGS*CB-CBS4CG)/S 



CBO - (CA$4CG-CGS*CA|/| 
CGO - (CBS*CA-CAS*CB)/S 



OK ■ (X1-XAJ4CA04(Y1-YA)*CB04( Zl-ZAl^CGO 
OK • ABS(OK) 



IF (OK 

xz 

YZ • YA*SZ*CBS 



LT.AH1 OK 
XA* 5Z*CAS 



AH 



0113 

8tis 

0116 

0117 

0118 

°oHo 9 

0121 

0122 

0123 

0124 

0125 

0126 
0127 
012B 
0129 

oil? 

0132 

0133 

0134 

0135 

0136 

0137 

mi 

0140 

0141 

mi 

0144 
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ZZ - ZA*SZ*CGS 

XP1 - X l-DK*CAD 

YPl ■ Y1-DK*C3D 

ZPl « Z l -DK*CGD 

CAP * CBS* CGO-CGS*CBD 

C BP « CGS*CAO-CAS*CGO 

CGP ■ CAS*CBO-CBS*CAO 

PI ■ CAP*<XPI-XZ)*CBP*(YP1-YZ1*CGP*(ZPI-ZZ1 
T l ■ Pl/SS 
SI - T1*CC-$Z 

CALL GGMM ( S l , S 1 *OS , Tl , Tl *OT * OK ,CGOS , SGOS t SCOT ,CC * ETA ,GAM , PU tP12 . 
1P21«P22) 

RETURN 

5 AHS - A M* AM 

RG * (Xl-XA)*<Xl-XA)* <Yl-VA)*( Yl-YA) 

IF (RG.LT.AMSJ RG » AMS 

OG * SORT ( < Z 1-Z A )*<Zl-ZA)fRG) 

C PH - ABS< Zl-ZAI/DG 
SSPH * RG/(OG*OGI 
RRI ■ C SORT! ERR- SSPH) 

RV1 - -<ERR*CPH-RRI)/<ERR*CPH*RR1) 

Pll •-Pil*RVi 

RG - <Xl-XBl*<Xl-XB)MYl-YB)*m-YBJ 

IF (RG.LT.AMS) RG - AMS 

OG - SQRTU Zl-ZBT*(Z1-ZB)*RG) 

C PH « ABS< Zl-ZBI/OG 
SSPH - RG/ I OG*OG I 
RRI ■ CSQRTI ERR-SSPH) 

RV1 ■ -<ERR*CPH-RRII/<ERR*CPH*RR1) 

P12 «-Pl2*RVl 

RG « <X2-XA)*(X2-XA|f ( Y2-YAI*! Y2-YA) 

IF (RG.LT.AMS) RG - AMS 

OG - SQRTUZ2-ZA)*(Z2-ZA)*RG) 

CPH « ABS1Z2-ZAI/OG 
SSPH * RG/(OG*OG) 

RRI » CSQRT(ERR-SSPH) 

RVI • -(ERR*CPH-RR1)/IERR*CPH*RR1) 

P21 — P2l*RVl 

RG - (X2-XBIM X2-XB)f I Y2-YB l*( Y2-YBI 

IF (RG.LT.AMS) RG - AMS 

OG » SORT ( ( Z2-Z B) * < Z2 - ZB) *RG) 



CPH • A 8$( Z2-Z 8 jf /OG 

■ RG/ (OG*d6| 

C$0 Rt I ERR- SSPH) 



SSPH 

Rvt • -(gRR*CPH-RRi)/(ERR*CPH*RRi) 
P22 — P22*RV1 
RETURN 
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014B 

0149 

0150 

0151 

0152 

0153 

0154 

0155 
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0157 

0158 
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0163 
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0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 
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m\ 

01B3 

01B4 
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GNF 



PURPOSE: to calculate the near-zone electric field of 

a sinusoidal electric monopole. 



z E 




1 



METHOD: An electric line source is located on the z 

axis with endpoints at z and z as shown in the above 

12 

figure. Let the electric monopole have the following current 
distribution : 



KD 



I sinh /(d - 1) + I sinh yi 

, 1 2 

sinh /d 



where and 1^ are the endpoint currents, y is the complex 



propagation constant of the medium, d = - z^ is the 

source length. The cylindrical components of the field are 
E (0) = 0 and 

77 -ys -ys 

E(e) = [ (I e 1 - I e 2) sinh yd 

4u6 sinh yd » 2 

-ys 

+ (I cosh yd - I ) e 1 cos 0 
12 1 

-y r 

+ (I cosh yd - I ) e 2 cos q t 
2 1 2 
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-/R 

I cosh yd) e 
2 R 

2 

-Y R 

I cosh yd) e 1 ] 

i R 

i 

where is the intrinsic impedance of the medium and where 
(6 , <J) , z) denote the cylindrical coordinates in a coordinate 

system centered at the endpoint of z . 

These expressions exclude the field contributions from 
the point changes at the endpoints of the line source, since 
these charges disappear when two monopoles are connected to 
form a dipole. 

Let the coordinate s measure distance along the test 
monopole with the origin at (XA,YA,ZA) . From any point 
X,Y,Z, a line is constructed perpendicular to the monopole. 
SZ denotes the s coordinate of the intersection of this line 
with the monopole. The length of the line is the radial 

2 

coordinate 6, and RS denote 6 . R1 and R2 are the distances 
from (X A , YA , ZA) and (XB,YB,ZB) to the point (X,Y,Z) . 

In the statements above statement 1, the above equations 
are solved; and after statement 1, the cartesian components 

(E ,E ,E ) of the field are determined. If a ground plane is 
x y z 

present (IGRD>0) the reflection coefficients are applied to 
the cartesian components before returning to the calling 
program. 

CALLED BY: GNFLD 

CALLS TO: NONE 



E (z) = [ (I 

4n sin h yd ‘ 



+ (I - 
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SUBROUTINE GNF ( X A, Y A , l A. X B, Y B . ZB * X , Y , l , AH, OS * CGOS , SGOS * E TA, GAM , EX 
1 1 ,EY1,EZ1 , EX2, E Y 2, EZ 2, I GftO. ERR I 
CCMPLEX ERR t RVl , RH1 , RV 2 , RH2 , RR 1 , R R2 * E 6 

COMPLEX E JA, E JB,E J1 *EJ2 .ER1 * E R2 . E SI * ES2 * SGO S , GAM . CST , CGO S * ET A 
COMPLEX EX1 ,EY1 ,EZi,EX2,EY2#EZ2 
DATA PI/3.14159/ 

CAS * < X6-XA ) /OS 
CBS * (YB-YAl/DS 
CGS = (ZB-ZAJ/OS 

SZ = ( X-XAMCASMY-YA) *CBS* U-ZA)*CGS 

III * SZ 

III *= SZ-OS 

XXZ * X-XA-SZ*CAS 

YYZ « Y-YA“SZ*CBS 

III * Z-ZA-SZ*CGS 

RS * XXZ**2*YYZ**2*ZZZ**2 

Rl * SQRT ( RS ♦ Z Z l * *2 1 

EJA * CEXP(-GAM*R1) 

EJ1 * EJA/R1 

R2 * SORT ( RS ♦ Z Z 2 * * 2 ) 

EJ6 * C E X P ( -GAM * R2 ) 

EJ2 * EJ8/R2 

E J2-E J 1 *CGDS 
E J i-E J 2*CGD S 
< . 0 , . 0 ) 

( .0, .0) 

AM* AM 

(RS.LT.AMS1 GO TO 1 



ESI 

ES2 

ER1 

ER2 

AMS 

IF 

CTH1 

CTH2 

ER1 

ER2 

CST 

E X 1 

E Y 1 

EZ1 

EX2 

EY2 

EZ2 



ZZ 1/Rl 
<* ZZ2/R2 

( EJA*SG0S*EJA*CG0S*CTH1“EJB*CTH2I/RS 
( -E JB*SGOS*E J8*CG0S*CTH2-E JA*CTHi )/RS 
ETA/(4.*P1*SG0S1 
CST*(ES1*CAS*ER1*XXZ) 
CST*(ES1*CBS*ER1*YYZ| 

CST*(ESl*CGS*ER i*ZZZl 
CST*(ES2*CAS*ER2*XXZ| 



CST* (ES2*C8S*ER2*YYZ 
CST*(ES2*CGS*ER2*ZZZl 
IF CIGR0.LE.01 RETURN 
RVl * (-1. ,01 
RH1 « C-l.,01 
RV2 * I - 1 . . 0 > 

RH2 • C - 1 . * 0 1 
IF (IGRO.EG.il GO TO 2 

K1 * SQRH<XA-X)MXA-X1MYA-Y)*<YA-Y|) 
R 2 * SQRT( < XB-X) MXB-X) MYB-Y )*< YB-YII 
TH1 * ATAN(Rl/( Z A- Z I ) 



1 

2 

3 

4 

5 

6 
7 
B 
9 

10 

11 

12 

13 

14 

15 

16 
17 
IB 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 
30, 

31 

32 

33 

34 

35 

36 

37 
3B 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 



TH2 
R R 1 
RR2 
RVl 
R HI 
R V2 
RH2 
RG * 
CPH 
SPH 
IF ( 
CPH 
SPH 
EE * 

set 

u i . 

EX2 

Ilf 

RETURN 
END 



ATAN(R2/(ZB-ZI I 
CSQRT(E*R-SIN(TH1)*SIN<THI J I 
CSORT(ERR-SlN(TH2l*SINiTH2 )) 
-<ERR*COS(THll-RRl)/(ERR*CO$ITHl J*RRl) 
(C0S(TH1 )-RRlJ/(COS( thii*rrIi 
-< ERR*COS< TH2l-RR2)/< ERR* cost TH2)*RR2) 
(COS(TH2)-RR2)/(COS(TH2l*RR2) 

SGRT< (XA-XBIMXA-XBIM YA-VB )*(YA-YBI) 



0 
0 

RG. LT . 



l.E-321 GO TO 3 



(XB-XAI/RG 
(YB-Y 



. YA 1 /RG 

( EX1 * SPH-E Y 1 *CPH ) *(RHl-RV i ) 
-EXI*RV1*EE*SPH 
-EYl*RVl-EE*CPH 
EZ1*(-RV1 1 

(EX2*$PH-EV2*CPH)*(RH2-RV2) 

-EX2*RV2*EE*$PH 



_Y2*RV2-EE*CPH 

EZ2*(-RV2l 
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GNFLD 



PURPOSE: to calculate the near-zone electric field 

intensity at a given point. 

METHOD: This subroutine calls GNF for the near-zone 

field of each wire segment, and sums over all segments to 
obtain the near-zone field of the wire antenna. FI is used 
in a manner similiar to FI of subroutine SGANT. CJ (I) is the 
loop currents calculated by subroutine GANT 1 . 



The use of the variables JFLAG and KFLAG are described 
in subroutine SGANT. 



CALLED BY: 



MAIN 



CALLS TO: 



GFF 



SUBROUTINE GNFLO < IA.IB, INM, 1 1 . I 2 ,13, MO, N, NO, 
l,CJ,D,X,Y »Z*XP» YP « ZP , EX « £Y , EZ, IGRO,£RR ) 

COMPLEX EX,EY,EZ,EX1,EvI,E*1,cX2,£Y2,E£2,ETA, 

COMPLEX ERR 

complex cjm.cGom , sgoiii 

DIMENSION lAill, I Bi I )« till). 12(1), 13(11, DU), X(l), Y(l), Hi 



N,ND,NM,AM,CGO,SGD,ETA,GAM 
GAM 



EX 

EV 

EZ 



DIMENSION MO ( INM ,41, ND(1I 
OATA PI .TP/3. 14159,6.28318/ 

( .0, .0) 

t.0,.0) 

( . 0 ,. 0 ) 

DO 2 K* l » NM 
KA * I A ( K ) 

KB - I B ( K ) 

NGRO « IGRO 

IF (K.LE.NM/2) IGRD— l 

i I . Y C K A » • Z(KA) .XIKBI.YiKBI.ZIKBl.XP.YP, ZP , AN, D( Kl ,CGD 
l iSil^GOI K * , ET A, bAM ,EXl ,EYl,EZl,EX2,EY2,EZ2,IGRO,£RR) 

IGRD « NGRO 
NDK « ND( K ) 

DO 2 1 1 *1 .NDK 

IF TkbI EQ.I2C II ) GO TO l 
IF (KB. EQ. IK III FI — 1 - 
EX « EX*FI*EXI*CJ 

EY ■ EYtF|*EYl*CJ 

| l • EZ*FI*EZl*CJ 

GO TO 2 

L IF (KA.EQ.mm FI — 1 
* §*tE!*§*2?cj( h 
EZ 

l CONTINUE 






0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 



RETURN 

END 



001 
001 
00 

881S 
0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 
003 5 
0036 
00 37 
003B 

0039 

0040 
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LEFT 



PURPOSE: to determine position (location) of the left 

paren symbol on the input data card. 

METHOD: The character search begins in the column 

passed to the subroutine. On returning to the calling 
program the argument passed is the column following the left 
paren symbol. 

CALLED BY: READ 

CALLS TO: NONE 



c 



c 



SUBROUT 1NE LEFT (Nl 
COMMON /A/ A < 80 I 
OATA PLEFT/M'/ 

K • N 

00 1 I •K* 00 
N ■ IM 

IF (All I « EQ.PLEFT ) GO TO 2 

1 CONTINUE 

N • l 

2 PE TURN 
END 



I 

. 2 
3 

A 

5 

6 
7 
0 
9 

10 

11 

13 
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LINECK 



PURPOSE: to insert grid characters on the polar plot. 

METHOD: The peroid character (ISM (2)) is inserted in 

the proper position in the statements above statement 4. In 
the statements after statement 4, the grid numbers labels 
are inserted on the horizontal axis. 

CALLED BY: POLPLOT 



CALLS TO: 



NONE 



SUBROUTINE LINECK <X,Y) 

THIS SU8R0UT I NE INSURES ALL GRIO CHARACTORS LIE ON THE POLAR GRID 

COMMON ISYMtL iNE 
INTEGER Y 

DIMENSION I SYMi 14 I • LINE(130) 

IF (Y.EQ.O) GO TO 5 
K = 0 

IF (X.LT.10.0) GO TO 5 

SET UP AREAS OF "PERIOD" POLAR GRIO POINT CHARACTERS 

I * I NT i XI 
I * I ABSi 1) 

Z « A8SIXI 

IF ((Z-D.GT.0.5) M»1 

1 IF (Z. LT. 10.0. DR. Z.GT. 111.0) GO TO 2 
L1NEIU - ISYM(2» 

L1NEI60) - ISYMI3) 

L INE I 62 ) • ISYMiil 
K « KM 

IF (K.EQ.2I GO TO 2 
1 * 122-1 
GO TO 1 

2 L INE (61 ) - I SYM( 2 ) 

IF (Y.NE.O) GO TO 5 

3 00 4 K= 1I« 111 
LINE(KI - I S Y M ( 2 ) 

4 CONTINUE 

FILL IN GRIP NJM8ER LABELS ON HORIZONTAL AXIS 



LlNEill ) 
L 1 NE ( 2 0 ( 
L INE t 2 1 1 
L 1 NE ( 2 2 » 
L I NE ( 301 
L 1 NE ( 3 I I 
LINE! 32 I 
L INE( 40 ) 
L I NE ( 4 1 ) 
L1NE<42) 
LINE(50I 
L1NEC51 ) 
LINEI52) 



I S YM ( 7 ) 

I SYMi 10) 
1SYM( 51 
l SYMi 1 1 I 
I SYMi 9 ) 

I SYMi 5» 

I SYMi 1 1 | 
ISYMi 8) 

1 S YM ( 5 I 
ISYMi 11) 
ISYMi 7) 

I S YM ( 5 ) 
ISYMI 11) 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
2i 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 



L INE ( 61 ) - 
L l NE ( 70 I - 
L INEi 71 I - 
L I NE ( 72 I « 
L I NF ( 80 I * 
L I NE ( 8 1 ) « 
L I NE ( 82 i - 
LINE (90) ■ 
L INEI 91 ) « 
L INEi 92 I « 
LINEi 100) 

L I NE ( 1 0 1 { 
LINE* 102) 
LINEI 1 1 1 1 
CONTINUE 
RETURN 
END 



ISYMi 1 ) 
ISYMi 7) 

I SYMi 5 » 
ISYMi II) 

I SYMi 81 
ISYMi 5 ) 
ISYMi 11) 

1 S YM ( 9 ) 

1 S YM ( 5 ) 
ISYMi 11 ) 

■ ISYMi 10 ) 

■ ISYMi 5 I 

» ISYMI 1 Li 

• ISYMI 7 I 



49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 
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NUMB 



PURPOSE: to place degree numbers on the polar plot. 

METHOD: The current line which is being printed is 

passed to the subroutine in the calling argument. If this 
line contains degree numbers, these numbers are placed in 
the correct position by the IF statements. 

CALLED BY: PTPLOT 



CALLS TO: NONE 



SU 8 RGUT IN 6 NUM 8 ( Y) 

THIS SUBROUTINE PUTS DEGREE NUMBERS ON POLAR GRID 

COMMON ISYM.LINE 
INTEGER Y 

0 1 MENS I ON ( SYM ( l 4 1 9 LINEU30) 

IF (Y.NE.37) GO TO I 



LINE I 33) 
1 1 NE ( 34 ) 
L I NE ( 3 5 ) 
LINE! 87 ) 
L INEI 68 ) 
LINE! 89) 



I S YM ( 7 ) 

I SYM ( 8) 
ISYMl 6 ) 

I SYM I 6 ) 

I SYM< 12) 
I S YM ( 6 ) 



IF IY.NE.2I) GO TO 2 
L INEI I ? ) * I S YM ( 7 ) 
LINE! 13) - ISYM(il) 

L INEI 14 ) - I S YM ( 6 ) 
LINE! 108) = ISYMI 6 ) 
LINE! 109) = I SYMI 9) 
LINE(IIO) = I SYM ( 6 ) 

IF (Y.NE.O) GO TO 3 
LINE! 7 ) = I SYM C 7 I 
LINEI 8 ) * I SYM ( 13) 
LINEI9) * ISYMI 6 ) 

L (NE( I I 3) = I SYMI 6 ) 

L I NE ( 114 ) = ISYMl 6 ) 

L INE III 5) = I S Y M ( 6 ) 

I F ( Y.NE.-21 ) GO TO 4 
LIMEU2) • ISYMl 8 ) 
LINE! 13) - ISYMl 7) 
l(NE<!4! • I SYM ( 6 ) 

L I NE ( l 0 6 ) - ISYM19) 
LINE! 109) « I SYM ( 9 ) 
LINE! 110) • ISYMl 6 ) 

IF (Y.NE.-37) GO TO 5 



LINE 33 
LINE! 34 \ 
L INF ( 35 ) 
L INEI 87 ) 
LINEI 88 ) 
L INEI 89) 
CONTINUE 
RETURN 
ENO 



I SYM I I 
I S YM ( 6 
ISYMl 9 
ISYMI 6 
I SYM i 6 
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9 
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22 
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35 

36 

11 

39 

40 

41 

42 

43 

44 
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NUMBER 



PURPOSE: to convert alpha-numeric numbers to floating 

or fixed point numbers. 

METHOD: After initially determining the sign of the 

number, the DO LOOP ending at statement 6 scans each 
character beginning at N1. The DO LOOP ending at statement 
3 terminates the outer DO LOOP if the character being 
compared is not an alpha-numeric number. The DO LOOP ending 
at statement 5 converts the alpha-numeric number to an 
actual number. Below statement 7, the multiplier correction 
is applied to the floating point number before returning to 
the calling program. 

CALLED BY: READ 

CALLS TO: NONE 



SUBROUTINE NUMBER (Ni,N2,X,IX) 
COMMON /A/ A10O) 

D I ME NS ION B ( 10 ) 

DATA 

DATA AMNUS, PLUS, POINT/ ,*♦•,*. '/ 
DATA AK,AM,AU/«K* ,*«* , •U*/ 

N = Ni 
NSIGN * 0 
l I * -1 

IX * 0 
ISET * 0 

IF (A(N).EQ.PLUS) N*NM 
IF (A(N) .NE.AMNUS) GO TO l 
NSIGN = l 
N * N«d 

DO 6 I * N , 80 

IF (Ad >.NE. POINT ) GO TO 2 
ISET * 1 
GO TO 6 

IF ( I SET . EQ« l ) II - ll*l 
DO 3 K= l, 10 

IF ( A ( I ) .EQ.8(K) ) GO TO 4 
CONT INUE 

GO TO 7 

DO 5 K=l,lO 
K K = K— 1 

IF (Ad J .EQ.B(K) ) NUM8*KK 
CONT INUE 

I X * NUM8M0MX 
N 2 = Id 
CONTINUE 

IF (NSIGN. EO.U IX - -IX 
Y * IX 

IF ( I I. IT. 01 II « 0 

X * Y/< 10**1 I ) 

IF (A(N2) .EQ. POINT! N2*N2M 
IF A N2I.E0.AK) X ■ X*L000. 

IF j A(N2) .EQ.AM 
IF (A(N2).E0.AU 
I F ( ( A ( N 2 ) . EQ. AK 
Nl * N2 
RETURN 



X • X*0.00I 
X - X*0. 000001 

0R.(A(N2) .EQ.AM) .OR. ( A ( N2 ) .EQ.AU ) ) N2-N2M 
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28 

29 
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31 
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POLPET 



PURPOSE: to control the plotting of the polar plot. 

METHOD: This subroutine is the main subroutine in the 

polar plot package and is responsible for calling the 
various subroutines of the package. 

The scale factor, S, must be changed according to the 
printer characteristics. The scale factor in this subroutine 
is set for ten, 10, characters per inch for the abscissa and 
eight, 8, characters per inch for the ordinate axis. 
Therefore S = 10. /8. 

After initializing DATAX, DATAY, and X, the input data, 
Y, is scanned to determine the normalizing factor. If this 
normalizing factor is less than 1.E-32, an error statement 
is printed and the plotting is abortted. 

In the DO LOOP ending with statement 8, each line of 
the polar plot is printed after a call is made to PTPLOT to 
establish the ploar grid information. The variable, DIM, is 
used to as a scaling factor for the polar plot. The value of 
1.0 will cause all of input data to be plotted, however, if 
only the values less than one-half of the normalizing factor 
are of interest, then DIM can be set to .5. This will 
enlarge of the center of the polar plot. 

CALLED BY: MAIN 

CALLS TO: PTPLOT 

SART 
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ooooooo o noono oooo non oo oooo oo o oooo oooo oooo 



SUBROUTINE POLPRT ( NAME * Y ) 

COMMON I SYM. LINE 

DIMENSION X i 360 ) * Y(360), DATAXI360), 0ATAY1360), LINEI130), ISYMI 

DIMENSION TI TLA (21* TITL2I2) 

DATA T I TLA/ • PH I f # • THET* A 

N ■ 360 

DIM - I .0 

NST * 1 

KST - I 



S IS SCALE FACTOR OF PRINTER* 

ABSCISSA CHAR . PER INCH / ORDINATE CHAR. 

S ■ 10.0/8.0 

ZERO DATAX AND DATAY 



PER INCH 



DO I I A * I » N 
0 * IA-1 

DATA X( I A) * 0.0 
DATA Y( I A I = 0.0 
l X(IA) = 0*3.1415927/180.0 



FACTOR IS THE NORMALIZING DIVISOR 
FACTOR = YII) 

DO 2 I A *2 * N 

2 IF (FACTOR.LT.YI I A ) ) FACT0R*Y ( I A ) 



IF (NAME.EQ.I) T1TU-TITLA( I) 

IF (NAME.EQ.2) T I TL1 = TI TLA( 2 I 

IF ( (NAME.EQ.3) .0R.(NAME.EQ.4).0R.(NAME.EQ.7).0R.(NAME.EQ.8 )) TITL 
1 2 ( 1 ) *T I TL A ( I ) 

IF ( ( NAME.EQ.5I.0R. ( NAME. EQ. 6 ) . OR . ( NAME . EO . 9 ) .OR. ( N AME. EO. 10 ) ) TIT 
I L 2( l)=TITLA( 2) 

IF ( (NAME.EQ.3 ). OR. (NAME. EO. 5 ).0R.( NAME. EO. 7). OR. (NAME. EO. 9 )) TITL 
l 2 ( 2) *T I TLAI1I 

IF ( (NAME.E3.4).0R. (NAME.E0.6) .OR. I NAME. EO . 8) .OR. I NAME . EO* 10 I ) TIT 
1L2C2) *T I TLA( 2 J 

IF (FACTOR. GT.l.E-32) GO TO 3 
IF (NAME. LE. 21 WRITE (6,9) TITLl 
IF ( NAME. GE . 3 ) WRITE (6,10) TITL2 
RETURN 
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29 
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NCRMALI ZE DATA TO ONE 



3 DO 4 IA*1.N 

4 Y( IA) « Y( I A) /FACTOR 



IF (NAME.LE.2) WRITE (6,11) TITLl, FACTOR TfTt , cirTno 

IF ( ( NAME. GE . 3) . AND. ( NAME.LE.6) ) WRITE (6,13) TITL2, FACTOR 
IF (NAME.GE.7) WRITE (6,12) TI TL2 ,FACTOR 
FILL DATAX AND DATAY ARRAY FROM X AND Y ARRAY 



DO 5 I A = l * N 

DATA X( IA) * Y ( I A ) *COS ( X ( IA) 
5 DATA Y4 I A ) * Y( I A 1 * SIN ( X ( I A ) ) 



SORT DATA BY ORDINATE MAGNITUDE 
CALL SART ( D A TAX , D AT AY ,N ) 

DATAX AND DATAY ARE SORTED BY DESENDING MAGNITUDE ON THE DATAY VAL 
SET UP FOR PLOTTING POLAR GRID WITH DATA 



C 



DO 8 I Y Y*1 ,81 
CALL PTPLOT ( IYY,S) 

LINE IS RETURNED WITH POLAR GRID INFORMATION 

SET UP # Y • BIN SIZE UPPER ANO LOWER LIMITS 
ULL IS THE LOWER BIN LIMIT 
UL IS THE UPPER BIN LIMIT 

BIN • DIM/80.0 

ULL - DIM-(2*IYY-1 )*BIN 

UL * UL L* 2* 8 1 N 



CYCLE THROUGH DATA TO FIND WHICH ONES FALL IN f Y* BINS 

IF (NST.GT.N) GO TO 7 
DO 6 J J *NS T , N 
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OOP 



8 

C 

C 

c 

c 



c 



c 

c 



IF (OATAYl JJl.LT.ULL) GO TO T 
KST - JJ 

AMAG - SQRT(DATAXIJJ)+DATAX(JJ)*DATAY1JJ)*0ATAYIJJ)) 

CHECK THAT MAGNITUDE IS NOT OVER OIM 
IF ( AMAG. GT. OIM) GO TO 6 

OK IS THE FINAL LINE POSITION FOR THE *** 

OK - DA TAXI JJ)AS*40.0/DIM*6i .0 
IF I0K.LT.10.0) GO TO 6 
K - I NT I OK) 

K - lABS(K) 

OK « ABS(OK) 

IF HOK-K) .GT.0.5 ) K«KH 

IF IOK.LT • 10.0.0R.0K.GT .1 11 .0) GO TO 6 

LlNE(K) * ISYM<4) 

6 CONTINUE 

7 CONTINUE 

NS T * K S T ♦ 1 

PRINT OUT ONE LINE OF PLOT 

WR ( TE (6,14) LINE 

8 CONTINUE 

RETURN 

•9 FORMAT ( 10X,1A4.' COMPONENT OF THE ELECTRIC FIELD IS LESS'/iOX, 
l 'THAN l.E-64, THEREFORE THIS FIELD WAS NOT ' / 10X , • PLOT T EO. EXEC 
2 UT ION WILL CONTINUE AS NORMAL.'//) 

10 FORMAT ( l OX » ' THE MAXIMUM VALUE OF THE BISTATIC PATTERN FOR •/ 

1 10X, 1A4.'-* ,1A4, • ( INCIOE NT-SCATTERED) IS LESS THAN */ 

2 10X, • l.E-30.) POLAR PLOT NOT CALLEO.'///) 

11 FORMAT I'1',IA4.' ELECTRIC FIELD ANTENNA PATTERN FOR SPECIFIED PLA 
INE.' • 9X • ' NORMAL I Z I NG FACTOR- '.E10.5) 

12 FORMAT «• 1BISTATIC SCATTERING PATTERN F OR • . 1 A4, • - • , 1 A4 , • | INC 1 DENT- 
1SCATTERE0) POLARIZATION. '.9X, 'NORMALIZING FACTOR*',Ei0.5l 

13 FORMAT ('1BACKSCATTERING PATTERN FOR • , l A 4 , • , l A4 , • I INC 1 OENT- SC AT T 
1ERED) POLARIZATION. ',9X, 'NORMALIZING FAC TOR - • , E 10 . 5 ) 

14 FORMAT 1 1 X , 1 30 A1 ) 

ENO 
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PTPLOT 



PUBPOSE: to establish the grid information for the 

polar plot. 

METHOD: In the DO LOOP ending at statement 1 the 

alpha-numeric characters are transferred to ISYN in order 
to pass via COMMON to other subroutines. In the statements 
following statement 2, the equations for the plotted 
concentric circles are established. Below statement 7 the 
grid marks on the 090-270 axis are inserted. 

CALLED BY: POLPRT 

CALLS TO: LINECK 

NUMB 
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ono o ooooo noon 



SUBROUTINE PTPLDT (IYY,SI 

THIS SUBROUTINE SETS UP POLAR GRID INFORMATION 
COMMON ISYH.UNE 

0IMENS10N LlNE( 130) * ISYMU4), 1SYNU4I 

OAT A 1SYN/ IHf , 1H. ,1H , 1H*,1H/, 1 HO . I H l , 1 H2 , 1 H3 1 1H4 * l H5, 1 H6 * l H0 , 1H9/ 
INTEGER Y , Y Y » W 

SET UP ISYM FROM 1 SYN FOR COMMON 



00 I K = 1 » 1 4 
ISYM(K) * ISYN(K) 
1 CONTINUE 



CLEAR LINE ANO SET TO BLANK 



00 2 1=1*130 

2 LINE! I I « ISYMI3) 

Y * 41-IYY 
IF (Y.EQ.O) GO TO 7 

SET UP EQUATIONS FOR CONCENTRIC CIRCLES 
YY ■ Y * Y 

Z = (YY*2.5/2)*S 
X = 61 .OfSQRT( 2500.0-Z) 
call LINECK <X,Y| 

IF tY.GT.32.0R. Y.LT.-32I GO TO 3 
X = 61 .0* SQRT( 1600.0-ZI 
CALL LINECK <X,Y) 

3 IF (Y.GT.24.0R. Y.LT.-24) GO TO 4 
X « 61 • 0* SQR T (900 • 0“Z ) 

CALL LINECK ( X * Y ) 

4 IF (Y.GT. L6.0R.Y.LT.-16) GO TO 5 
X - 61 .0*$QRT<400.0-Z F 

CALL LINECK ( X » Y ) 

5 IF (Y.GT.8.OR.Y.LT.-0I GO TO 6 
X - 61 .OfSQRTt IOO-Z) 

CALL LINECK (X,Y) 

SET UP EQUATIONS FOR MULTIPLES OF 30 DEGREES 

6 X * 61 .0*1.732051 *Y*S 
CALL LINECK (X,YI 

X - 61 .0*Y*S/1. 732051 
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CALL LINECK ( X * Y ) 
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PUT IN POLAR PLOT NUMBER LABELS 


51 








52 




CALL NUMB <Y) 




53 




W - I ABSt Y) 




54 








55 




fill in polar plot 


AT 000, 090, 180, AND 270 


56 








57 




IF (W.NE.40) GO TO 


0 


50 




L I NE ( 5 5 ) - I S Y M ( 2 ) 




59 




LINEI57) « I SYH ( 2 ) 




60 




L 1 NE ( 59 ) - I SYM ( 2 I 




61 




L I NE ( 63 ) « I SY M ( 2 ) 




62 




L I NE ( 65 ) = ISYMI2) 




63 




L I NE ( 6 7 ) * ISYM(2) 




64 


0 


IF (W.NE.32) GO TO 


9 


65 




LINEC56) ■ ISY^I2) 




66 




L INE (58 ) ■ I SYM< 2 ) 




67 




LINE! 60 ) * I SY H ( 2 ) 




60 




L I NE ( 62 ) * I SY M ( 2 1 
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L I NE ( 6 4 ) * I SYM ( 2 ) 




70 




L I NE ( 66 ) = I S YM ( 2 ) 




71 


9 


IF (W.NE.24) GO TO 


10 


72 




L I NE ( 5 7 ) ■ 1 SYM( 2 ) 




73 




L 1 NE ( 59 ) • I SYM ( 2 1 
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L I NE ( 60 ) - I S Y M ( 2 ) 




75 




L I NE ( 62 ) ■ I SY M( 2 ) 




76 




L INE ( 63 ) - I SYM ( 2 | 




77 




L I NE ( 65 ) - I SY M( 2 1 




70 


10 


IF (M.NE.I6) GO TO 


11 


79 




L I NE ( 58 ) = I S YM( 2 | 




00 




L I NE ( 60 ) » I SY M ( 2 1 




01 




L I NE ( 62 ) - I SY M( 2 | 




02 




LINEI64) - ISYM(2) 




03 


11 


IF (W.NE.O0) GO TO 


12 


04 




LINE159) * 1SYMI2) 




05 




L I NE ( 63 1 = 1 SYM ( 2 i 




06 


12 


CONTINUE 




07 




RETURN 




00 




END 




09 
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READ 



PURPOSE: to interpret and translate the input data 

cards. 

METHOD: The program utilizes free format for the data 

cards, that is, the program uses character recognition to 
determine which parameters are being read. In the IF 
statements containing A(1), A (2), A (3), and A (4), the first 
four characters on the data card are compared to the first 
four letters of the key words. This will determine the type 
of parameters that card contains. The other IF statements 
determine which parameters are being read. 

Subroutine BLNK is called to remove the blank spaces on 
the parameter cards. Subroutines EQUAL and LEFT are called 
to determine the position of the equal character and the 
left paren, respectively. Subroutine NUMBER is called to 
convert the alpha-numeric characters to numbers, either 
fixed or floating point. This numerical value is assigned 
to the parameter just determined. 

A detailed explanation of the data cards is found in 
appendix II titled "USERS MANUAL". 

CALLED BY: BLNK 

EQUAL 

LEFT 

NUMBER 
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non 



SUBROUTINE RE AO ( I A . I 8 , 1 B 1 S C . 1 C A RD. 1 G A I N , 1 G RO, 1 NE A R, IN T , I SC A T , 1 WR , 
l lFLAG,KFLAG,KGEN,LOAD,LZD,MSG,NBAP,N8IP ,NFFP,NGEN,NM f NP, A8AP . A8AT. 
2AFFP.AFFT.ABIP, AB IT , AM.BM, C MM . E R 2 , E R 3 . £ R 4 . F MC . HGT . PHAF . PHAI.PHIF 
3HI 1 , PHSF ,PHSI , THAF.THAI.THIF, TH1 l.THSF, T hM , stG2»$lG3»5lG4,TD2,Tu3 
4,VOLT,X»XNP,Y,YNP,Z»ZLLD*ZNP*STEPl 
CCMMON /A/ A ( 8 0 I 



j I ON I A ( l ) • I B I 1 * X(l). Y(l), Zlll, KGEN(l), KFLAG(l) 

> I ON XNP(l), YNP(l). ZNP( l I . LZDU) 

)A • A8 , ACfAO, AE , A F«A g,AH,A I .AKtAL,AMA,AN,AO,AP,AQ,AR • A S » A T * A U . 
f * A* , * B * » * C * * , D , **E** , F** , G , **H , ,*1 , , , K*, , L , ,'M , f'N*, , 0 , ,'P , f 



U 



COMPLEX VOLTU ) # ZILD< l I 

0 1 ME NS I ON 1 A ( 1 ) . IBID, XU). YU), Z ( 1 I , KGEN( D • KFL AG ( 1 I 

DIMENSION * * 

DATA AA, ‘ 

I AW , AX / * 

2 , Q , »*R* » , S , , , t , DU* ,*W , , , X*/ 

DATA BL ANK, COMMA, MINUS, PLEFT.POI NT, R IGHT, SLANT/* 

I,*)',*/*/ 

RAD * 5 7.295 779 
INT * 4 
IBISC = -I 
I GAIN * -I 
INEAR = -I 
ISCAT = -1 
IWR = -I 

IF UFLAG.E0.6I GO TO 2 
IF (MSG.NE.O) GO TO 4 

1 READ ( 5.76,ENO=72) A 

2 IF ((A(i).NE.AC) .OR. ( A< 2) .NE.8L ANK ).OR. ( A ( 3 I . NE. 8L ANK 1 . OR . ( A ( 4 ) . NE 
l . BLANK) ) GO TO 3 

WRITE (6,74) A 
GO TO I 

3 WRITE (6,75) 

GO TO 5 r ' 

4 READ (5, 76.ENDU2) A 

5 ICARD = ICAROM 
WRITE (6,77) ICARD, A 

IF U MSG.NE.O) . AND. (( A(l ).EQ.AE) .AND • ( A ( 2) . EG. AN) «AND.( A( 3) .EQ.AO) 

D) GO TO 7D 

I F ( (MSG.NE.O). AND. ((AID .EQ . AS ) .AND. ( A ( 2 ).EQ.AT ) . AND. ( A ( 3 ) . EQ. AO ) 
I . AND. (A(4) .EQ.AP) ) ) GO TO 69 

IF ((AU).EQ.AC) .AND. ( A ( 2 ) . EQ • 8L ANK ) • AND • ( A ( 3 I • EQ. 8L ANK ) • ANO . ( A ( 4 ) 
1.EG.8LANK)) GO TO 73 
IF (MSG.GT.O) GO TO 4 
CALL 8LNK (A) 

N * 4 

INSULAT ION 

IF ( ( A(l) .NE.AI ) .OR. ( A( 2) .NE.AN) .OR.( A( 3) .NE.AS).OR. ( A (4) .NE. AU) ) 

1 GO TO ID 
KFLAGI20) - 1 
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21 
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CALL LEFT (N) 

6 IF { (A(N).NE.AR) .OR.(A(N+l).NE.AA).OR.(A(N+2).NE.AD).OR.(A(N+3) .NE 
• •Ail) Gu TO 7 

K F L AG ( 4 ) - I 

CALL EQUAL (N) 

CALL NUMBER (N,N2,X1,1X) 

8M ■ XI 

IF ( A(N2) .EQ. RIGHT ) GO TO 4 
IF (A(N2).NE. SLANT) GO TO 71 
N ■ N2* I 
GO TO 6 

7 IF ((A(N).NE.AD).0R.(A(NM).NE.AI).0R.(A(N»2).NE.AE).0R.(A(N+3)«NE 

1 . AL) ) GO TO 8 
KFLAGI6) * 1 
CALL EQUAL IN) 

CALL NUM8ER (N,N2 ,Xi, IX) 

ER2 * XI 

IF ( A(N2) .EQ. RIGHT) GO TO 4 
IF (A(N2).NE. SLANT) GO TO 71 
N * N2U 
GO TO 6 

8 IF ((A(N).NE.AC) .OR. (A(N*1) .NE.A0).0R.(A(N*2). NE.AN) .OR . (A(N+3).NE 
1 .AD) ) GO TO 9 

KFL AG ( 5 ) * l 
CALL EQUAL (N) 

CALL NUM8ER (N,N2,X1,IX) 

SIG2 * XI 

IF ( A(N2» .EQ. RIGHT ) GO TO 4 
IF (A(N2).NE. SLANT) GO TO 71 
N * N2U 
GO TO 6 

9 j}A( JO.NE.AL) .OR.CAINUDNE. AQ).OR.< A(N*2i.NE.AS I.OR.( AIN*3).NE 
1 . A S ) ) GO TO 71 

KFLAGI7) - 1 
CALL EQUAL (N) 

CALL NUMBER (N,N2,X1, IX) 

T 02 * XI 

IF (A(N2).EQ. RIGHT ) GO TO 4 
IF (A(N2).NE. SLANT) GO TO 71 
N « N2U 
GO TO 6 
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onno 



( A(2).NE.Al) .0R.( A( 3) .NE.AR) , 
(A(NM) .NE.AAl.OR.I A(N*2) .NE, 
tlX) 



0 R • ( A ( 4 ) .NE • AE ) ) 
AD).OR.4A(N+3).NE 



10 IF ( l At 1 ) .NE .AW) .OR. 

1 GO TO l 3 

CALL LEFT (N) 

1 1 JF ( ( A( N) .NE.AR) .OR. 

1 . AI ) ) GO TO 12 

• KFLAGI2) = 1 
CALL EQUAL (N) 

CALL NUMBER (N«N2«X1 
AH * XI 

IF (A(N2) .EQ. RIGHT) 

IF (AIN2I.NE. SLANT) 

N * N2* 1 
GO TO 11 

12 IF ( (AIN) .NE.AC) *OR.(A(NM) .NE . AO) .OR. ( A ( N* 2 ) .NE . AN ) .0 R. ( A ( N A 3) .N E 
1 • AO ) ) GO TO 71 

KFLAGI3) * 1 
CALL EQUAL (N) 

CALL NUMBER <N,N2fXl, 

C MM * XI 

IF ( A(N 2) • EQ.R IGHT ) 

IF ( A( N2 ) * NE • SLANT ) 

N * N 2 ♦ 1 
GO TO 11 

EXTERNAL MEDIUM 

13 IF ( ( A( l ).NE.AE) .0R.( A (2) .NE.AX) .0R.( A4 3 ) .NE . AT ) .OR. ( A ( 4 ) .NE . AE ) ) 

1 GO TO l 7 

KFLAG(B) « 1 
CALL LEFT (N) 

14 IF (( AIN). NE.AC) .OR. I A(N*1) .NE.AO).OR.4A(N*2) .NE* AN) .OR •(A(NO).NE 
1 .AO) ) GO TO 15 

KFLAG49) * 1 
CALL EQUAL (N) 

CALL NUMBER ( N f N2 * XI % 

S1G3 « XI 

IF (A(N2).E0. RIGHT) 

IF (A4N2) .NE. SLANT) 

N * N 2* 1 
GO TO 14 



GO TO 4 
GO TO 71 



IX) 



GO TO 4 
GO TO 71 



It IX) 

GO TO 4 
GO TO 71 



15 IF ( ( A(N) .NE.AD) .OR 
1 . AL) ) GO TO 16 
KFLAG(IO) * l 
CALL EQUAL (N) 



( A ( N+ 1 ) . NE • A I ).0R« (A4N+2) .NE.AE) «0R.(A(N+3) .NE 
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CALL NUMBER (N,N2,Xl,IX) 

ER 3 * XI 

IF ( A(N2) .EQ. RIGHT ) GO TO 4 
IF (A(N2) .NE. SLANT) GO TO 71 
N « N2M 
GO TO 14 

l6 l! AS) ) A G0 , T0 E 71 U #0R#( A(NM) •NE.AOl.OR.I A4N»2) .NE • AS ) .OR • 4 A 4 N + 3) • NE 
KFLAG4 11) * l 
CALL EQUAL (N) 

CALL NUMBER <N,N2,X1,1X) 

T 03 ■ XI 

IF (A(N2).EQ. RIGHT) GO TO 4 
IF ( A( N2) .NE. SLANT ) GO TO 71 
N * N2M 
GO TO 14 

LOAO 

17 1G0 T0 A 1 £ ,#NE * AL) A(2) -NE.AO) .0R.(A(3)«NE.AA).0R.(A(4).NE«AD)) 

KFLAG414) « l 
GO TO 19 

l8 l l GO ( TO < 22 ?NE * An • 0R - (A( 2)* NE * AMA >* 0R .<A( 3 ).NEfAP , .0R.(A(4).N E .AE) ) 
KF LAG (24) * l 

£all 1 left (N) 

20 CALL NUMBER (N.N2.XI.IX) 

IF (IX.LE.O) GO TO 21 
LZO(I) = IX 
N * N2U 

CALL NUMBER <N,N2,X1,IX) 

RM AG * XI 
N « N2M 

CALL NUMBER (N v N2fXl,IX) 

RDEG * XI 

R RE AL - RMAG*CO$<RDEG/RAD) 

R 1 MAG ■ RMAG*S 1N( ROEG/RAD) 

ZLLD4 I ) « CMPLX4 RREALfRIMAG) 

LOAD - I 

IF (A4N21.EQ. RIGHT) GO TO 4 
IF 4A4N2) .NE. SLANT ) GO to 71 
l - l ♦ l 
N « N2M 
GO TO 20 



19 



145 

146 

147 
14B 

149 

150 

151 
15 2 

153 

154 

155 

156 

157 
15B 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

lit 

174 
1 75 

176 

177 

178 

179 

180 
181 
IB 2 

183 

184 

185 

186 
187 
I 88 
IB9 

190 

191 

192 



75 



21 KFLAG(2*) - -l 
LOAO - *1 

GO TO * 

FREQUENCY 

22 IF (( A( l) .NE.AF) .OR.( AI2) .NE.4R) .OR.(A( 3).NE.AE) .OR.(A(*J .NE.AQ) ) 

IGO TO 23 

KFIAGU) * 1 , 

CALL LEFT <N) 

CALL NUMBER (N,N2,X1, IX) 

FMC - XI 
GO TO * 

PLOT 

23 IF ( ( A( 1) .NE.AP) .OR. ( A ( 2 ) .NE. AL ) .OR. ( A ( 3 ) .N E . AO) . OR. ( A ( * ) .NE . AT ) ) 
IGO TO 31 

KF L AG ( 2 2 ) * l 
CALL LEFT (N) 

2 A IF (<A(N) .NE.AF). OR. (A(N*l).NE.AA).OR.( A (N» 2) .NE .AR) .OR. ( AlN*3) .NE 
l.AF)) GO TO 25 
IGAIN * 1 

NFFP » 1 

GO TO 21 

25 IF ( (A( N) .NE *AB) «OR.(A(N+l) .NE.Al ) • OR. ( A ( Nf 2) .NE.AS) .OR. ( A ( N ♦ 3) .NE 
1 . AT) ) GO TO 26 

1BISC = 1 
NBIP * l 
GO TO 27 

26 IF (( A(N).NE.AB) . OR . ( A (N» 1 ) . NE . AA ) .OR. I A ( N A 2 ) .NE . AC ) .OR . ( A I Nf 3 ) . NE 
1 • AK ) ) GO TO 71 

1 SC AT - l 
NBAP - 1 



27 DO 28 ( “N t 80 
K » 1*1 

IF (A(I ).EQ. SLANT) GO TO 29 

28 CONTINUE 



GO TO 71 
29 N - K 

IF (( A(N) .NE.AT) . OR. ( A ( Nf 1 ) .NE . AH ) .OR. ( A ( Nf 2 ) .NE.AE).0R.(A(N+3).NE 
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206 
207 
20B 

209 

210 
211 
212 
213 
21 * 

215 

216 

217 

218 

219 

220 
221 
222 
223 
22 * 

225 

226 

227 

228 

229 

230 

231 

232 

233 
23* 

235 

236 

237 

238 

239 
2*0 



l . AT) ) GO TO 30 
CALL ECUAL (N) 

CALL NUMBER <N,N2,Xl,lX) 

IF (NFFP.EO.l) A F F T = X 1 
IF (NBtP.EO.O A 3 I T = X l 
IF (NBAP.EQ.l) ABATcXl 
IF ( A(N2) .£0. RIGHT) GO TO * 

(F ( A( N2) .NE. SLANT) GO TO 71 
N « N2*i 

GO TO 2* , 

30 IF ( (A(Ni .NE.AP) .0R.( A(NM ).NE.AH) .0R.( A(N*2).NE.AI )) GO TO 71 
CALL EQUAL <Nl 

CALL NUMBER (N f N2.Xl,lX) 

IF (NFFP.EO.l) AFfP=Xl 
IF (NB(P.EQ.l) ABIP-Xl 
IF (NBAP.EQ.l) ABAP-X1 
IF ( A(N2) .EU.ft (GHT) GO TO * 

IF ( A(N2) .NE. SLANT ) GO TO 71 
N * N2»l 
GO TO 2* 

OUTPUT 

31 IF (( A I D.NE.AO) ,0R.(A(2) .NE.AUI.OR. ( A ( 3 ) .N E . AT ) . 0 R. ( A ( * ) .NE . AP) ) 

IGO TO ** 

K F L AG ( 2 2 ) = 1 
CALL LEFT C N ) 

32 IF ( ( A( N) .NE.AB) *0R.( A(N»1 ) .NE.Al ).0R.( A(N*2) .NE .AS) .0R.(A(N>3) .NE 
1 . AT) ) GO TO 33 

K F L AG ( 18) - 1 
I B < SC = l 
CALL ECUAL (N) 

CALL NUMBER (N,N2»X1» IX) 

PHSI - XI 
N « N2» l 

CALL NUMBER (N,N2»X1»(X) 

PH SF ■ XI 
N - N2M 

CALL NUMBER IN»N2, Xl» IX ) 

THSl • XI 
N - N2M 

CALL NUMBER ( N*N2 » XI * IX) 

THSF - XI 

IF (A(N2) .EQ. RIGHT) GO TO * 

IF (A(N2) .NE. SLANT ) GO TO 71 
N • N2»l 
GO TO 32 



2*1 
2*2 
2*3 
2 ** 
2*5 
2*6 
2*7 
2 *B 
2*9 

250 

251 

252 

253 
25* 

255 

256 

257 

258 

259 

260 
261 
262 
263 
26* 

265 

266 

267 

268 

269 

270 

271 

272 

273 
27* 
275 
2 76 

277 

278 

279 

280 
281 
282 
283 
28* 

285 

286 
287 



76 



ooo non 



33 l !af! ) A G0 J T0 E 3S F * * 0R * u<Nf 11 •NS.AAJ .OR. I AIN* 2) # NE. AR) .OR . I A(N*3 ) .NE 
KFL AG ( 16) • l 
1GAIN = l 
CALL ECUAL < N ) 

CALL NUMBER (N.N2.X1, 1 X9 
PHAI * XI 
N = N2*l 

CALL NUMBER (N.N2.X1, IX) 

PHAF * XI 
N = N2*l 

CALL NUMBER (N,N2 f XI, IX) 

THAI = XI 
N = N2*l 

CALL NUMBER (N,N2,XI,1X) 

I HAF = XI 

IF (AIN2) .EQ. RIGHT) GO TO 4 
IF (AIN2) .NE. SLANT) GO TO 71 

GO TO 32 

3 \! arI ! a go ! To E io NI • OR * U<N * l, * NE * A E).OR.(A(N*2).NE.AA).OR.(A(NO).NE 

KFL AG ( 19) * l 
(NEAR = 2 
CALL EQUAL (N) 

IF (A(N).EQ.PLEFT) GO TO 35 
I NEAR a l 
1 * l 
GO TO 36 



35 DO 37 L • l » 50 
I * L 

N = N*1 

36 CALL ^ NUMBER (N,N2,X1,1X) 

N » N2*l 

CALL NUMBER (N,N2,Xl f (X) 
YNP(l) -XI 
N » N2*l 

CALL NUMBER (N,N2,Xl,lX) 

ZNPI ( ) -XI 

IF ( INEAR.EQ.l) GO TO 39 
INEAR * L ♦ 1 

IF ( A(N2) .EQ. RIGHT ) GO TO 38 



289 

290 

291 

292 

293 

294 

295 
2 96 
29 7 
298 
2 99 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 



N * N2 
37 CONTINUE 



GO TO 71 - 

38 N2 ■ N2*l 

1 NEAR = (NEAR-1 

39 IF (A(N2) .EQ. RIGHT) GO TO 4 
IF (A(N2J .NE. SLANT) GO TO 71 
N * N2* 1 

GO TC 32 f 

40 IF ((A(N).NE.AB) .OR.IAIN+1) .NE . AA ) .OR. ( A ( N* 2) .NE.AC) .OR • I A (N* 3) .NE 
l.AKJ) GO TO 41 

KF L AG ( 17) * 1 
1SCAT = 1 

CALL EQUAL (N) 

CALL NUMBER IN,N2»Xl 9 lX) 

PHI I * XI 
N * N2f 1 

CALL NUMBER (N,N2 f Xl*lX) 

PHIF * XI 
N = N2* I 

CALL NUMBER (N*N2.Xl,lX) 

THI1 « XI 
N * N2* l 

CALL NUMBER (N»N2 v Xlt IX) 

• TH1F * XI 

IF I A I N 2 ) .EQ. RIGHT) GO TO 4 
IF (A(N2).NE. SLANT) GO TO 71 
N ■ N2* 1 
GO TO 32 

41 IF (( AIN) .NE.AC ) .0R.(A(N*1).NE.AU).0R.(A(N*2) .NE.AR).0R.tA(N*3).NE 
l.ARI ) GO TO 43 

KFLAG( 15) • l 

IWR - l 



00 42 K-N.80 

IF (A(K). to. RIGHT) GO TO 4 
n ■ k*i 

IF (AIK). EQ. SLANT ) GO TO 32 
42 CONTINUE 

GO TO 71 
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346 
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348 
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350 

351 

352 

353 

354 

355 

356 
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358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 
3 70 
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372 

373 
3 74 
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3 76 

377 

378 

379 

380 

381 

382 

383 
3 84 
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no 



*3 IF ( ( AIN) .NE.AS) .0R.(A(N*1).NE.AT).0R.( A ( N* 2 ) .NE . A E ) .OR. ( A ( N* 3) .NE 
1 • AP ) ) GO TO 71 
CALL EQUAL (N) 

CALL NUMBER (N,N2,X1, IXI 
STEP = XI 

(F (A(N2).EQ. RIGHT) GO TO 4 
IF ( A ( N2 ) .NE . SLANT ) GO TO 71 
N = N2*I 
GO TO 32 

FEEO POINT 

44 IF KAI U.NE.AF) .OR. ( A( 2) .NE.AE ) .OR. ( A( 3 ) .NE.AE ) .OR. ( A( 4 ) .NE .ADI ) 

1G0 TO 45 
KFLAGU3) * I 

4 5 IF ((All).NE.AG) .OR. < A( 2 ) .NE. AE ) .OR. ( A( 3 ) .NE. AN) .OR. ( A( 4 ) .NE • AE ) ) 

1 GO TO 4 9 
KFLAG123) * 1 

46 NGEN * 0 
CALL LEFT IN) 

47 CALL NUMBER (N,N2,X1,1XI 
NGEN = NGEN* 1 

KGEN 1 NGEN) = IX 

IF ( A(N2) .EQ.R13HTI GO TO 4 

N = N2* l 

CALL NUMBER (N,N2 ,X( , I X I 
VMAG * XI 
N = N2* L 

CALL NUMBER (N,N2,Xl, IX) 

VOEG * XI 

VREAL * VHAG*COS( VOEG/RAOI 
VIMAG = VMAG • S ( N ( VDEG/RAOI 
VOLTINGEN) * CHPLX (VREAL , VI MAG) 

IF ( A(N2) .EQ. RIGHT GO TO 4 
IF < A( N2) .NE. SLANT) GO TO 71 

IF ( (A(N2) .EQ. SLANT) .AND. (A(N2*l) .EQ. BLANK) ) GO TO 48 
N * N 2 ♦ 1 
GO TO 47 

48 RE AO (5,76) A 
ICARO = ICAROH 
WRITE (6,77) ICARO, A 
N * 1 

CALL 8LNK (A) 

GO TO 47 
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386 

387 
380 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 
4(3 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 



OESCP IPT10N 

49 IF ((A( 1).NE«A0).0R.(A(2).NE.AE). OR. (A( 3) • N E. AS ) .0R.(A(4).NE.AC) ) 

1 GO TO 52 
KFLAG112) * 1 
J * 0 

CALL LEFT (N) 

50 CALL NUMBER (N,N2,X1,1X) 

J = JU 

NM = J 
IA(J) * IX 
N = N2*l 

CALL NUMBER ( N , N2 , X 1 , IX) 

I 8 ( J ) = IX 

IF ( A(N2) .EQ. RIGHT) GO TO 4 
IF ( A(N2) .NE. SLANT) GO TO 71 

IF l(A(N2I.EQ. SLANT). ANO . ( A ( N* 1 ) . EQ . BLANK)) GO TO 51 
N = N 2 * 1 
GO TO 50 

51 R E AO (5,76) A 

(caro ■ (Caro* I 

CALL 8LNK (A) 

WRITE (6,77) ICARO, A 
N * I 
GO TO 50 

GEOMETRY 

52 IF (( A( 1J.NE.AG) .0R.( A(2).NE.AEI .OR.(A( 3).NE.A0) •OR. ( A ( 4 ) .NE . AMA) ) 

I GO TO 55 
KF LAG (12) « 1 
JJ = 0 

CALL LEFT (N) 

53 CALL NUMBER (N,N2,XI,IX) 

JJ = JJ*1 

NP = JJ 
X(JJ) * XI 
N * N 2* 1 

CALL NUMBER (N,N2,X1, IX) 

Y1JJ) * XI 
N « N 2* 1 

CALL NUMBER (N,N2,X1,IX) 

Z I J Jl « XI 

IF ?A(N2 ) • EQ. RIGHT ) GO TO 4 
IF (A(N2).NE. SLANT) GO TO 71 

IF ((A(N2) .EQ. SLANT). ANO. ( A (N* 1 ) . EQ. BLANK)) GO TO 54 
• N - N2* 1 
GO TO 53 
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461 

462 

463 

464 

465 
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4 79 
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54 R E AO (5,76) A 

1 CARO * ICAROM 
WRITE (6,77) 1CAR0.A 
CALL 8LNK (A) 

N * l 
GO TO 53 

INTERVAL FOR CALCULATION 

55 IF ( ( A( I) .NE.Al) .OR. (A(2) .NE.AN) .OR.(A( 3) .NE.AT ) .OR. (A (4) .NE.AE) ) 
1 GO TO 56 

KF L AG ( 2 1 ) * 1 
CALL LEFT (N) 

CALL NUM8ER (N,N2,Xl, (X) 

(NT = IX 

IF ( A(N2) .EQ. RIGHT) GO TO * 

GO TO 71 

GROUNO 

56 (F ( ( A ( D.NE.AG) .0R.(A(2) .NE .A R) . OR. ( A( 3 ) .NE. AO ) , 
l GO TO 66 

KF L AG ( 2 5 ) = 1 
K F L AG ( 2 6 ) * 1 
1GRO = 2 
CALL LEFT (N) 

57 IF ((A(N),NE.AP).OR.(A(NM).NE.AE).OR.( A(N*2).NE. 

1 .AF) ) GO TO 58 

IGPO * l 
GO TO 64 

58 IF ((A(N).NE.AGI.0R.(A(N*i).NE.A0).0R.(A(N*2).NE, 

1 • AO ) ) GO TO 59 

ER4 = 30. 

S1G4 = .02 
GO TO 64 

59 IF ((A(N).NE.AP).0R.(A(NM).NE.A0).0R.(A(N*2).NE, 

1 .AR) ) GO TO 60 

ER4 * 4. 

S1G4 = .001 
GO TO 64 

60 IF ((A( N) .NE.AS) .0R.( A(N*1) .NE. AE) .0R.( A(Nf 2) .NE< 

ER4 » 80. 

S1G4 « 4. 

GO TO 64 

61 IF ( (A( N) .NE.AH) .0R.(A(N*ll .NE. AE ) .OR • ( A(N*2) .NE< 

1. AG) ) GO TO 62 

CALL EOUAL (N) 

CALL NUM8ER (N,N2,XI,IX) 



OR. (A(4) .NE.AU) ) 

AR).0R.(A(N*3).NE 
AO) *QR.(A(N+3).NE 

A0).0R.(A(N+3).NE 

AA)) GO TO 61 
Al ) .OR. ( A ( N + 3) .NE 
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490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 
50 5 

506 

507 

508 

509 

510 

511 

512 
5(3 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

III 



HGT * XI 

(F ( A ( N 2 ) • EQ . R 1GH T ) GO TO 4 
IF (A(N2) .NE. SLANT) GO TO 71 
N * N2» 1 
GO TO 57 

62 IF ( (AIN) .NE.AC) .OR.(A(N+l) .NE.A0).0R.( A(N*2). NE.AN) *OR.(A(NO).NE 
l.ADll GO TO 63 

CALL EQUAL (N) 

CALL NUMBER (N,N2,X1,!X) 

SIG4 = XI 

IF ( A(N2) .EQ.R(GHT) GO TO 4 
IF (A(N2).NE. SLANT) GO TO 71 
N * N2U 
GO TO 57 

63 IF ((A(N1 . NE.AO) «0R.(A(N+1) .NE.Al). OR. ( A (N+2) .NE.AE) .0R.(A(N*3).NE 
1. AL) ) GO TO 71 

CALL EQUAL (N) 

CALL NUM8ER (N,N2,X1,IX) 

ER4 = X( 

IF ( A(N2) .EQ.R(GHT) GO TO 4 
IF ( A(N2) .NE. SLANT) GO TO 71 
N * N2M 
GO TO 57 



64 00 65 K S N* 80 

(F (AIK). EQ. RIGHT) GO TO 4 



N * K*l 

(F (A(K).EQ. SLANT) GO TO 57 
65 CONTINUE 



GO TO 71 



66 (F ( ( A( 1) .NE.AS) .0R.( A( 2) .NE.AT) .OR. ( A( 3 ) . N E. AO) . OR . ( A ( 4 ) . NE . AP) ) 
1 GO TO 67 

1FLAG * 2 
RETURN 

67 (F ( ( A ( l) .NE.AC) *0R.( A(2) .NE.AH) .OR . ( A( 3 ) .NE. AA ) . OR. ( A ( 4 ) . NE . AN) ) 
1G0 TO 68 

I FLAG * 3 
RETURN 

68 IF (( A( 1). NE.AE) .0R.( A(2). NE.AN) .OR. (A( 3) .NE.AD) ) GO TO 71 
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530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 
5 70 

571 

572 

573 

574 

III 



(FLAG * 1 

return 

69 (FLAG * 5 
return 

70 (FLAG * 4 
R C TURN 

71 MSG * l 

KF L AG ( 30 ) * (CARO 

72 ^ F UFLAG.NE.5) WR(TE (6,78) 
IFLAG = 5 

RETURN 

73 IFLAG = 6 
(CARO * I CARO- 1 
RETURN 



74 FORMAT (5X,80A1) _ 

75 FORMAT J//7/5X, *0ATA CAROS*//) 

76 FORMAT I 80 A I > 

77 FORMAT (6X , 12 ,2X , 80AI ) 



78 FORMAT ( • 
ENO 



J»sis°END CARO/STOP CARD *C SSING**** 9 I 
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RITE 



PURPOSE: to generate a list of branch currents from 

the input loop currents. 

METHOD: The generation of branch currents is 

accomplished in the DO LOOP ending at statement 2. The 
branch currents are stored in CJ (I) by the latter part of 
the DO LOOP ending at statement 3. If the branch currents 
are reguested for output (IWR positive ) , the DO LOOP ending 
at statement 5 accomplishes this. 

CALLED BY: GANT1 

GFFLD 

CALLS TO: NONE 
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on on 



SUBROUTINE RITE ( l A , I 6 . INM, IWR, I 1 , 12. I 3, MD.ND.NM.C J ,CG, IGRO) 
COMPLEX C J 1 1 1 »CG ( 1 ) tCJA.CJB 

DIMENSION I A ( I ) , I8(l), 11(1)» 12(1). 13(1). HD(INM,4), ND( 1 ) 

AMAX - .0 



DO 3 K- I , NM 
KA * I A ( K ) 

KB * I 8 ( K ) 

C JA * ( .0, .0) 
C JB * ( .0, .0) 
NDK * ND( K ) 



DO 2 I I * l * NDK 
I * MD(K,ll) 

FI * I. 

IF (KB. EQ. 12(1) ) GO TO 1 
IF (KB. EQ. II ( l) ) FI—l. 

C JA = CJA*FI*CJ( I ) 

GO TO 2 

1 IF (KA. EQ. 1 3 1 I)) FI--I. 

C JB * CJB*Fl*CJ( I) 

2 CONTINUE 



CG(K) = CJA 
KK = K* NM 
CG(KK) = C JB 
AC J * CABS(CJA) 

8CJ = CABS ( CJ8 ) 

IF (ACJ.GT.AMAX) AMAX-ACJ 
IF (BCJ.GT .AMAX) AMAX*6CJ 
3 CONTINUE 



IF (IWR.GT.O) GO TO 4 
RETURN 

IF (AMAX.LE.O.) AMAX* 1 • 
WRITE (6,8) 

NMG * NM 

IF IIGRD.GT.O) NMG ■ NM/2 

DO 5 K= 1 * NMG 
CJA * CG( K ) 

KK * K*NM 
C JB * CG(KK) 

CC JA - CABS(CJA) 



1 

2 

3 

4 

5 

6 
7 
B 
9 

10 

11 

12 

13 

14 

15 

IS 

IB 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
3B 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 



CC JB * CABS(CJB) 




49 


AC J * CCJA/AMAX 




50 


BCJ * CCJ8/AMAX 




51 


PA = .0 




52 


PB = .0 




53 


IF (ACJ.GT.O.) PA - 57.2957B*ATAN2(A1MAG(CJA) ,REAL(CJA ) ) 


54 


IF (8CJ.GT.0.) PB ■ 57.29578*ATAN2( AIMAG(CJB) ,REAL(CJ8) ) 


55 


WRITE (6,7) K.IA(K) ,CJA«CCJA. AC J .PA.IB(K) , 


CJB.CCJB.BCJ ,P8 


56 

57 


WRITE (6.6) 




58 

59 


RETURN 




60 

61 


FORMAT ( I HO ) 




62 

63 


FORMAT (2X,I2,2(2X. I2.2X.E1 1.5,1X»E11.5* IX 


• E 1 1 .5 • 1X.E11.5, 1X.F6.1 ) 


64 


1) 




65 


FORMAT (/2(46X. ’NORMALIZED*. 5X)/’ SEG*,2(* 
1GINARY* ,3X,’MAGNI TUDE * , 3X , • MAGN ITUOE • , 3X , * 


NODE* .4X, ’REAL * ,6X , * I MA 


66 


PHASE*) ) 


67 


END 




66 
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SART 



PURPOSE 



to sort data for polar plot 



METHOD 



This subroutine sorts the values of the 



points to be plotted by the polar plot package starting with 

the greatest positive value of y to the greatest negative 

value. In the DO LOOP ending at statement 1 , the value of 

(x ,y ) is interchanged with the value of (x ,y ) if y is 
i i D 3 3 

greater than y . 

i 

CALLED BY: POLPRT 

CALLS TO: NONE 



iUDKUUUINt i AK I IUAI AAiUAI AT*N| 

0 I ME NS I ON DA T AX ( 500 I t DAT AY t 500 ) 

THIS ROUTINE SORTS DATA IN DATAY BY MAGNITUDE 

NN - N-l 

DO 2 I 3 1 1 NN 
NM » IM 

nn i I. him kt 



SUBROUTINE SART I DAT AX .DATAY.N I 



C 

c 

c 



* 

5 

6 
7 



C 



6 



C 

C 



C 



2 CONTINUE 




RETURN 

ENO 
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SGANT 



PURPOSE: to calculate the mutual impedance between 

filamentary monopoles. 



1 

2 



B 



SOURCE 

DIPOLE 




RECEIVING 

DIPOLE 



Z = P11 + P12 + P21 + P22 



METHOD: In the induced emf formulation, the mutual 

impedance of coupled dipoles is 



Z = - 



j I 2 (t) E t (t) dt 



where I (t) denotes the current distribution (normalized to 
unit terminal current) on dipole 2, and E (t) is the field 



of dipole 1 when it transmits with unit terminal current. 
Distance along the axis of dipole 2 is denoted by the 
coordinate t. E may be expressed as the sum of the fields 



from each of the monopoles comprising dipole 1. Furthermore, 
the integral is the sum of the integrations over each of the 
monopoles comprising dipole 2. Thus, the dipole-dipole 
mutual impedance may be expressed as the sum of four 
monopole-monopole impedances. 



It may be convenient to draw the above figure in terms 
of monopoles with the current distribution shown as dotted 
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lines. (The monopole letters remain the same.) 



\ 






\ ~ 
\ 1 


\ P11 / 

J I 


c 


A 


\ PI 2 \ 

l \ 

l \ 


I / 

' / 






r 

■ l 


; P2i : 

/ i 

L-J 


c 


B 


; P 22 > 

V v] 



The surface impedance is calculated just above statement 

2. B01 denotes J / J where J and J are the Bessel 

0 1 0 1 

functions of order zero and one with complex argument, ZARG. 
It is assumed that all the wire segments have the same 
radius, conductivity and surface impedance. 

In the DO LOOP ending with statement 3, SGANT calculates 
the segment lengths D(J). DMIN and DMAX denote the lengths 
of the shortest and longest segments. If the wire radius or 
the segment lengths are clearly beyond the range of 
thin-wire theory, N is set to zero at statement 4 followed 
by RETURN to the main program to abort the calculation. 

At statement 5, the program selects a segment K, and a 
few statements below this it selects another segment L. K is 
a segment of test diple I, and L is a segment of expansion 
mode J. The mutual impedance between segments K and L is 
obtained by calling subroutine GGS or GGMM. In statement 18, 
this impedance is lumped into C (MMM) . The mutual impedance 

Z between dipoles I and J is the sum of four 

ij 

segment-segment impedances. 

The variables IFLAG and JFLAG are used if a ground 
plane is present for the calculation of the mutual impedance 
elements. If IFLAG is equal to JFLAG, the mutual impedance 
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terms will not have the effects of a ground plane since both 
monopoles lie on the same side of the ground interface. If 
the monopoles are on the opposite sides of the interface 
(IFLAG not equal to JFLAG) , the reflection coefficient 
correction must be applied to the mutual impedance elements. 
This same technique is applied in subroutines GNFLD and 
GFFLD . 

In SGANT, segment K has endpoints KA and KB, and segment 
L has endpoints LA and LB. It is convenient to think of KA 
and KB as points 1 and 2 on segment K, and LA and LB as 
points 1 and 2 on L. The four segment-segment impedances 
can be defined as P(IS,JS). The first subscript IS refers 
to the terminal point on segment K, and the second subscript 
JS refers to the terminal point on L. Thus IS=1 or 2 if 
dipole I has its terminal point 12(1) at KA (point 1) or KB 
(point 2) , respectively. Similarly, JS=1 or 2 if mode J has 
its terminal point 12 (J) at LA or LB. The impedances 
P(IS,JS) are defined with the following reference directions 
for current flow: from point 1 toward point 2 on each 
segment. If dipole I has this same reference direction on 
segment K, FI=1; otherwise FI=-1. Similarly FJ=1 or -1 in 
accordance with the reference direction for mode J on 
segment L. In statement 18, P(IS,JS) is multiplied by FI and 
FJ before its contribution is added to Z 

ij 

Subroutine GGMM calculates the impedances Q(KK,LL) which 
are like the P(IS,JS) but have different conventions for 
reference directions and subscript meaning. The 
transformation from the Q impedances to the P impedances is 
accomplished in the DO LOOP ending with statement 13. 

If the wire has finite conductivity , the appropriate 
modification is applied to the impedance matrix just above 
statement 15. The terms arising from the dielectric shell 
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on an insulated segment are obtained from subroutine DSHELL 
just above statement 16. Finally, the lumped loads, ZLD, are 
added to the diagonal elements of the impedance matrix in 
the DO LOOP ending at statement 23. 

K is a segment of test dipole I, and L is a segment of 
expansion mode J. When the segment numbers K and L are 
egual, SGANT calls GGMM to obtain the mutual impedance 
between two filamentary electric monopoles. These monopoles 
are parallel and have the same length. Monopole K is 
positioned on the axis of the wire segment, and monopole L 
is on the surface of the same wire segment. Thus, the 
displacement is egual to the wire radius. The two monopoles 
are side-by-side with no stagger. 

When segments K and L intersect, SGANT again calls GGMM 
for the mutual impedance between the two filamentary 
monopoles. Monopole K is situated on the axis of wire 
segment K, and monopole L is on the surface of wire segment 
L. The axes of segments K and L define a plane P,and 
monopole K lies in this plane. Monopole L is parallel with 
plane P and is displaced from it by a distance egual to the 
wire radius. 

CALLED BY: MAIN 

CALLS TO: CBES 

DSHELL 

GGMM 

GGS 
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SUBROUTINE SGANT (I A , I 6 , I NM . I^T , I SC . I 1 , I 2 * I 3, JA . J8 ,M0, N, NO ,NM .NP , A 
1M ,BM,C,CGO,CMM,O.EP2, E P J , E T A , FHZ , GAM , $G0 ,X , Y , Z , ZLD , ZS, ERR, IGRO) 

C CM PL EX ERR 

COMPLEX ZG,ZH» Z S « EGO,GO,CGDS, SGDS.SGOT *801 

COMPLEX Pit .Pl2*P2i,P22 f 01I ,Q l 2 , Q2 1 , Q22 » EP2 , E P ,E TA , GAM , E P3 

Complex epsi la.Cwea.betA. zarg 

COMPLEX P(2, 2), Q(2,Z),CGD( l » , SGD( I ) .C ( II • ZLOIII 
0 I MENS I ON XII). veil. zm r 0(1). I A ( 1 ) i 18(1), MO ( l NM * 4 ) 

DIMENSION IMlI, 12(1). 1 3 1 l ) • JA(I), J&(1), N0(1), 1SC(I) 

OAT A EO, TP, U0/B.B54E- 12,6.2831 8, 1 .2566E-6/ 

EP - EP3 

ICC - (N*N»N)/2 

00 i I-I.ICC 
L C ( I ) = ( .0,.0) 

ZS * (.0,.0) 

IF (CMM.LE.D.) GO TO 2 
OMEGA = TP*FHZ 

EPSILA - CMPLX(EO,-CMMM.E6/OHEGA) 

CWEA - ( .0,1 . )*OMEGA*EPSILA 

BETA • OMEGA*SQRT(UO) ♦CSQRTIEPSILA-EP) 

ZARG ■ SETA* AM 
CALL CBES (ZARG, 801) 

ZS - BfTA*BOI/C*EA 
l ZH - ZS/( T P* AM*G AM I 
OMIN - 1 • 1 30 
OMAX - .0 

00 3 J » 1 • NM 
K • I A ( J I 
L - I 8 ( J ) 

D(J) - SQRT( (X(K)-X(L) ) ♦♦ 2* ( Y( K )— Y( L ) ) *♦ 2 ♦( ZCK)-Z(L) •♦♦21 
IF (O(J).LT.OMIN) OMIN-O(J) 

IF (0( J).GT.OMAX) OMAX-O(J) 

EGO * C E X P ( GAM*0 ( J ) ) 

CGD(J) * ( EG0*i./EG0)/2. 

5 SGD(j) - ( EGO- I • / EGO )/ 2 • 

IF (OMIN. LT. 2. ♦AM) GO TO 4 
IF (CA8S( GAM*AM) .GT.0.06) GO TO 4 
IF ( C A8S( GAM* OMAX ).GT.3.J GO TO 4 
IF (AM.GT.O.) 50 TO 5 
t CONTINUE 
N=0 

WRITE (6.24) AM, OMAX, OMIN 
WRITE (6,25) 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
IB 

19 

20 
21 
22 

23 

24 

25 

26 
27 

IS 

30 

31 

ll 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 



00 19 K-I.NM 

1 FLAG * 0 

IF (( I5R0.GT.0) . AN0.IK.GT.NM/2) ) IFLAG* 1 



NDK 
KA - 
KB - 
OK * 
CGOS 
SGOS 



NO ( K > 

I A ( K ) 

I B ( K ) 

D( K) 

* CGO(K) 
- SGO(K) 



00 19 L-i.NM 
JFLAG = 0 

IF ( ( IGRO.GT.O) .AN0.( L.GT.NM/2) ) JFLAG-l 
NOL * NO ( L ) 

LA * I A ( L ) 

L 8 • I B ( L ) 

OL * 0(1) 

SGOT - SGO(L) 

NIL - 0 



00 19 

1 



I ( = 



I , NDK 
M0( K. I I ) 

= ( I-l)*N-( IM-I )/2 
= l . 

(KB. EO. I2( I ) ) GO TO 6 
(KS.EQ.IHI)) FI — I. 

1 



1 3 ( m fi — i. 



MM 
FI 
IF 
IF 
IS 

GO TO 7 
IF ( KA. EQ 
IS ■ 2 



00 19 J J= I »NDL 
J = MO ( L, JJ) 

MMM * MM* J 
IF (I .G T . J ) GO TO 
FJ - I. 

IF (L8. EQ.I2( J) ) 

IF (L8. EQ. IK J) ) 

JS * I 
GO TO 9 

IF (LA.EQ.I3I J) ) 

JS • 2 
IF (NIL.NE.O) GO TO 18 
NIL - I 

IF (K.EQ.L) GO TO 14 

INO - (LA~KA)*IL8-KA)*(LA-KB)*(te-KB) 
NGRO • IGRO 

IF (IFLAG.EQ. JFLAG) IGRO — I 



19 

GO TO 8 
FJ— i. 



FJ— l. 
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50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 
7B 

79 

80 
BI 
82 

83 

84 

85 

86 

87 

88 
B9 

90 

91 

92 

93 

94 

95 

96 
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IF (1ND.EQ.0) __ 
SEGMENTS K AND L 



CALL GGS 
I.Y(L“ 



GO TO 10 

SHARE NO POINTS 



1 LB) * Y (LB) 
2Pi2.ll .PI 



igrO - ngr6 

GO - - 






TO 18 

SEGMENTS K AND L SHARE ONE POINT (THEY I NTERSECT ) 

10 KG « 0 
JM * KB 
JC » KA 
KF - 1 

INO * <KB-LA)*(K8-LB) 

IF ( INO. NE *0 1 GO TO 11 
JC - KB 
KF - -1 
JM - KA 
KG - 3 

11 LG • 3 
JP « LA 
LF • -l 

IF IL8.EQ.JC) GO TO 12 
JP - LB 
LF - 1 
LG * 0 

12 £GN - KF*LF 



(IX(JP)-X(JC))*lX(JMl-XI JCJ )HY( JP)-Y(JC) )*(V(JN)-V(JCI'IMZ 

n jp)-z_(jci>*izjjM)riiJC|ii/i 



:psi 

_ ! JP)- 

CALL 



I JC1 

GGMH ( • O.DK . • O.DL * Al 

1 ) *Q( 2. 1 ) .01 2 . 2 ) ) 



m;c£os?§&d$ » SGDT .CPS(.ETA.GAM.Q(1.1).Q(1.2 



00 13 KK*i *2 
KP ■ 1ABSIKK-KG) 



00 13 LL-1.2 
LP * I ABSI LL-LG ) 

PC KP • L P) - SGM-QUK.LL) 
13 CONTINUE 



I GRD*NGRD 
GO TO ‘ 

K- 






t TO 10 
L (SELF 

|: 1:8;: 



REACTION OF SEGMENT K) 



■■ l tVI tvl 

(CNM.Lt *0 •) GO TO 15 
GO * GAM*DK 
ZG - ZH/I SGD$**2 j 
Oil - ZG*($GD$*CGD 



S0S-G01/2. 



012 ■ ZG*(GD*CGD$-SGDS)/2. 



15 ISCK - I SC (K | 
Pll • I .0 .ol 
P 12 ■ I .0*.G) 



"U - I .U« .Ol 

IF (ISCK.EQ.O) GO TO 16 



IF IBM.LE.AM) GO TO 16 
CALL 1 " “ “ 



OSHELL i AM« BM. OK .CGOS » SGDS . EP2. EP. E T A* GAN* PI 1 * P12 1 
16 Oil - PlUQll 
012 - P12*Q12 

CALL GGMM ( .0,0K , .0. OK .AM, CGOS . SGDS.SGOS , 1. .ETA.GAM, PI l # P12, P21 ,P2 

Oil » P 1 1 *01 l 
P 12+012 



17 



P I 2*2 1 - Oil 
IF (KA.NE.LA) GO TO 17 
GO TO IB 
Pll.l) • -012 
Pll. 2) » -Oil 
PI2.ll - -Oil 



PI 2.2) - -012 
IB CtMMM) - C(MMi)*Fl*FJ*PC(S»J5) 
19 CONTINUE 



00 23 I -1, N 

MM • I I ~ 1 1 *N-( I*I-I)/2 
IJ * MM*I 
JJA ■ JA(l) 

J1 - JJA 

Ilf I!!!! 

IF (II2.EQ.IBU1 )) Jl»Jl*NN 
JJB - JBII) 

n uit •EQ.I B( J2 ) ) J2« J2*NM 

CIIJ) • Cl IJ)*ZLD< JlFfZLDI J2) 
JJJ - JJA 



po » K 5 }.? 



NOJ 



JJJ) 



00 21 J J* 1 .NOJ 

. .. 

I Jl «NE .112) GO TO 21 



, MOi J JJ, J _ _ 

|F U.EQ.lJ GO TO 21 

If ( 121 . * 
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98 

99 
100 
loi 

.81 

104 

05 

06 

07 

08 

09 

10 



3 

.4 



17 

B 



9 

20 

21 

ll 

24 

25 

26 
27 
2B 

29 

30 
y 3l 
1 32 

33 



It 



36 

37 



18 



4 

SI 

43 

44 



45 

46 

47 
4B 

149 

. 1 ? 

1! 

55 

156 

57 

18 

60 



*64 

65 

66 
67 
6B 

ft 

71 

72 

73 

74 

75 

76 

77 
7B 
75 
80 



ll 



B3 

84 

85 



18 

90 

‘SI 
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ij « Mrttj 
FI - 1. 

IF (K.EQ.2) GO TO 20 
IF Ul( JI.NE.IIII FI — 
cuji « C< IJ)»FI*ZL0U1> 

GO TO 21 

20 if mui.NE.m in fi— i. 

CUJ) - C< IJUF I*ZLD< J2) 

21 CONTINUE 

22 JJJ • JJB 

23 CONTINUE 
RETURN 

24 FORMAT !3Xt 9 AM - # lE 10.3.3X t?0MAX • •,£ 10, 3,3X,» OMIN • ’.E10»3> 

25 FORMAT I 9 WARNING A******** ♦♦•♦•••••♦♦♦♦•A********************** • / 
^AR^fc^CORkiM | CEE0 LI " IT 0F THIN WIRE CONOITIONt THE RESULTS 

ENO 
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SORT 



PURPOSE: to define the set of dipole modes. 

o 

METHOD: In the DO LOOP ending at statement 3, the set 

of dipoles is defined by filling the vectors 11(1) and 13 (I) 
(the endpoints of dipole I) ; 12 (I) ( the terminal point of 

dipole I); and the vectors JA(I) and JB(I) (the monopoles 
comprising dipole I) with the node numbers and segment 
numbers, respectively. The DO LOOP ending at statement 8 
determines MD(J,K) (the list of dipoles sharing segment J) 
and ND (K) (the number of dipoles sharing segment J) . 

CALLED BY: MAIN 

CALLS TO: NONE 
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SU8R0UT INE SORT ( IA«IB»I1»I2» 1 3 * J A* JB* MD • ND t NM« NP» N» MAX « M IN » 1CJ» IN 
1 M I 

DIMENSION JSP(20) 

D I ME NS I ON 11(1), 12(1), 13(1), JA(|i, JB(1) 

DIMENSION I A ( 1 1 « I 8 ( 1 » , ND(1), M0(INM,4) 

I * 0 

00 3 K = I » NP 
NJK * 0 

DO 1 J * 1 » NM 

IND * ( IACJ)-K)*( 18(J)-K) 

IF < IND.NE.O) GO TO 1 
NJK = NJKM 
JSP(NJK) ■ J 

1 CONTINUE 

MOD = NJK- I 
IF CMOD.LE.O) GO TO 3 

DO 2 I MD* 1 » HOD 

1 = IM 

IF ( I.GT.ICJ) GO TO 2 
I PD * IMDM 
J A I * JSPC 1M0) 

JA(I) = J A I 
J8I * JSP(IPO) 

jam * j8i 
inn - i a t j a 1 1 

IF I I A ( JA I ) .EO.K) IMI)-IB(JAl) 

12(1) « K 
13(1) - IA(JBI) 

IF ( IA( JBD.EO.K) 1 3( 1 ) a l B( JB1 ) 

2 CONTINUE 



CONTINUE 
N - I 

DO 4 J-l.NH 
ND( J ) ■ 6 

DO 4 K* 1*4 
MD(J*K) • 0 

III - N 
IF (N.GT.ICJ) 



111 



ICJ 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 



DO 8 1*1,111 
J * J A ( I i 



DO 7 L* 1 > 2 
ND(J) * ND( J ) ♦! 

K * 1 
M * 0 

MJK * MO(JtK) ^ , 

IF (MJK.NE.O) GO TO 6 
M * 1 

MD(JtK) * I 
K * KM 



IF (K.GT.4) 
IF (H.EQ.O) 

7 J - J8(l) 

8 CONTINUE 



GO TO 7 
GO TO 5 



MIN 

MAX 



100 

0 



DO 9 J* I , NM 
NDJ * NDU) 

IF (NDJ «GT .MAX ) MAX-NOJ 
IF (NDJ.tT.HIN) MIN-NDJ 

RETURN 

END 
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SQfiOT 



PURPOSE: to solve the set of simultaneous equations to 

determine the currents on the thin wire structure. 

METHOD: This subroutine considers the matrix equation 

ZI = V which represents a system of simultaneous linear 
equations. NEQ denotes the number of simultaneous equations 
and the size of the matrix Z. 

On entry to SQROT, S is the excitation column V. On 
exit, the solution I is stored in S. Z(I,J) denotes the 
symmetric square matrix. Also on entry, the upper-right 
triangular position of Z(I,J) is stored by rows in C(K) with 

K = (I - 1) * NEQ - (1*1) / 2 + J . 

If 112 = 1, SQROT will transform the symmetric matrix 
into the auxiliary matrix (implicit inverse) , store the 
result in C (K) and use the auxiliary matrix to solve the 
simultaneous equations. If 112 = 2, this indicates that C(K) 
already contains the auxiliary matrix. 

The transfroma tion from the symmetric matrix to the 
auxiliary matrix is accomplished in the DO LOOP ending at 
statement 5. The solution of the simultaneous equations is 
accomplished in the remainder of the program. 

CALLED BY: GFFLD 

CALLS TO: NONE 
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on 



SUBROUTINE SQROT ( C , S , I WR , I 12, NEQ) 
COMPLEX C< l),S(l),S$ 

N - NEO 

IF ( 112. EQ. 2) GO TO 6 
CUJ - CSQRUC(U) 

00 l K-2.N 
1 C(K) - C<KI/CU I 



00 5 1*2. N 
IMO * I-l 
IPO - !♦! 

ID * < I -1 1 *N-( I * I -I 1/2 
II • I 0» I 

00 2 1*1, IMO 

LI * <L-l)*N-(L*L-L)/2M 

cun * cmi-ciu )*ccli i 

cun * CSQRUCC III) 

IF ( IPO.GT.N) GO TO 5 

00 4 J* I PC* N 
IJ - IOM 

00 3 MM, IMO 
MO - (M-l ) *N-I M* M-H ) /2 
MI * MOM 
MJ - MOM 

CUJI * CUJI-C1HJ )*C(MI I 



CUJ) - CUJ)/C( I I ) 
CONTINUE 

sin - sm/cm 



00 0 
IMO - 



I *2,N 
I-l 



00 7 L-l.IMO 

LI - (L-1)*N-(L*L-L)/2M 
SU) * SIIK(LII*S(ll 



I I - 
SU) 



11-11 ♦N-I 1*1 — I ) /2M 
- S(II/C(III 



1 

2 

3 

4 

5 

6 
7 
0 
9 

10 

11 

\l 

14 

15 

I; 

10 

19 

20 
21 
22 

23 

24 

25 

26 
27 
20 

29 

30 

31 

li 

34 

35 

36 

37 
30 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 



9 

10 

11 



12 



13 



14 

15 



NN * ( ( NM ) *N ) /2 
SIN) ■ S(N)/C(NN) 

NMO * N-l 

00 10 1*1, NMO 
K = N-I 
KPO * K ♦ I 

KO ® CK-l )*N-(K*K-K)/2 

00 9 L*KP0,N 
KL = KD*l 

SI Kl - S ( K ) -C ( KL ) • SC L ) 

KK * K0»K 
S I K ) - S( K ) /C( KK ) 

IF ( IMR.LE.0I GO TO 13 
CNOR * .0 

00 II I M.N 
SA * CABS(SCI)) 

IF (SA.GT.CNOR) CNOR-SA 

IF (CNOR^LE.O.) CNOR-1. 

00 12 I * 1 * N 
SS • S< I) 

SA * CABS ( SS ) 

SNOR ■ SA/CNOR 
PH * .0 

IF (SA.GT.O.) PH - 57. 29570* AT AN2 (AIM AG C SS) «REAL( SS) ) 
WRITE (6,14) I ,SNOR,SA,PH,SS 

WRITE (6,15) 

RETURN 

FORMAT ( 1 X • I I5,IF10.3,1F15.7, 1F10.0.2F15.6) 

FORMAT ( 1 HO ) 

ENO 
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.READ • I 



. • BLNK 
.EQUAL 

{ • « A 

• • * * • 'LEFT 

* * NUMBER 



. SORT 





. *CBES 
DS HELL 

. .SCANT* 1$: 

. • • '*.*•• GGS • • < 

'"© 

GANT1 @ 

* # . 

* GNFLD 

* . # 

GFFLD* •::**** G II 
'* © 










POLPRT* 



•POLPLOT* 

•SART 










LINECK 

NUMB 



. . • .EXP J 
• • 

GGMM 

. . GDISS 

• | ••• • .RITE 
• • . 

* SQROT 

CALLING SEQUENCE OF THE SUBROUTINES 
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SYMBOL DICTIONARY 



A 

ABAP 

ABAT 

ABIP 

ABIT 

ACSP 

ACST 

AFFP 

AFFT 

AM 

BM 

C 

CG 

CGD 

CJ 

CMM 

D 

ECSP 

ECST 

EFF 

EP 

EPP 

EPPS 

EPTS 

EP2 

EP3 

EP4 

ERR 

ER2 



characters of the input data cards 

backscattering phi plane angle for plotting 

backscattering theta plane angle for plotting 

bistatic scattering phi plane angle for 
plotting 

bistatic scattering theta plane angle for 
plotting 

absorption cross section for phi polarization 

absorption cross section for theta 
polarization 

far-zone phi plane angle for plotting 

far-zone theta plane angle for plotting 

radius of the thin wire of the structure 

outer radius of the dielectric shell of the 
insulation of the wire 

elements of the open-circuit impedance matrix 

branch currents for the structure 

cosh yd for a given segment 

loop currents for the structure 

conductivity of the wire 

length of a given segment 

extinction cross section for phi polarization 

extinction cross section for theta 
polarization 

radiation efficiency 

loop currents induced by a phi polarized wave 

phi-polarized far-zone feild of the dipole 
mode 

scattered electric field in the phi direction 
due to a phi polarized wave 

scattered electric field in the theta 
direction due to a phi polarized wave 

complex permittivity of insulation 

complex permittivity of ambient medium 

complex permittivity of ground 

EP4/EP3 

relative dielectric constant of insulation 
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ER3 

ER4 

ET 

ETA 

ETPS 

ETT 

ETTS 

EX 

EY 

EZ 

EO 

FHZ 

FMC 

GAM 

GG 

GPP 

GTT 

HGT 

IA 

IB 

IBISC 

ICABD 

ICJ 

IFLAG 

IGAIN 

IGRD 

INC 

INEAR 

INM 

INT 



relative dielectric constant of the ambient 
medium 

relative dielectric constant of the ground 

loop current induced by a theta polarized 
wave 

intrinic impedance of ambient medium 

scattered electric field in the phi direction 
due to a theta polarized wave 

theta polarized far-zone field of the dipole 
mode 

scattered electric field on the theta 
direction due to theta polarized wave 

near-zone electric field in x direction 

near-zone electric field in the y direction 

near-zone electric field in z direction 

8.8 54E-12 

frequency in hertz 
frequency in megahertz 

intrinic progration constant of the ambient 
medium 

time-average power input 

power gain associated with the phi polarized 
component 

power gain associated with the theta 
polarized component 

height of the structure above ground plane 

first node of a given segment 

second node of a given segment 

indicator for bistatic scatter calculations 

indicator for the data cards 

dimension corresponding to the number of 
simultaneous linear equations 

indicator for program termination 

indicator for antenna gain calculations 

indicator for presence of the ground plane 

indicator for the type of far-zone 
calculations 

indicator for near-zone calculations 

dimension corresponding to the number of 
mon poles 

number of integration steps 
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ISC 

ISCAT 

IWR 

11 
1.12 

12 
13 
JA 
JB 

KFLAG 

KGEN 

LOAD 

LZD 

MAX 

MD 

MIN 

MSG 

N 

ND 

NGEN 

NM 

NPL 

OMEGA 

PH 

SCSP 

SCST 

SGD 

SIG2 

SIG3 

SIG4 

SPPM 

SPTM 

STPM 

STTM 



indicator for the insualtion 
indicator for backscatter calculations 
indicator for current distribution output 
endpoint node of a given dipole 
indicator for auxiliary matrix 

terminal node number of a given dipole 

endpoint node number of a given dipole 

first segment number of a given dipoile 

second segment number of a given dipole 
print indicator 

list of generator/feed locations 

indicator for structure load 

list of impedance/load locations 

maximum of the number of segments connected 
to any one given node 

list of dipoles sharing a given segment 

minimum of the number of segments that 
connected to any one given node 

indicator for error printout 

number of simultaneous linear equations 

total number of dipoles sharing a given 
segment 

indicator for antenna calculations 
number of segments 
indicator for polar plot 
angular frequency 

phi angle for far-zone calculations 

scattering cross section for phi polarization 

scattering cross section for theta 
polarizatxon 

sinh /d of a given segment 
conductivity of insulation 
conductivity of the ambient medium 
conductivity of ground 

echo area phi incident-phi scattered wave 
echo area phi incident-theta scattered wave 
echo area theta incident-phi scattered wave 
echo area theta incident-theta scattered wave 
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TD2 loss tangent of the insulation 

TD3 loss tangent of ambient medium 

TH theta angle for far-zone calculations 

TP 2n (6.28318) 

UO 1 . 2566E-6 

VG antenna complex driving voltages 

VOLT list of VG's 

X x-coordinate of each node 

XNP list of XP's 

XP x-coordinate for near-zone calculations 

Y y-coordinate of each node 

YNP list of YP's 

YP y-coordinate for near-zone calculations 

Y11 complex power input 

Z z-coordinate of each node 

ZLD complex load at a given node 

ZLLD list of ZLD 1 s 

ZNP list of ZP's 

ZP z-coordinate for near-zone calculations 

ZS surface impedance of the wire 

Z11 antenna input impedance 
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APPENDIX B 
USER'S MANUAL 



99 



TABLE OF CONTENTS 



PROGRAM LIMITS 101 

MINIMUM DATA 101 

OUTPUTS 102 

DATA CARDS 102 

1* WIRE 104 

2* INSULATION 105 

3* EXTERIOR MEDIUM 106 

4* DESCRIPTION 107 

5* GEOMETRY 108 

6 • FEED 109 

1 * LOAD 110 

8 • OUTPUT 111 

9 • PLOT 113 

10* GROUND 114 

11* INTERVAL FOR CALCULATION 116 

12* GENERATOR 117 

13* IMPEDANCE 118 

14* FREQUENCY 118 

15* CHANGE 119 

16« END 119 

17* STOP 119 

SAMPLE PROBLEM * 120 



100 



USER'S MANUAL 



The Antennas-Scatterers Analysis Program (ASAP) for thin 
wire structures in a homogenous conducting medium performs a 
frequency domain analysis of antennas and scatters. The 
program is applicable in the presence of either a perfect or 
a finite ground. This appendix will describe and explain 
the data cards necessary to execute the compute program. 
Although the program was written for the IBM 360 computer 
system, it can be executed on another system with minor 
modifications. 

The program utilizes piecewise sinusoidal expansion for 
the current distribution with Kirchhoff Current Law enforced 
everywhere on the structure. If the structure contains end 
points, the currents at these points are assumed to vanish. 

I. Program Limits 

The thin wire assumptions are questionable and the 
accuracy and convergence deteriorate if the radius of wire 
utilized for the structure exceeds 0.01 of a wavelength, if 
the longest segment is greater than one-fourth of a 
wavelength, if the length ratio of the longest and shortest 
segments exceeds 100, or if the total wire length is less 
than 30 times the wire diameter. If a wire is bent sharply 
to form a small acute angle (less than 30 degrees) , the thin 
wire model is questionable. It is assumed that the wire 
conductivity greatly exceeds the conductivity of the ambient 
medium, for insulated wires, the dielectric layer is assumed 
to be electrically thin. 

II. Minimum Data 

The minimum data necessary to execute the program is: 
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a. description of structure 

b. radius of wire used for the structure 

The program will default to the other parameters necessary. 
The default parameters are: 

a. wire for the structure is copper 

b. frequency of operation is 300 mhz 

c. homogneous medium is free space 

A more detailed explanation of the defaults will be 
discussed when the data card for the parameter is described. 

III. Outputs 

In antenna problems, the output includes structure 
currents, impedance(s) of feed (s) , gain, polar radiation 
plots, and near field calculations. In bistatic scattering 
problems, the output includes structure currents, complex 
elements of the polarization scattering matrix, polar 
reradiation pattern plots, and echo areas produced by a 
plane wave. For backscatter ing problems the output includes 
absorption, scattering and extinction cross sections in 
addition to the outputs of bistatic scattering. Most of the 
outputs are suppressed and must be requested. Since the 
program can produce a large volume of output, care should be 
exercised until the user is familiar with the outputs. 

IV. Data Cards 

The Analysis Program utilizes free format for the data 
cards, that is, the program utilizes character recognition 
to determine which parameters are being read. Data placement 
(location) on the input card is not critial. Blank 
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characters, on all input cards but the COMMENT data card, 
are ignored and may be used at the discretion of the user. 
Since character recognition is used, only the first four 
characters of the key words must be present and correct. 

The format for the COMMENT CARD utilizes standard 
FORTRAN format (i.e. 'C' in column 1 followed by at least 

four blanks) . The COMMENT CARD is the only type of input 
card that position in the data block is critial. This 
(these) card (s) must be placed at the beginning of a data 
block. A data block is a series of related data cards. 
Several data blocks may be used to define an analysis 
problem. This will become clear when the termination cards 
(END, STOP, or CHANGE) are discussed. There is no limit to 
the number of comment cards that may be used. As a check for 
the user, all input data cards will appear on the output as 
they appear in the input deck. 

The format of other data can be of one of two forms: 

a. type of card (option 1/option 2/ ) 

b. parameter (value) . 

The type of format to use will be apparent as the individual 
data cards are discussed. 

The numerical values for the parameters may be stated in 
any one of the following forms. The program will translate 
the number to the proper form for the specified parameter, 
either fixed or floating point. All of the following 
examples have the same value. 

0.0001 or .0001 or 100. U or 100U or .1M or 0.1M or .0000001K 

-6 -3 3 

U = 10 M = 10 K = 10 
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1. WIRE This card is used to define the parameters 
associated with the wire utilized by the thin wire 
structure. Two options are available and are defined as: 

RADIUS=value of the radius of the wire in meters 

CO NDUCTIVITY= value in megamhos per meters . 

The wire data card must appear in the first data block to 
define wire radius. The default value of the conductivity 
is 50 megamhos/meter (copper) . 



'WIRE ( RADIUS= . 00 1/ CONDU=28.5) 



2. INSU L A TION This card is utilized to define the parameters 
associated with the insulation of the wire used for the 
structure to be analyzed. If this card is omitted, the 
program assumes that the structure is uninsulated. Four 
options are available and are defined as: 

RADIUS=value of outer radius in meters 

CONDUCTIVITY=value in micromhos per meter 

DIELECTRIC=value of relative dielectric constant 

LOSS TANGE NT=value . 

The conductivity and either the relative dielectric constant 
or the loss tangent (but not all three) options may be 
stated. 



INSULATION ( RADIUS=. 0 15/ C0ND=7 . /DIEL=5) 



3. EXTERIOR MEDIUM This card is utilized to describe the 
homogeneous medium surrounding the structure. If the medium 
is free space, this card may be omitted. Three options are 
available and are defined as: 

DIELECTRIC=v alue of relative dielectric constant 

CONDUCTIVITY=value in micromhos per meter 

LOSS TANGE NT=value . 

As with INSULATION card state either conductivity or loss 
tangent. 



/EXTE (LOSS=. 4 5) 
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4. DESCR IPTI ON This card is utilized to describe the shape 
of the wire structure to the program. The user must divide 
the wire structure into segments of the appropriate length 
and number each node starting at one. A node is a point 
where a segment begins or ends. A maximium of four segments 
can meet at any given node. An isolated wire must contain 
at least two segments and three nodes. Thus the DESCRIPTION 
CARD must show at least 3 consecutive nodes for all portions 
of the wire structure. The structure is described by 
stating the node numbers that each segment connects. The 
description of a square loop might appear as; 

DESCRIPTION ( 1-2/2-3/3-4/4-1) . 

The description of a dipole and reflector might appear as; 

DESCRIPTION (1-2/2-3/3-4/4-5/6-7/7-8/8-9/9-10) . 

If the description will not fit on one data card continue on 
the next card as if the previous card were longer. The 
dipole example might appear as; 

DESCRIPTION ( 1-2/2- 3/3 -4/4- 5/ 

6-7/7-8/8-9/9-10) . 

Note that the last character on the card to be continued is 
a slant (/) . As many cards as necessary may be used. The 
maximum number of nodes permitted is fifty. If ground plane 
is present, the maximum number is twenty-five. If a ground 
plane is present and the structure touches the ground plane, 
the lowest node numbers MUST be used for the touching nodes. 
That is, if the structure touches the ground plane at two 
points, node numbers 1 and 2 MUST be assigned to these 
nodes. 



DESCR (1-2/2 -3/3 -4/4-1) 
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5. GEOMETRY This card is used to state the physical location 
in rectangular coordinates of each node of the DESCRIPTION 
CARD . The rectangular grid is in units of meters. If node 1 
is located at x1,y1,z1 and node 2 at x2,y2,z2 and node 3 at 
x3,y3,z3,etc. , the GEOMETRY CARD might appear as: 

GEOMETRY (xl , yl r z1/x2, y2,z2/x3,y3,z3/ ) 

As with the DESCRIPTION CARD, continuation cards are 
permitted. 



'GEOM (. 1 ,0, 



1 /-. 1 , 0 . 1 /-. 1 , 0 -. 1 /. 1 , 0 ,-. 1 ) 
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6. FEED For antenna analysis the feed point (s) and 
voltage (s) must be stated. In the foremention dipole and 
reflector example if the feeds were at node 2 with a voltage 
source of .5 at an angle of -90 degrees and at node 4 with a 
voltage source of .5 at an angle of +90 degrees the FEED 
CARD might appear as: 

FEED ( 2,. 5 # -90/4, .5, +90) 

The order of the information for each voltage source is node 
number, magnitude, and phase angle. This order is repeated 
until all sources are stated. If the source information will 
not fit on one card, use another card similiar to the 
initial one; that is, repeat the word '•FEED”. If only one 
voltage source is applied to the structure, only the node 
number must be stated. In the dipole example, if the drive 
is at node 3, the FEED CARD might appear as: 

FEED (3) 

A default source of one volt at zero degree phase is 
assumed. Voltage sources should only be stated for nodes 
with only two segments. 



FEED (2,. 5, -9 0/4,. 5, +90) 



7. LO AD This card is used to describe the loads to be placed 
at various locations on the structure. The format for this 
card is similiar to that of the FEED CARD, that is, the 
word "LOAD" is used in the place of "FEED” . The order of 
the information on the card is the same. Since this card is 
frequency dependent, it must be changed if the frequency of 
operation is changed. No default parameters are available. 
The structure is assumed unloaded unless this card is used. 
Once the structure is loaded, it will remain loaded for the 
remainder of the data block series. To unload the structure 
the following card may be used: 

LOAD (- 1 ) 
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8. OUTPU T This card is used to request output data. Most of 
the output is in tabular form. More than one OUTPUT CARD 
is permitted per data block, but not for the same type of 
output. If only the antenna input impedance, antenna 
efficiency, or time-average power input is of interest, no 
OUTPUT CARD is necessary. These parameters are automatically 
printed if a FEED CARD or GENERATOR CARD is utilized. One 
or more of the following options may be used to request the 
various outputs available. 

FAR FIELD=phi initial, phi final, theta initial, theta 

final 

This option gives the components of the electric field 
intensity in the far field as phi and theta varies between 
limits specified in one degree divisions. 

BA CKSC ATERING=phi initial, phi final, theta initial, 
theta final 

This option gives the absorption, scattering, and extinction 
cross sections, and the complex elements of the polarization 
scattering matrix for an incident plane wave illuminating 
the structure from the spherical direction of phi, theta as 
both vary between limits specified in one degree divisions. 

BISTATIC=phi inital, phi final, theta initial, theta 

final 

This option gives echo area and the complex elements of the 
polarization scattering matrix for an incident plane wave 
illuminating the structure from the spherical direction phi, 
theta final of the backseat tering output option, reradiated 
in the phi, theta direction as both vary between limits 
specified in one degree divisions. A bistatic output request 
must be accompanied with a backscattering request in the 
same data block. 

STEP=value in degrees 
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This option will cause any of the above output options to be 
stepped at a different interval size. That is, if one of the 
above options is to be stepped at ten degrees intervals, use 
this option. This option overrides the one degree stepping. 

CURRENT 

This option gives the currents on the structure which are 
produced by the feed/generator voltages and/or the incident 
plane wave of the backscattering request. 

NEAR=x1 ,y1,z1 

or 

NEAR= (Xl,y1,z1/x2,y2,z2/x3,y3,z3/etc ) 

This option gives the value of electric field components in 
the near field for the antenna at the point or points 
specified. 



/OUTPUT (FARF=45,50,25, 50) 



% 
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9* PLOT This card will produce normalized polar plots in the 
specifed plane for the stated option. The plane is 

specified by stating either "PHI= " or “THETA 3 The 

PLOT CARD overrides the limits of the OUTPUT CARD for the 
same option. If only a normalized pattern is of interest, 
only a PLOT CARD is necessary. If a table of values and a 
normalized pattern is desired, both a PLOT CARD and OUTPUT 
CARD must be used. Only one PLOT CARD is permitted per 
data block. The following pattern plots are available: 

FAR FIELD/plane 

This option will plot the far field intensity for each 
component of the electric field. 

BACKSCATTERIN G/plane 

This option will plot the normalized magnitude of each of 
the elements of the polarization scattering matrix. 

BIS T ATI C/ plane 

This option will plot the normalized magnitude of each of 
the elements of the polarization scattering matrix produced 
by the incident plane wave stated by final limits of the 
backscatter ing option of the output request. 



'PLOT (FARF/THET=90) 
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10. GROU ND This card is used to describe the ground 
parameters if a ground plane is present. If no ground plane 
is present, the structure is assumed to be in free space or 
the homogeneous medium of the EXTERIOR MEDIUM data card. 
Seven options are available and are defined as: 

PERFECT 

This option will analyze the structure over a perfect ground 
plane. 

GOOD 

This option will analyze the structure over a good ground 
plane where the conductivity of the ground is .02 mhos/meter 
and the relative dielectric constant is 30. 

POOR 

This option will analyze the structure over a poor ground 
plane where the conductivity of the ground is .001 
mhos/meter and the relative dielectric constant is 4. 

SEA 

This option will analyze the structure over salt water where 
the conductivty of the water is 4. mhos/meter and the 
relative dielectric constant is 80. 

HEIGHT=value in meters 

This option will analyze the structure with origin of the 
GEOMETRY card this height above the ground plane. The 
lowest point of the structure must not lie below the ground 
plane. It may lie on the ground plane. 

CONDUCTIVITY= value in mhos/meter 

This option is used to state the value of conductivity of 
the ground plane if the default values mentioned above are 
not utilized. 
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DIELECTRIC 3 value 



This option is used to state the relative dielectric 
constant of the ground plane if the default values mentioned 
above are not utilized. 



/ground (HEIG=1 O/COND 3 . 002/DI EL= 10) 



11. INTER VAL FOR C ALC ULATION This card is used to state the 
number of intervals to be used for calculating the elements 
of the impedance matrix with Simpson • s-rule integration. A 
large value for the number improves the accuracy at the 
expense of greater execution time. For most problems a 
suitable combination of speed and accuracy is obtained with 
a value of four, the default value. If the rigiorous 
closed-form impedance expressions in terms of the 
exponential integrals is desired, set this value to zero. 

INTERVAL=value 



INTE (6) 



12. GENERATO R This card is similiar to the FEED CARD in use, 
except that the segment numbers are stated instead of the 
node numbers. This is useful if three or four segments meet 
at a node. The positive terminal of the generator is 
connected to the specified segment such that current is 
forced in the the positive direction. The positive direction 
of current flow is from the first stated node number of that 
segment toward the second stated as ordered on the 
DESCRIPTION CARD. 



/GENE (2 r . 5, -90/4, 



5, +90) 
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13. I MP E DAN CE This card is similiar to the LOAD CARD in use, 
except that the segment numbers are stated instead of the 
node numbers. As with the GENERATOR CARD, this is used if 
three or four segments are connected to a node. The 
impedance will be connected to the positive terminal of the 
specified segment. The format of this card is the same as 
the LOAD CARD. 



/iMPE (1,120,-45/3, 120, +45) 



14. FR EQ UEN CY This card is used to specify the frequency in 
megahertz if it is to be other than the default value of 300 
MHz. 



/ FREQ ( 1 2. 5) 
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15. CHAN GE This card at the end of the data block signals 
the program that the following data cards are changes to the 
previously read data, for the next run. If a "CHANGE CARD" 
is used, the outputs must be requested again in the next 
data block. The "CHANGE CARD" cannot be used to change 
"DESCRIPTION CARD" or "GEOMETRY CARD" data when operating 
with a "GROUND CARD". Use an "END CARD" to make changes when 
a "GROUND CARD" is used. 



16. END This card signals the program that this is the end 
of a data block series and to reinitialize data for the next 
problem. An "END CARD" cannot be used with a "CHANGE CARD". 



17. STOP This card signals the program that all of the data 
cards have been read and to terminate itself when execution 
is completed. This card must be used as the last card in 
place of the "END CARD" of the last data block series. A 
"STOP CARD" cannot be used with an " END CARD" in the same 
data block. 
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C AN EXAMPLE PROBLEM 

C 

C V ANTENNA 

C 

WIRE (RADIUS=1 M) 

GEOMtO,-. 18, +.18/0,-. 09,*. 09/0, 0,0/0, 0.09,. 09/0,. 18,. 18) 
DESC(l-2/2-3/3-4/4-5) 

FEED! 3) 

OUTPUT (FARF=45»50»65»80/STEP=5) 

CHANGE 

OUTPUT! BIST=45, 45, 45, 45 /BACK=0 , 0 , 10 , l 2 ) 

OUTPUT ( CURRENT J 

CHANGE 

C 

C CHANGE STRUCTURE SHAPE TO DIPOLE 

C 

GEOM (0,-. 25, 0/0, -.125, 0/0,0, 0/0, . 125, 0/0,. 25,0 ) 

PLOT (FAR F/ PHI =90) 

GROUND (HEIGHT = .2 5/ GOOD ) 

STOP 



THE ABOVE DATA DECK WILL PRODUCE THE OUTPUT ON THE 
FOLLOWING PAGES. 
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FAR FIELD 
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CONTINUE EXECUTION WITH THE FOLLOWING ADDITIONS AND/OR CHANGES 
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CONTINUE EXECUTION WITH THE FOLLOWING ADDITIONS AND/OR CHANGES 
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PHI ELECTRIC FIELD ANTENNA PATTERN FOR SPECIFIED PLANE. NORMALIZING FACTOR* .28833E-06 
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